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CHAPTER  1 


INTRODUCTION 

This  document  is  the  final  report  on  the  AFOSR  supported  research,  "Design  of 
Robust  Optimal  Control  Systems  and  Stability  Analysis  of  Real  Structured  Uncertainties," 
(AFOSR  Contract  F33615-90-C-3613).  The  report  summarizes  our  research  progress  in 
the  following  problems:  design  of  robust  controllers  via  H“  optimization,  investigation  of 
the  properties  of  die  parameter-dependent  H*°  Riccati  equation,  fast  computation  of  the 
optimal  H~  norm,  design  of  robust  controllers  via  ^.-Synthesis,  controller  reduction  by 
structured  trunction,  a  parametrization  approach  to  reduced-order  H*°  controller  design, 
controller  reduction  via  observer-based  controller  parametrization,  and  dominant  modes  of 
mechanical  systems. 

H~  optimization  theory  plays  a  key  role  in  many  robust  control  problems  such  as 
tracking  and  disturbance  attenuation,  mixed  sensitivity  minimization,  model  matching,  and 
^.-synthesis  [1-10].  One  of  the  most  elegant  and  efficient  solutions  to  the  H”  optimization 
problem  is  the  state  space  approach  proposed  by  Doyle,  Glover,  Khargonekar,  and  Francis 
[7]  and  Glover  and  Doyle  [8].  As  shown  in  [7,8],  die  major  computation  involved  in  this 
approach  is  solving  two  y-dependent  FT  Riccati  equations.  The  parameter  y  serves  as  an 
upper  bound  of  the  HP  norm  of  the  closed  loop  system.  If,  for  a  given  y,  the  two  H°° 
Riccati  equations  have  positive  semidefinite  stabilizing  solutions  and  the  spectral  radius  of 
the  product  of  the  two  solutions  is  less  than  y2,  then  all  stabilizing  controllers  can  be  easily 
constructed  from  the  plant  data  and  the  Riccati  solutions.  The  dimension  of  the  resulting 
controller  is  not  higher  than  that  of  the  generalized  plant  and  the  H~  norm  of  the  closed  loop 
system  is  less  than  y. 

In  Chapter  2,  we  will  discuss  how  to  formulate  a  robust  control  problem  as  an  H~ 
optimization  problem  and  how  to  construct  a  state-space  realization  of  the  generalized  plant 
For  convenience,  Glover  and  Doyle's  suboptimal  H*°  controller  formulas  [8]  will  be  briefly 
reviewed.  These  suboptimal  H°°  controller  formulas  can  be  easily  transformed  into  a 
descriptor  form  such  that  an  optimal  H*°  controller  can  be  constructed  without  numerical 
difficulties  if  the  optimal  H~  norm  is  known.  The  optimal  H“  controllers,  with  very  few 
exceptions,  have  direct  feed  through  terms  and  therefore  infinite  bandwidth.  Hence,  control 
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engineers  may  prefer  strictly  proper  suboptimal  H~  controllers  to  the  optimal  ones. 
However,  knowing  the  optimal  H*°  norm  is  important  in  determining  which  suboptimal 
controller  to  be  chosen  in  practical  design. 

Unlike  the  well-known  LQG  Riccati  equations  [1 1],  the  H~  Rice  ad  equations  are 
parameter  dependent  and  may  not  have  solutions  for  some  values  of  the  parameter  y.  The 
H**  optimization  design  can  only  be  achieved  by  iteratively  searching  for  the  optimal  H~ 
norm,  which  is  the  smallest  y  such  that  an  H**  controller  exists.  This  search  is  the  major 
computational  burden  in  H*°  design,  and  bisection  method  has  been  widely  used  in  the 
design  procedure.  As  it  is  known,  the  bisection  method  converges  slowly  and  therefore  is 
inefficient  To  improve  the  convergence  rate,  we  need  to  develop  a  more  efficient  search 
scheme  which  may  require  a  better  understanding  of  the  H**  Riccati  equations. 

To  investigate  the  properties  of  the  two  H*°  Riccati  equations,  first  of  all  in  Chapter 
3  we  will  consider  a  special  case  with  which  the  Dn  matrix  of  the  generalized  plant  is 
assumed  zero.  With  this  assumption,  the  two  H*°  Riccati  equations  are  greatly  simplified. 
We  will  investigate  the  structure  of  the  y-domain  where  the  H“  Riccati  solutions  exist 
and/or  positive  semidefinite,  and  reveal  some  useful  properties  of  the  Riccati  solutions  on 
the  y-domain,  such  as  continuity,  monotomcity,  and  convexity  with  respect  to  y.  By  these 
properties,  quadratically  convergent  algorithms  can  be  easily  developed  to  compute  the 
optimal  H"  norm  of  the  closed-loop  system. 

The  assumption  that  the  Du  matrix  of  the  generalized  plant  is  zero  can  hardly  be 
satisfied  by  many  practical  problems.  In  the  general  case  with  Du  nonzero,  the  H~  Riccati 
Hamiltonains  are  much  more  complicated  functions  of  y  than  those  in  the  special  case  and 
hence  make  the  investigation  of  the  properties  extremely  difficult.  In  Chapter  4,  several 
decompositions  and  some  key  intermediate  variables  are  employed  to  tackle  the 
complexities  and  to  show  that  the  monotonicity  and  convexity  properties  also  hold  for  the 
general  case.  Based  on  these  properties,  a  quadratically  convergent  algorithm  can  be  easily 
developed  to  compute  the  optimal  H°°  norm.  The  algorithm  is  presented  in  Chapter  5. 

In  Chapter  6,  a  fighter  aircraft  longitudinal  flight  control  design  problem  is  used  to 
demonstrate  an  application  of  p-syn  thesis  technique.  With  the  robust  controller  obtained, 
simulations  and  analyses  are  performed  on  the  closed  loop  system.  For  the  example  we 
considered,  the  order  of  the  p-syn  thesis  controller  can  be  reduced  to  that  of  the  plant 
without  degrading  the  performance. 
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From  Chapters  7  to  9,  three  controller  reduction  approaches  are  presented.  A 
controller  reduction  approach  based  on  the  closed-loop  controllability  and  observability 
gramians  is  proposed  in  Chapter  7.  The  motivation  is  to  retain  the  input/output  relationship 
of  the  closed-loop  system  while  the  order  of  the  controller  is  reduced.  A  structured 
transformation  on  die  gramians  is  employed  to  achieve  this  objective. 

In  Chapter  8,  a  parametrization  method  is  proposed  for  reduced-order  H*°  controller 
design.  By  using  the  state-space  H“  optimization  approach,  one  can  characterize  the  set  of 
all  (sub)optimal  controllers  in  terms  of  a  parameter  matrix.  It  is  shown  that  there  always 
exists  a  suitable  parameter  matrix  such  that  a  reduced-order  H*°  controller  can  be 
constructed.  The  procedure  for  obtaining  such  a  parameter  matrix  is  also  presented  in  this 
chapter. 

In  Chapter  9,  the  reduced  compensator  design  issue  was  tackled  from  an  observer- 
based  compensator  point  of  view.  It  develops  two  properties  related  to  observer  based 
controller  parametrization  and  pole  placement  It  shows  that  the  poles  of  the  closed-loop 
system  with  the  observer-based  controller  parametrization  are  the  regulator  poles,  the 
observer  poles,  together  with  the  poles  of  the  added  stable  parameter  matrix.  If  the 
controller  is  realized  by  a  minimal  realization,  the  closed-loop  poles  will  include  all  the 
poles  of  the  added  stable  parameter  matrix  and  a  subset  of  the  regulator  and  the  observer 
poles.  It  also  points  out  that  there  exist  parameter  matrices  which  would  render  the 
controller  non-minimal,  thereby  yielding  reduced  order  controllers.  Such  order-reducing 
parameter  matrices  -  both  static  and  dynamic  -  have  been  characterized  herein.  With  such  a 
parametrization  available,  one  could  then  choose  a  K(s)  in  order  to  best  approximate  large 
order  controllers,  such  as  the  H°°  compensators,  by  a  lower  order  controller. 

Incorporating  such  a  parametrization  in  a  design  procedure  requires  the 
identification  of  poles  that  are  to  be  rendered  uncontrollable  and  unobservable.  Part  of  the 
research  has  focused  on  this  issue.  In  particular,  we  have  developed  a  means  of  identifying 
dominant  modes  of  mechanical  systems  based  on  an  L2  norm  which  is  summarized  in 
Chapter  10. 

Chapter  11  is  the  conclusion.  The  work  for  future  research  will  also  be  briefly 
described  in  Chapter  11. 
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CHAPTER  2 


FORMULATION  OF  H“  CONTROL  PROBLEMS 


2.1  The  H~  Optimization  Problem 


Most  control  problems  can  be  formulated  as  the  following  H“  optimization 
problem.  A  detailed  procedure  of  the  formulation  will  be  discussed  in  Section  2.3.  In  the 
H“  optimization  problem  formulation,  the  system  under  consideration  is  described  as 
follows. 


z(s)' 

-Gn(s) 

°12(s)l 

r*»i 

v(s)‘ 

yOO. 

V5> 

[u(s)j  °(5) 

u(s) 

(2-1) 


where  Gn(s)  e  R(s)plxml,  G12(s)  €  R(s)plxm2,  G21(s)  e  R(s)p2xml,  and  G22(s)  e 
E^y’2""2.  B(Sr*  is  the  set  of  pxq  proper  rational  matrices  with  real  coefficients.  In  (2- 
1),  z,  y,  v,  and  u  are  the  controlled  output,  the  measured  output,  the  exogenous  input,  and 
the  control  input  respectively.  The  controlled  output  vector  z  usually  includes  the  error 
signal  and  a  weighted  control  input  The  exogenous  input  v  contains  the  disturbances,  the 
noises,  and  the  commands.  The  measured  output  vector  y  consists  of  all  the  signals  which 
can  be  measured  and  available  for  feedback.  Through  the  control  input  u,  the  behavior  of 
the  system  can  be  modified.  The  vector  y  will  be  used  as  the  input  to  a  controller  K(s)  and 
the  output  of  K(s)  will  be  connected  to  the  control  input  u.  That  is, 


u(s)  =  K(s)y(s).  (2-2) 

The  standard  H°°  optimization  problem  is  the  problem  of  finding  a  proper  controller  K(s) 
such  that  the  closed-loop  system  is  internally  stable  and  II 7  j(GJQ 11^  is  minimized  where 

7  ,*(0,10  (s)  =  G„(s)  +  G12(s)  K(s)  [  I  -  G22(s)  K(s)  ]  -1  G21(s) .  (2-3) 

That  is,  7  2(G,K)  (s),  the  lower  linear  fractional  transformation,  is  the  transfer  function  of 
the  closed-loop  system  from  v  to  z. 

Let 
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'  A 

B.  V 

G(s)  = 

C! 

DU  D12 

5 

^21  ^22_ 

(2-4) 


be  a  realization  of  the  generalized  plant  G(s)  and  A  €  R0*". 


In  [8],  Glover  and  Doyle  assume  that  the  realization  of  the  generalized  plant  G(s)  is 
given  by  (2-4)  with  the  following  assumptions. 

(i)  (A,  Bj)  is  stabilizable  and  ( C 2,  A)  is  detectable. 

(ii)  rank  D12  »  m2,  rank  D21  =  p2. 


r°i  r 

(iii)  D12  *  ,  D21  *  10  IJ,  and  Dn  is  partitioned 


as 


Dllll  D1U2 
L  ^1121  D1122 


with  D1122  €  Rp2xm2, 


(iv)  D22  *  0  (this  can  be  removed,  for  details  please  refer  to  [8] ). 
joal  -  A  B2 

'  *  n  +  nij  VtoeR. 


(v)  rank 


(vi)  rank 


^1  ~i2j 
jcol  -  A  Bj 

C2  D21 J 


=  n  +  p2  VtoeR. 


Define  two  Hamiltonian  matrices  as  follows, 

H_(y)  := 


’  A 

0 " 

B 

T 

T 

_ 

T 

R'1 

DlCi  BT 

L-cjq  *a  j 

L-cXJ 

and 


Uy)  := 


at  0 


-Bib;  -a 


T 

•Bl< 


_-l 

R 


[DX  c] 


where 


and 


[Dn  d12]' 

D-  ■ 

D.r 

D 

L21 J 

(2-5a) 


(2-5b) 


(2-5c) 
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Then  the  following  theorem  shows  an  easy  way  to  construct  a  suboptimal  stabilizing 
controller  such  that  HIT  ^(G  JC)Hm  <  y  where  IT  ^(GJC)  is  the  closed-loop  transfer  matrix 

from  v  to  z. 


Theorem  2.1:  [8] 


There  exists  a  stabilizing  controller  such  that  HIT  i(G,K)ll-> 
following  three  conditions  hold. 

<  y  if  and  only  if  the 

(i)  Y>max(^[Diin  Dni2j  ®[Dnn  Dn2i  ) 

(2-6a) 

(ii)  HL(Y)  €  dom(Ric)  and  X«(y)  :=  Ric[FL(y)]  2:  0. 

(2-6b) 

(iii)  J-(Y)  e  dom(  Ric )  with  Y«(Y)  =  Ric[J„(7)]  2>  0. 

(2.6c) 

(iv)  p[X_(y)Y..(Y)]  <V2- 

(2-6d) 

Moreover,  when  these  conditions  hold,  one  such  controller  is 

Ak  Bk  * 

K  As)  = - 

ck 

(2-7a) 

where 

Dk='DU21Dllll^1'  D11UD1111^  D1112  "D1122 

(2-7b) 

Ck=  {F2-Dk(C2+F12)}Z 

(2-7c) 

Bk=-H2  +  (B2+H12)Dk 

(2*7d) 

Ak=  A+HC+(B2+H12)Ck 

(2-7e) 

Z  :=  (I-r2YX)4 

(7-7f) 

pT=[FL  f5]"(XB  +  CiDi-)R‘! 

(2-7g) 

H  =  [H„  Hj2  H2]  =  -CYCT+B,D^)R 

(2-7h) 

and  F„  €  R(mi'P2)xn,  p12  €  BP2Xn,  F2  e  Rm2*n,  Hn  e  Rnx(Pl'm2),  h12  e  Rnwn2, 
H2  €  Rnxp2. 
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In  the  above  theorem,  condition  (ii)  means  that  there  exist  positive  semi-definite 
solutions  X  and  Y  to  the  algebraic  Riccati  equations  corresponding  to  the  Hamiltonians 
H«.(y)  and  JJy)  respectively.  Condition  (iii)  means  that  the  spectral  radius  of  XY  is  less 
than  y2. 

The  above  theorem  provides  an  easy  way  to  construct  a  stabilizing  suboptimal 
controller  such  that  HIT  i(G,K)ll->  <  y.  The  order  of  the  suboptimal  controller  can  be  the 

same  as  that  of  the  plant  G(s).  The  major  computation  involved  is  the  solution  of  two  H°° 
Riccati  equations  which  are  easy  to  solve  if  solutions  exist 

Theorem  2.1  can  also  be  used  to  compute  the  optimal  H“  norm  and  to  construct  an 
optimal  H~  controller.  Algorithms  for  computing  the  optimal  H~  norm  will  be  discussed  in 
Chapter  4  and  the  construction  of  an  optimal  H~  controller  will  be  given  in  Section  2.2. 
Section  2.3  will  explain  how  to  formulate  a  robust  control  problem  as  a  standard  H~ 
optimization  problem. 

23.  Construction  of  Optimal  H~  Controllers 

The  optimum  can  occur  in  three  cases.  In  case  (1),  the  optimum  occurs  at  the 
smallest  y  such  that  the  two  H~  Riccati  equations  have  stabilizing  solutions  X  and  Y,  i.e., 
these  X  and  Y  happen  to  be  positive  semi-definite  and  p(XY)  <  y2.  Case  (2)  occurs  when 
Y=0  (or  X=0)  for  all  y  and  the  optimal  H~  norm  is  the  smallest  y  such  that  X  (or  Y)  is 
positive  semi-definite.  The  most  likely  one  to  happen  most  of  the  time  is  case  (3)  in  which 
the  optimal  FT  norm  is  the  y  such  that  the  two  H“  Riccati  equations  have  positive  semi- 
definite  stabilizing  solutions  X  and  Y  and  p(XY)  =  y2  where  p(XY)  is  the  spectral  radius  of 
XY. 

From  Theorem  2.1,  a  suboptimal  H“  controller  can  be  easily  constructed. 
However,  as  y  approaches  to  the  optimum  we  will  encounter  the  inversion  of  a  singular 
matrix  except  case  (1)  which  seldom  occurs.  To  eliminate  the  numerical  difficulty,  Safonov 
et  al.  [12]  rederived  the  optimal  controller  formulas  in  a  descriptor  form  (or  generalized 
state-space  representation). 

The  formulas  in  (2-7a)  -  (2-7h)  can  also  be  written  in  a  descriptor  form  after  slight 
rearrangement.  We  will  firstly  consider  case  (3)  which  occurs  much  more  frequently  than 
the  other  two  cases.  When  y  reaches  the  optimum,  y0,  which  satisfies  =  p[X(y0)Y(y0)], 
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the  matrix  Z  in  (2-7f)  will  become  infinity  since  the  matrix  I-)oY(y0)X(y0)  is  singular.  If  we 
try  to  apply  the  formulas  (2-7a)  -  (2-7h)  directly  to  construct  an  optimal  H“  controller,  a 
numerical  difficulty  will  arise  in  the  implementation  of  the  Ak  and  Ck  matrices.  We  will 
rearrange  these  formulas  such  that  an  optimal  H“  controller  can  be  constructed  without  any 
numerical  difficulty. 


The  dual  system  of  the  realization  in  (2-7a)  can  be  easily  rewritten  in  a  descriptor 
form.  The  state  equation  (generalized  state  equation)  of  the  descriptor  representation  can  be 
split  into  two  set  of  equations:  one  involves  first  derivative  of  some  state  variables  and  the 
other  is  just  an  algebraic  equation.  The  state  variables  which  have  no  derivative  in  the 
equations  can  be  eliminated  and  then  we  have  a  lower  order  state  space  representation  for 
the  dual  system.  The  dual  of  the  dual  system  is  identical  to  the  original  and  therefore  we 
have  an  optimal  H°°  controller  as  follows: 


where 


(2-8) 


\  =  [  v]adUi  *  V[ADU2  (V^ADU2)t  V^AdUj  ]  I'/  (2-9a) 

Bc  =  vX  -  v]adU2  (vX vj  vJbd  (2-9b) 

Cc  =  [  CDU!  -  CDU2  (vX U2)t  vXUi  1  Z]1  (2-9c) 

Dc  =  Dfc  -  CDU2  (vXU2)f  V2Bd  (2-9d) 

and 

Bd  =  -H2  +  (B2  +  H^)!^  (2- 10a) 

Cd  =  F2  -  Dk(C2  +  F12)  (2- 10b) 

Aj)  =  (B2  +  H12)  Cj)  +  (A  +  HC)  Ed  (2- 10c) 

Ed  =  I  -  Yo2  X(Yo)Y(Yo)  (2-10d) 


Zlt  Uj.Uj,  Vj,  and  V2  are  obtained  from  the  singular  value  decomposition  of  ED,  i.e„ 


U*1 


z,  0 
0  0 


(2-11) 


Similar  procedure  can  be  used  to  construct  optimal  H~  controllers  for  Case  (2)  and 
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die  formulas  in  Theorem  2. 1  can  be  used  to  construct  optimal  H~  controllers  for  Case  (1). 


2.3  Formulation  of  H“  Optimization  Problems 

Many  control  problems  can  be  formulated  as  the  standard  H“  optimization  problem. 
For  the  purpose  of  demonstration,  two  examples  are  given  in  the  following.  The  first  is  a 
mixed-sensitivity  optimization  problem;  the  second  is  a  disturbance  reduction  problem  with 
measurement  noise. 


A.  Mixed-Sensitivity  Optimization  Problem 


Consider  the  following  system: 

y(s)  =  P(s)u(s)  +  v(s)  (2- 12a) 

u(s)  =  K(s)y(s)  (2- 12b) 


where  v(s)  is  disturbance,  y(s)  is  output  and  K(s)  is  controller  to  be  designed.  It  is  well 
known  that  a  smaller  ll(I-PK)'1llos  means  a  better  disturbance  attenuation,  whereas  a  smaller 
IIPKCI-PK)'1  II  ^  implies  a  better  robust  stability.  Unfortunately,  the  H~  norms  of  (I-PK)1 
and  PK(I-PK)4  may  not  be  made  small  at  the  same  time.  If  we  make  one  of  them  smaller 
then  the  other  will  become  larger.  To  have  a  trade-off  between  these  two  quantities, 
Kwakemaak  [3]  formulated  the  mixed-sensitivity  problem  as  the  problem  of  finding  a 
controller  K(s)  which  stabilizes  the  closed-loop  system  and  minimizes  I  WML  where  d>  is 
given  by 


Wjd-PK)'1 

<D  * 

w2pk(i-pk)  1 


(2-13) 


Wj  and  W2  are  weighting  matrices  chosen  by  the  designer  according  to  the  concrete 

situation.  In  other  words,  they  depend  on  the  characters  of  the  disturbances  and  system 
uncertainties.  Usually,  the  disturbances  occur  most  likely  at  low  frequency,  therefore 
Wt(s)  is  chosen  to  be  a  low-pass  filter  to  emphasize  the  error  energy  at  low  frequency.  The 
plant  uncertainty  is  also  frequency-dependent;  the  higher  the  frequency  is,  the  larger  the 
uncertainties  become.  Hence,  W2(s)  is  usually  chosen  to  be  a  improper  transfer  function 
(but  W2P(s)  has  to  be  a  proper  transfer  function),  which  is  analytic  in  closed  right  half 
plane.  In  the  following,  we  assume  that  W^s)  is  strictly  proper,  W2(s)  is  a  polynomial 
such  thatW2P(s)  remains  proper  and  both  of  them  are  analytic  in  closed  right  half  plane. 
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The  problem  of  finding  a  K(s)  which  stabilizes  the  closed-loop  system  and 
minimizes  IHWI^  can  be  rearranged  into  the  standard  H~  optimization  problem.  Consider  the 

following  system: 


u  =  Ky  (2-14b) 

It  is  easy  to  show  that  the  matrix  <£>  defined  by  (2-13)  is  just  the  transfer  function  from  v  to 
[ZjT  z^]T  of  the  closed-loop  system  (2-14).  Comparing  (2-14a)  with  (2-1),  we  can  see 

that 


If  P,  W2P,  and  W1  have  state-space  realizations  as  follows 


Then  the  generalized  plant  G(s)  has  a  state  space  realization  as  shown  in  (2-4)  with 


c2  =  [CP  °] .  D21  =  I ,  D22  =  Dp  (2-17) 

Note  that  because  W2  is  a  polynomial,  the  A-matrix  of  W2P  is  same  as  that  of  P. 


B.  Disturbance  Reduction  Problem 
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Fig.2.1  A  disturbance  attenuation  problem 


Consider  the  feedback  system  shown  in  Fig.2.1.  P(s)  is  a  given  plant,  W^s), 
i=l,2,3,4  are  weighting  matrices,  and  K(s)  is  the  controller  to  be  designed.  The 
disturbance  and  noise  are  the  outputs  of  W3  and  W4  driven  by  d  and  n  respectively,  zj 
the  weighted  error  response;  z2  is  the  weighted  control  input  Let  zT  =  [ZjT  z*]1,  vT  =  [g 

nT]T  and  assume  that  v  is  unknown  but  with  its  energy  bounded  by  unity.  The  objective  is 
to  find  a  controller  K(s)  which  stabilizes  the  closed-loop  system  and  minimizes  the  worst 
l(zll2,  i.e.,  minimizes  the  H“  norm  of  T^,  the  closed-loop  transfer  function  from  v  to  z.  Tw 
is  given  by 


W1Pa-KP)'1W3  W1PKa-PK)‘1W4 
W2KPa-KP)1W3  W2K(I-PK)'1W4 


Note  that  W1PK(I-PK)  lW4  and  W2KP(I-KP)1W3)  are  the  output  and  input 
complementary  sensitivity  functions.  Their  H*°  norms  indicate  the  stability  robustness  of 
the  closed-loop  system  for  the  multiplicative  plant  uncertainty  introduced  at  the  output  and 
input  respectively.  W2K(I-PK)1W4  is  the  control  complementary  sensitivity  function 
whose  H~  norm  indicates  the  stability  robustness  of  the  closed-loop  system  for  additive 
plant  uncertainty.  Hence,  reducing  IIT^IL  will  also  improve  the  robust  stability  of  the 
closed-loop  system. 

It  is  easy  to  verify  that  the  generalized  plant  of  the  system  can  be  expressed  as: 


V 

'WjPWj 

0 

W,p" 

‘d" 

*2 

= 

0 

0 

W2 

n 

_y . 

[PW3 

W4 

p  . 

_  u  _ 

That  is. 


(2-19) 
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'11 


WjPWj  0 

0  OJ 


G12  = 


WjP 

w„ 


G21=[PW3  w4] 


G22  =  P 


If  P,  Wif  i=l,2,3,4  have  state-space  realizations  as  follows: 

i=  1,2,3, 4 


P  = 


bp1 

.  cp 

D-J 

W  = 
1 


■^wi 

Bwi 

C* 

Dwi_ 

Then  the  generalized  plant  G(s)  has  a  state  space  realization  as  shown  in  (2-4)  with 


bp  ' 

BwlDP 


["  Ap  0  0  BpC.3  0-1 

|BpDw3  0  1 

r 

A  = 

BwlGP  Awl  ®  BwlDPGw3  0 

Bi- 

B.,DrDw3  0 

b2  = 

o 

o 

Jf 

o 

o 

0  o 

o 

o 

© 

£ 

© 

0 

<? 

© 

o 

o 

© 

_ 1 

♦ 

- 1 

1 

CQ 

o 

_ 1 

9  W 

cl= 


DwiCp  cwi  0  DwlDpCw3  0 


0  0  c 


w2 


0 


0 


Dll= 


^wl^P^w3  B 

0  OJ 


D12  = 


Bw2 

0 

0 

DwlD 

L  DW2 


C  =rcn  0  0  D  C  .  C  J  D,  =rD„D, 

2  L  P  P  w3  w4j*  21  (_  P  w3 


Dw4], 


D22  ~  DP  • 


(2-20) 


(2-21) 


p 

(2-22) 


Above  {A3,CJJ}  is  the  state-space  representation  for  the  generalized  plant  G(s). 


CHAPTER  3 


STRUCTURE  OF  THE  y-DOMAIN  FOR  THE 
DGKF  H**  RICCATI  SOLUTIONS 


3.1  Introduction 

The  major  computational  burden  involved  in  the  two  y-dependent  H*°  Riccati 
equation  approach  is  the  computation  of  the  optimal  H~  norm.  To  develop  efficient 
algorithm  for  computing  the  optimal  H**  norm,  it  is  necessary  to  investigate  how  die  Riccati 
solutions  vary  as  a  function  of  the  paramter  y. 

There  have  been  several  attempts  to  address  these  issues.  Pandey  et  al.'s  hybrid 
gradient-bisection  method  [13]  and  Chang  et  al.'s  double  secant  and  bisection  method  [14] 
were  proposed  for  the  computation  of  the  optimal  H*°  norm.  In  these  two  papers,  a 
conjecture  was  made  that  the  spectral  radius  of  the  product  of  the  two  Riccati  solutions  is  a 
convex  function  of  y2.  Base  on  this  conjecture,  improved  search  schemes  were  developed. 
However,  since  there  was  no  proof  for  this  conjecture,  bisection  was  used  in  these  two 
algorithms  as  a  supplement  to  guarantee  convergence.  In  [IS],  Scherer  considered  the 
inverse  of  the  H“  Riccati  solutions,  defined  a  new  independent  variable  p.  =  y‘  ,  and 
showed  that  these  inverses  are  concave  functions  of  p  in  matrix  sense.  Based  on  this  fact,  a 
Newton-like  algorithm  was  proposed  to  compute  the  optimal  H*°  norm  when  an  appropriate 
starting  point  is  available.  The  convexity  of  the  spectral  radius  of  the  product  of  the  two 
Riccati  solutions  with  respect  to  y  was  first  shown  by  Li  and  Chang  [16].  With  this 
convexity  property,  the  optimal  H~  norm  can  be  easily  computed  by  a  quadratically 
convergent  algorithm  if  an  appropriate  starting  point  is  available.  However,  it  is  not  trivial 
to  obtain  such  a  starting  point  that  can  guarantee  the  convergence.  Furthermore,  no  general 
and  rigorous  analysis  is  available  for  the  y-dependent  H~  Riccati  equations. 

In  this  chapter,  we  investigate  the  structure  of  the  y-domain  where  the  FT  Riccati 
solutions  exist  and/or  positive  semidefinite,  and  reveal  some  useful  properties  of  the  H~ 
Riccati  solutions  on  the  y-domain,  such  as  continuity,  monotonicity,  and  convexity  with 
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respect  to  y.  These  results  provide  a  better  understanding  of  the  state  space  approach  to  H~ 
optimization  problems.  Beside  their  theoretical  value,  these  properties  can  be  employed  to 
develop  efficient  algorithms  for  partitioning  the  Y-domain,  finding  the  smallest  y  such  that 
the  two  Riccati  solutions  are  positive  semidefinite,  and  computing  the  optimal  H~  norm. 
The  chapter  is  organized  as  follows.  Section  3.2  lists  the  notations  and  reviews 
preliminaries.  In  Section  3.3,  we  investigate  the  structure  of  Y-domain  of  the  IT*  Riccati 
solutions.  The  continuity,  monotonicity,  and  convexity  properties  of  the  Riccati  solutions 
are  discussed  in  Section  3.4. 


3.2  Notations  and  Preliminaries 


E 

R+ 

inf(S) 
sup  (5 ) 
£mxn 
E^s) 
I 
In 

At  or  A' 
A* 

G-(s) 

A£0 

A>0 

MA) 

Xd(A) 

^ma(A) 

^min(A) 

P(A) 


The  set  of  real  numbers. 

The  set  of  positive  real  numbers. 

The  infimum  of  the  set  S. 

The  supremum  of  the  set  S. 

The  set  of  mxn  matrices  whose  entries  are  complex  numbers. 
The  set  of  mxn  transfer  matrices. 

Identity  matrix, 
n-dimensional  identity  matrix. 

Transpose  of  a  matrix  A. 

Conjugate  transpose  of  a  matrix  A. 

GT(-s). 

A  is  a  positive  semidefinite  matrix. 

A  is  a  positive  definite  matrix. 

An  arbitrary  eigenvalue  of  A. 

The  ith  eigenvalue  of  A. 

The  maximal  eigenvalue  of  A. 

The  minimal  eigenvalue  of  A. 

The  spectral  radius  of  A. 
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CW  or  a 

Mi) 

My) 

Mi) 

A1  (7) 


The  maximal  singular  value  of  A. 

The  matrix  A  is  considered  as  a  function  of  ye  R+. 


The  first  derivative  of  A(y)  with  respect  to  y,  Le., 


Hie  second  derivative  of  A(y)  with  respect  to  y,  i.e.. 


dA(y) 

dy 

d2A(y) 


dy2 


The  first  derivative  of  A'^y),  i.e.,  —(A  J(y)). 


d2 


A1  (Y)  The  second  derivative  of  A'^y)  with  respect  to  y,  i.e.,  ^(A  (y)). 


Throughout  die  chapter,  the  notation  P(s)  = 


or  P(s)  =  {A,  B,  C,  D}  is 


used  to  represent  a  state  space  realization  of  a  system  whose  transfer  function  is  P(s)  =  C(sl 
-  A)_1B  +  D.  The  HT  norm  of  P(s)  is  defined  as 


II  Pll  :=  sup  o  [P(jo)>] 
os 

Many  control  problems  can  be  formulated  as  the  following  standard  H*° 
optimization  problem.  In  the  standard  H“  optimization  problem  formulation,  the  system 
representation  is  rearranged  as  follows. 


'z(s)' 

'Vs) 

G12(S)1 

w(s)' 

_y(s) 

Vs) 

/-"■v 

& 

u(s) 

:=  G(s) 


w(s) 

u(s) 


(3-D 


where  z,  y,  w,  and  u  are  the  controlled  output,  the  measured  output,  the  exogenous  input, 
and  the  control  input  respectively.  The  controlled  output  vector  z  usually  includes  the  error 
signal  and  a  weighted  control  input  The  exogenous  input  w  contains  the  disturbances,  the 
noises,  and  the  commands.  The  measured  output  vector  y  consists  of  all  the  signals  which 
can  be  measured  and  available  for  feedback.  Through  the  control  input  u,  the  behavior  of 
the  system  can  be  modified.  The  vector  y  will  be  used  as  the  input  to  a  controller  K(s)  and 
the  output  of  K(s)  will  be  connected  to  the  control  input  u.  That  is, 

u(s)  =  K(s)  y(s)  (3-2) 

The  standard  H~  optimization  problem  is  the  problem  of  finding  a  proper  controller  K(s) 
such  that  the  closed-loop  system  is  internally  stable  and  II  Tw !!_  is  minimized  where 
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(3-3) 


Tw(s)  =  G n(s)  +  G12(s)  K(s)  [  I  -  G^s)  K(s)  ]  4  G21(s) 


That  is,  Tw(s)  is  the  transfer  function  of  the  closed-loop  system  from  w  to  z. 
Let  a  realization  of  die  generalized  plant  G(s)  be 


G(s)  = 


On(S)  GU(S)_ 

.<*»«  <y». 


with  the  following  assumptions: 

(i)  Both  G12(s)  and  G21(s)  have  no  transmission  zeros  on  the  jeo  axis. 

(ii)  (A,  B2)  is  stabilizable  and  (Q,  A)  is  detectable. 

(HODylCi  D12]  =[0  I] 


(3-4) 


(3-5a) 

(3-5b) 

(3-5c) 

(3-5d) 


The  two  Riccati  equations  involved  in  die  state  space  approach  are: 

A’X«+X^A  +  XM(Y'2BI^'-B2B2)X«+qCi  =  0  (3-6a) 

and 

A  Y„  + Y„A'  + Y„(  i1  qCi  -  Cf2  )  Y.  +  BiBj  =  0.  (3-6b) 

The  following  theorem  by  Doyle  et  al.  characterizes  all  suboptimal  stabilizing 
controllers  such  that  II  Tw  ll_  <  y. 


Theorem  3.1:  [7]  There  exists  a  stabilizing  controller  such  that  IITjJL  <  Y  if  and  only  if 
the  following  three  conditions  hold: 

0)  There  exists  a  positive  semidefinite  stabilizing  solution  X_(y)  to  (3-6a). 

(ii)  There  exists  a  positive  semidefinite  stabilizing  solution  Y„(y)  to  (3-6b). 

(iii)  p[X«(Y)Y.(Y)]<V2. 

When  these  conditions  hold,  one  such  controller  is 


where 
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At  =  A  +  y  2B1BjX_  +  B2F  +  ZLC, 

F  =  -%X_,  L>-Y_q,  E  =  ( I  -  y2Y«X.  ),  Z  =  ( I  -  f2Y_X_ 


The  major  computational  task  of  H~  design  is  to  solve  the  two  Riccati  equations, 
eq.(3-6a&b).  From  [7],  the  stabilizing  solution  X„  and  Y„  can  be  solved  via  their 
corresponding  Hamiltonian  matrices 


H Jy)  = 


'^BlB'l'  ®2B2 


(3-7a) 


Ui)  = 


A' 


-B  B 
11 


Y  C  C  -  C  C 
1  11  2  2 


(3-7b) 


Namely,  if  we  let 


(3-8) 


be  a  modal  matrix  for  HL,  corresponding  to  all  the  eigenvalues  on  the  open  left  half  of  s- 
plane,  then 

X_  =  T2T'11,  (3-9) 

is  the  unique  stabilizing  solution  to  eq.(3-6a)  such  that  (y^BjB^  -  B2B2  )X„  is  stable  [1,7]- 

It  is  well-known  [7]  that  the  stabilizing  X«,  exists  if  and  only  if  (1)  H„  has  no  j co-axis 
eigenvalues,  and  (2)  T}  is  invertible.  These  two  conditions  are  called  stability  condition  and 
complementary  condition  respectively  in  [7].  We  have  the  same  arguments  for  Y„,  the 
stabilizing  solution  to  eq.(3-6b). 


3 3  Structure  of  the  y-Domain  of  H°°  Riccati  Solutions 

It  is  well-known  that  for  the  standard  LQG  Riccati  equations  [17],  positive 
semidefinite  stabilizing  solutions  exist  if  the  system  is  stabilizable  and  detectable. 
Unfortunately,  the  H~  Riccati  equations  are  more  complicated  because  they  are  parameter 
dependent  and  their  quadratic-term  coefficient  matrices  may  not  be  sign  definite. 

For  a  sufficient  large  y,  a  suboptimal  H~  controller  always  exists.  This  is  due  to  die 
fact  that  as  y— *»,  the  H~  Riccati  equations  eq.(3-6a&b)  become  the  standard  LQG  Riccati 
equations  (the  quadratic  terms  become  sign-definite),  therefore  both  X.  and  Y„  exist  and 
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are  positive  semidefinite.  Meanwhile,  since  y  is  sufficiently  large,  the  third  condition  of 
Theorem  3.1  is  also  satisfied.  However,  as  y  decreases,  the  three  conditions  of  Theorem 
2.1  may  not  hold  any  more.  For  instance,  X„  may  not  exist  or  fail  to  be  positive 
semidefinite  for  some  small  y.  Thus  the  following  problems  arise: 

1.  How  small  the  y  could  be  such  that  the  Riccati  equations  have  stabilizing  solutions? 

2.  What  is  the  structure  of  y-domain  where  the  solutions  exist? 

3.  How  to  compute  the  optimal  y,  which  is  defined  as 

Yop,  =  inf  (y:  all  the  three  conditions  of  Theorem  3.1  hold}. 

We  will  answer  these  questions  by  starting  the  investigation  of  the  structure  of  y- 
domain  for  XJi). 

Theorem  3.2:  If  we  define 

ctx  :=  inf  {y :  ye  R+  and  XJy)  exists) 

Px  :=  inf  [y :  ye  R+  and  XJy)  is  positive  semidefinite}, 

then  we  have  the  following  results: 

a)  On  (GLx,  +*»),  XJy)  is  well  defined  almost  everywhere; 

b)  On  (Px,  -h»),  X„(y)  is  well  defined  and  positive  semidefinite  everywhere; 

c)  ax  £  Px* 

Fig. 3.1  gives  an  illustration  for  the  theorem.  One  can  see  that  the  y-domain  for 
X^(y)  consists  of  three  parts  which  are  bounded  by  a,  and  px.  For  any  y  in  Region  c, 
X„(y)  is  positive  semidefinite;  whereas  in  region  a,  XJy)  does  not  exist  at  all.  In  Region  b, 
although  XJy)  is  not  positive  semidefinite,  it  exists  everywhere  except  for  some  possible 
isolated  points. 


Fig.  3.1  y-domain  of  X^Cy). 


In  order  to  simplify  the  proof  of  Theorem  3.2,  we  will,  without  loss  of  generality, 
strengthen  die  assumption  on  (Cj,  A)  from  detectable  to  observable.  To  show  that  this  can 
be  done,  we  assume  that  (Cj,  A)  is  not  observable.  Then  one  can  always  restrict  the 
problem  into  the  observable  subspace  by  finding  an  orthogonal  matrix  [18] 

U=[Uj  U2]  (3-10) 
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such  that 


and 


U'AU  = 


U’jAUj  0 

u’2au1  u'2au2 


U’B  = 


u;b1  u;b2 
UiBl  U2B2 


CjU  =  [CjUi  0] 


(3-11) 

(3-12) 

(3-13) 


with  (CiUi,  UjAUi)  observable.  It  is  easy  to  verify  that  the  solution  to  eq.(3-6a)  can  be 
expressed  as 


X-(y)  =  U 


X(Y)  0 
0  0 


U\ 


(3-14) 


where  X(y)  is  the  stabilizing  solution  to  the  eq.(3-6a)  with  (A,  Bi,  B2,  CL)  replaced  by 
(UjAUi,  UjBi,  U'^,  C1U1).  Note  that,  since  the  matrix  U  in  eq.(3-14)  is  independent 
of  y,  Xm(y)  exists  if  and  only  if  X(y)  does,  and  X_(y)  is  positive  semidcfmitc  if  and  only  if 
X(y)  is.  Hence,  in  the  following,  we  will  assume  that  (C]t  A)  is  observable  and  use  X_(y) 
and  X(y)  interchangeably  without  loss  of  generality. 

To  prove  Theorem  3.2,  we  need  die  following  lemmas. 

Lemma  3.1  [7]:  The  stabilizing  solution  to  the  H~  Riccati  equation  eq.(3-6a)  can  be 
written  in  terms  of  two  LQG  Riccati  solutions: 

XJY)  =  X^  X2  -  W(y)  ]%,  (3-15) 

where  X2  is  the  stabilizing  solution  to: 

A'X2  +  X2A  -  X2B2B2  X2  +  Ci  Ci  =  0  (3-16) 

and  W(y)  is  the  stabilizing  solution  to: 

(A-B2B'2X2)^  +  W(A-B2B’2X2)  +  WX2  C’CjX'1 w  +  y'XBjBjX^  0  (3-17) 


whose  corresponding  Hamiltonian  matrix  is 


(3-18) 


(a-b2b2x2)  X^C’CjX^ 

Hw(Y)  := 

-yXb.b.x,  -(A-BBXV 
'  2  1  1  2  2  2  2  ] 

Since  X„(y)  can  be  expressed  in  terms  of  W(y)  and  a  constant  matrix  X2,  the 
properties  of  X„(y)  can  be  derived  from  those  of  W(y).  Lemma  3.2  describes  the  y-domain 
of  W(y),  and  Lemma  3.3  shows  some  useful  properties  of  W(y). 

Lemma  3.2:  There  exists  a  positive  number  aw  such  that  W(y)  is  well  defined 
everywhere  on  (a*,  +<»),  and  W(y)  does  not  exist  on  (0,  a*). 


Proof:  First  of  all,  it  is  claimed  that  all  the  jo>-axis  eigenvalues  of  H^y)  are  identical  to  the 
j to- axis  transmission  zeros  of  die  transfer  function 

r(s,y):=Iy2-G;(s)Gw(s)  (3-19) 

where 

Gw(s)  =  B\x2[  Is  -  (A-B2B2X2)  t[ls~X  ft.  (3-20) 

This  can  be  shown  by  the  following  derivations: 
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By  the  similarity  transformation 


iy1 

T  = 

we  have 


A  66' 

y1! 

r'(s,  D  - 

-Y =tt  -A' 

0 

_  o  Y‘fi' 

iy2 

(3-22) 


The  last  equation  shows  that  {  H^y),  [y^CjXj1  0]\  [0  yk^X’1],  Iy2  }  is  a  realization  of 
r\s,  y).  Since  A  =  (A-B  B'XJ  is  stable,  by  applying  PBH  test  [19],  it  is  clear  that  on  the 

*  *  Z 

j  co-axis  we  have 


A-Is 

A  A 

BB' 

y^B 

A -Is 

0  y^B 

rank 

-y2C’C 

A 

-A’-Is 

0 

=  rank 

-y2C'C 

-A’-Is  0 

A -Is 

A  A 

BB* 

A 

A-Is 

0 

rank 

-y'2C’C 

-A’ -Is 

=  rank 

9  A  A 

-y2C  C 

-A’  -Is 

0 

y'16' 

0 

y  1B' 

This  implies  that  the  above  realization  is  both  controllable  and  observable  on  the  jeo-axis. 
That  is,  the  realization  does  not  have  pole-zero  cancellation  on  die  joo-axis.  Hence,  the  set 
of  the  eigenvalues  of  H^y)  on  the  jeo-axis  is  identical  to  the  set  of  the  jo-axis  poles  of 

r  \s,  y),  which  in  turn  is  identical  to  the  set  of  die  jeo-axis  transmission  zeroes  of  T(s,  y). 

Next,  we  claim  that  if  there  exists  a  y  such  that  T(s,  Yj)  has  transmission  zeros  on 
die  jeo-axis,  then  r(s,  y)  will  have  transmission  zeros  on  the  jeo-axis  for  any  ye  (0,  Yj].  To 
prove  this,  one  just  needs  to  note  that  Gw(s)  is  a  strictly  proper  transfer  function  and 
therefore  it  goes  to  zero,  as  co  -»  -h».  Hence,  for  any  y  <  ylf  there  always  exists  an  co  such 
that  r(jeo,  y)  is  singular.  For  instance,  in  Fig.3.2,  we  can  see  that  if  r(s,Yj)  has  a 
transmission  zero  at  joy,  then  for  any  y2  <  yP  there  exists  an  co2  such  that  jco2  is  a 
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transmission  zero  of  r(s,  y). 


0  ©i  ® 

Fig.  3.2  The  j co-axis  transmission  zeros  of  r(s,  y)  occurs  at  the  frequencies 
such  that  y =  o(Gw(jco)). 

These  two  claims  imply  the  fact  that  if  there  exists  a  yx  such  that  H^Xy)  has  j  co-axis 
eigenvalues,  then  H^y)  has  jco- axis  eigenvalues  for  any  yless  than  yr  Therefore,  if  we 
define 

a*  =  inf  {y.  H^(y)  has  no  j  co-axis  eigenvalues},  (3-23) 

then  it  is  easy  to  see  that  H^y)  has  jco-axis  eigenvalues  for  any  ye  (0,  c^)  and  has  no  jco- 
axis  eigenvalues  far  any  ye  (aw,  +*»). 

Finally,  since  A-B2B^X2  is  stable  and  H^Cy)  has  no  jco-axis  eigenvalues  for  any 
Ye  (a*,  +°°)» then  W(y),  the  stabilizing  solution  of  H^y),  exists  for  any  ye  (txw,  +*»)  [20]. 
Thus,  we  complete  die  proof. 

From  the  above  proof,  it  is  easy  to  see  that  the  aw  defined  in  (3-23)  can  be 
expressed  by  die  H*°  norm  of  a  known  transfer  function: 

Ow  =  II  Gw(s)  ll„,  (3-24) 

where  Gw(s)  is  defined  in  (3-20). 

Lemma  3.3:  On  (aw,  -h»), 

a)  WY  :=  dW/dy  ^0,Ww  :=  d2W/dy2  SO; 

b)  all  eigenvalues  of  W(y)  are  analytic,  nonincreasing  functions  of  y; 

c)  W(y)  S  0. 


Proof:  Since  W(y)  is  an  analytic  function  of  yon  (oc*,  -h»),  it  is  well-known  [21]  that,  by 
appropriate  ordering,  the  eigenvalues  {X^}  and  eigenvectors  {vi}  of  W(y)  are  analytic 
functions  of  y  on  (aw,  -h»)  such  that  (IXi-W(y))vi  =  0,  i  =  1, 2,—,  n.  Rewrite  Riccati 
equation  (3-17)  as: 

A*W  +  WA  +  wfcfcw  +  p£'6  =  0  (3-25) 

where  p  =  y*2,  A  =  A-B  B’X  ,  §  =  X^C'  and  6=  B*  X2.  Differentiating  the  above 
equation  with  respect  to  p,  we  have  die  following  Lyapunov  equation: 

(A  +  +  W^A  +  fc&’W)  + 1 1  =  0  (3-26) 


Q  W  A  A  A 

where  Because  A  +  BB'W  is  stable,  it  can  be  inferred  that  W^>  0  and  is 

the  unique  solution.  Continuously  differentiating  eq.(3-26)  with  respect  to  p,  we  have 


d*W 

0,  where  =  — 5".  By  using  the  chain  rule  of  differentiation,  it  yields 


m 


dp 


w„  -  ■ +  4y\/mi0. 

Now  let  us  consider  equation 
(IX  -  W)v  =  0, 


and 


(3-27) 

(3-28) 


(3-29) 


where  X  is  any  eigenvalue  of  W  and  v  a  corresponding  eigenvector.  Differentiating  eq.(3- 
29),  we  have 

vW  v 

X(Y)=-^0,  (3-30) 

since  ^  0.  Thus  we  proved  part  b). 

Part  c)  is  a  direct  consequence  of  part  b).  Since 


Hw(oo)  = 


(a-bb'xj  x^c'cx;1 
v  2  2  2  2  1  1  2 


0 


the  stabilizing  solution  W(<»)  =  0.  Hence,  it  can  be  concluded  that  W(y)  £  0  on  (c^,  +°°) 
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from  the  fact  that  all  the  eigenvalues  of  W(y)  are  nonincrcasing. 


In  the  proof  of  Theorem  3.1  we  need  to  show  that  there  exist  at  most  n  isolated 
values  of  y  such  that  [X2  -  W(y)]  is  singular,  let  Xj(y)  be  the  ith  eigenvalue  of  [X2  -  W(y)],  i 
=l,2,...n,  where  n  is  the  dimension  of  W(y).  Then  we  have  the  following  lemma. 

Lemma  3.4:  If  there  is  a  Yie  (a*,  +®°)  such  that  XjCyj)  =4),  then  Xi(y1)>0. 

Proof:  Define  a  function  on  (a*,  +«)  as 

f(y):=  x’[X2-W(y)]x,  (3-31) 

where  x  =  v^)  is  an  eigenvector  corresponding  to  X^).  By  Lemma  3.3(a),  it  is  easy  to 
see  that  f(y)  is  a  nondecreasing  and  concave  function  of  y  on  (ax,  -k»).  Furthermore,  we 
have  f(+~)>0,  since  [X2  -  W(+«>)]  =  X2  >0  under  the  assumption  of  (Ci,  A)  being 
observable  [17].  Hence,  if  there  is  a  Y]€  (a,,  -h»)  such  that  f(yj)  =  0  then  f(Y2)>0.  By  this 

property,  we  can  see  that  XjfyjH)  implies  f(Yj)>0,  which  in  turn  implies  A..(y1)>0- 

Together  with  the  monotonicity  of  Xj(y),  it  comes  the  conclusion  that  [X2  -  W(y)] 
has  at  most  n  singularity  points  on  (a*,  +°°). 

With  these  lemmas,  now  we  are  ready  to  prove  Theorem  3.2. 

The  Proof  of  Theorem  3.2: 

a)  Since  XJi)  ^XJ  X^-  W(y)  ]_1X2  (from  Lemma  3.1),  it  is  obvious  that: 

OxsOw  (3-32) 

where  a*  was  defined  in  Theorem  3.2  and  a*  was  defined  in  Lemma  3.2.  It  is  possible  for 
[Xj  -  W(y)]  to  lose  rank  at  some  points  on  (a*,  -k»),  therefore  XJy)  is  not  well  defined  at 

these  points.  In  Lemma  3.4,  we  show  that  there  are  at  most  n  such  isolated  points  where  n 
is  the  dimension  of  XM. 

b)  Suppose  that  there  exists  a  Yj  such  that  X<o(y1)  >  0,  and  hence  [Xj-Wfy)]  >  0. 
Because  all  the  eigenvalues  of  W(y)  are  nonincreasing  (from  Lemma  3.3),  [  X2  -  W(y)  ] 
will  not  become  singular  for  ally  ^y^  which  implies  that  XJy)  keeps  positive  definite  for 
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all  y  €  fy,  +°o).  Therefore,  if  we  define 

Px  =  inf  {y.  XJj)  is  positive  definite},  (3-33) 

then  die  proof  is  completed. 

The  following  theorem  characterizes  P*  and  shows  how  to  compute  it. 

Theorem  33:  Define  f(y)  :=  Xnjax[W(7)X'21]  on  (o^,  -h»),  dien 

a)  either  f(Px)  =  1  or  Px  =  o^; 

b)  f(y)  is  a  convex  decreasing  function  of  y; 
where  denotes  the  maximal  eigenvalue. 

Proof :  The  proof  of  part  a)  is  trivial,  which  is  simply  due  to  the  fact  that  (Xj  -  W(y))  >  0 
if  and  only  if  Xn^fWO/JX.1]  <  1.  To  prove  part  b),  we  consider  the  following  equation 

[XX2-W(y)]u  =  0.  (3-34) 

Note  that  X  in  the  above  equation  is  an  eigenvalue  of  W(y)X2\  Taking  derivative  with 
respect  to  yon  cq.(3-34),  we  have 

(XX2-WT)u  +  lXX2-W(|)]u  =0  (3-35) 

It  is  easy  to  show  that 

.  u’W^i 

(3‘36) 

on  (ax,  +oo),  since  W^,<  0  from  Lemma  3.3.  If  we  replace  X  in  eq.(3-35)  by  X,,^,  the 
maximal  eigenvalue  of  W(y)X2\  then  we  have 

u  X  imnX2  -  Wy )  u  =  -  u  '[Xm„X2  -  W(y) )  u  (3-37) 

where  u  is  the  eigenvalue  corresponding  to  the  Xn,„.  Note  that  eq.(3-35)  implies 

u'(  XmaxX2-W7)u  =  0.  (3-38) 

Taking  derivative  on  eq.(3-38)  yields 

u  X  XmixX2  -  WT  )  u  +  uX  imMCX2  -  Wy  )u  +  uxU  X2-Wyr)u  =  0.  (3-39) 

Combining  eq.(3-37)  and  eq.(3-39),  we  have 


25 


(3-40) 


uW  u  +  2u'(Xm  X  -  W(Y)]u 
Wy)=  ^  u.x-u-2- - 

Since  (XmaxX2  -  W(y))  ^  0  and  0  from  Lemma  3.3,  we  have  X^Cy)  £  0,  which, 
together  with  eq.(3-26),  implies  that  f (y)  is  a  convex  decreasing  function  of  y  on  (a,,  -h»). 

Corollary  3.1:  The  derivative  of  f(y)  can  be  computed  by 
dy  w  v 

where  w  and  v  are  right  and  left  eigenvectors  of  W(y)Xj1  corresponding  to  its  maximal 
eigenvalue  respectively,  and  can  be  computed  from  eq.(3-26)  and  eq.(3-27). 

Fig.3.3  illustrates  two  possibilities  mentioned  in  Theorem  3.3.  As  shown  in  the 
figure,  jix  can  be  obtained  either  by  finding  the  value  of  y  such  that  f(y)=l  or  simply  by 
computing  a*. 


HD 


Fig.  3.3  Two  possible  locations  of  Px. 


This  Corollary  suggests  that  one  can  develop  an  algorithm  to  compute  (i*  by  making  use  of 
these  properties. 

3.4  Properties  of  the  H~  Riccati  Solutions 

In  this  section,  some  important  properties  of  X^Cy)  on  the  interval  (jJ*,  -k»)  will  be 

given.  We  are  only  interested  in  this  interval,  since  the  H~  problem  has  no  solution  for  any 
y  less  than  J3*  ( see  Theorem  3.1 ).  The  properties  are  summarized  in  the  following. 

2 

Theorem  3.4:  On  (px,  +°°),  Xjy)  :=  ■JjCX-CY))  £  0,  XJj)  :=  >  0. 
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Proof:  LetS(Y):=[X2-W(Y)],  S:=^{S(tf}  and S^tSCy)1},  we have 

s  *-wr  S'^-S^SS^-S-'WjS4. 

From  Lemma  3.1, 

x_(Y)  -  ^  (Xj[  Xj  -  W«  ]‘x2)  =  Xj  S1  Xj 

«  X.s'w^s'1  X,  (3-42) 

This  implies  that  Xj(y)  £  0,  since  WY  £  0  (see  eq.(3-27)).  Continuously  taking  the 
derivative  on  eq.(3-42),  we  have 

xjy) = x2s'1wyS‘1x2+x2s-1wtyS‘1x2+x2s-1wys'1x2 

=  x2s-1wYs*1wYs‘1x2+x2s-1wYys-1x2+x2s-1wYs-1w1s*1x2 

=  X2S-1[2WTS'1WY+WYylS-1X2  (3-43) 

This  equation,  together  with  eq.(3-28),  indicates  that  Xjy)  -  0.  Thus  the  proof  is 
completed. 

From  these  inequalities,  we  have  the  following  properties  for  Xjy). 

Theorem  3.5:  On  (Px,  +«»), 

a)  all  eigenvalues  of  XJy)  are  analytic,  nonincreasing  functions  of  y  ; 

b)  the  maximal  eigenvalue  of  X^Cy)  is  a  nonincreasing,  convex  function  of  y. 

Proof:  a)  Since  XJy)  is  an  analytic  function  of  y  on  (P*,  +«>),  it  is  well-known  [21]  that, 
by  appropriate  ordering,  the  eigenvalues  {A*}  and  eigenvectors  {vj}  of  X„( y)  are  analytic 
functions  of  y  on  (pj,  +*»)  such  that 

(iAi  -  X„(Y»Vi  =  0  i  =  1, 2,—,  n.  (3-44) 


To  prove  the  monotonicity,  we  take  derivative  on  eq.(3-44): 
(IXj  -  X«)vi  +  (IXj  -  Xco(Y))vi  =  0, 


(3-45) 


.  v!  X^Vj  • 

which  implies  X.(Y)  =  — L  <,  0,  since  X^<,0.  Thus  we  proved  part  a). 

1  v.  v. 

!  I 

b)  From  eq.(3-45)  with  Xi  =  X^,*,  the  maximal  eigenvalue  of  X„(y),  and  Vj  *  v, 
the  eigenvector  corresponding  to  we  have  the  following  equation: 


vXjr  +  ti'VX^-XJv 

v’v 


(3-46) 


which  implies  K^fy)  £  0,  since  3c»  ^  0  and  (IX^-X^)  £  0.  Because  X^^Cy)  £  0  and  all 
eigenvalues  of  X„(y)  under  appropriate  ordering  are  analytic  and  nonincreasing  on  (Px, 
-h»),  it  is  easy  to  see  that  X„,u(y)  is  convex  and  nonincreasing  on  ({lx,  +»),  even  though 
Xn,n(Y)  may  not  be  smooth  at  some  isolated  points.  Fig.  3.4  gives  an  example  showing  the 
possibility. 


Fig.  3.4  Smoothness  and  convexity  of  the  largest  eigenvalue. 

In  the  figure,  both  Xi  and  X2  are  smooth  functions  of  y,  but  Xmtx  is  not  smooth  at  Yi- 
However,  it  is  still  a  convex  function. 


Since  Riccati  equation  eq.(3-6b)  is  dual  to  eq.(3-6a),  we  have  similar  results  for 
Y„(y).  with  the  following  corresponding  notations: 


and 


oty  :=  inf  {y:y<=  R+  and  YJy)  exists} 

Py  :=  inf  { Y :  Y€  R+  and  YJy)  is  positive  semiderinite}, 


yjy)=v 


Y(y)  O' 

0  oJv’ 
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(3-47) 


Y(y)  and  V  are  defined  in  a  similar  way  as  those  for  Xjtf)  in  eq.(3-14). 

Recall  that  the  condition,  p[X«.(y)Y..(y)]  <  y2,  is  required  in  Theorem  3.1.  To 
investigate  die  properties  of  X»(y)Y„(y),  we  define 

a  :=  maxta*,  <Xy} 

P  :=  max{(3x,  Py}. 

Then  it  is  easy  to  infer  that  a)  X^(y)Y^(y)  exists  on  (a,  -h»)  almost  everywhere;  b) 
X-(y)Y„(y)  has  no  negative  eigenvalues  on  0,  +»),  since  both  X^(y)  and  Y«,(y)  are 
positive  semidefinite  on  0,  +«•).  Moreover,  the  eigenvalues  of  X-(y)Y^(y)  have  the 
following  properties. 

Theorem  3.6  [16]  :  On  0,  +®°), 

a)  all  eigenvalues  of  X.Y.  are  analytic,  nonincreasing  functions  of  y ; 

b)  p(y):=p[X«,(y)Y„(y)]  is  a  nonincreasing,  convex  function  of  y. 
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CHAPTER  4 


MONOTONICITY  AND  CONVEXITY  OF  THE  GD 
PARAMETER  DEPENDENT  H~  RICCATI  SOLUTIONS 

4.1  Significance  of  the  monotonicity  and  convexity  properties 

In  this  chapter,  we  will  investigate  some  important  properties  of  the  H**  Riccati 
solutions  which  play  a  key  role  in  the  solution  of  the  H“  optimization  problem.  It  is  well 
known  that  Doyle,  Glover,  Khargonekar,  and  Francis  (abbr.:  DGKF)  [7]  presented  a 
celebrated  two-Riccati-equation  type  solution  to  a  standard  H~  control  problem.  Through 
solving  two  H**  Riccati  equations,  an  optimal  (or  suboptimal)  stabilizing  H*°  controller  can 
be  easily  constructed.  Unlike  the  constant  Riccati  solutions  in  Linear  Quadratic  Gaussian 
(LQG)  or  H2  problem  [17],  the  HT  Riccati  solutions  are  functions  of  a  parameter  y  which 
is  an  upper  bound  of  the  optimal  HT  norm  of  the  closed-loop  system.  To  construct  a  better 
suboptimal  or  an  optimal  HT  controller,  it  is  necessary  to  search  for  the  optimal  H“  norm  of 
the  closed-loop  system,  i.e.,  the  smallest  y  such  that  the  two  HT  Riccati  solutions  X«(y) 
and  Y„(y)  are  positive  semidefinite  and  p[X«(y)Y«(y)],  die  spectral  radius  of  X»( y)Y„(Y), 
is  less  than  y2. 

Recently,  efficient  algorithms  for  computing  the  optimal  H*°  norm  were  proposed 
by  Scherer  [15]  and  Li  and  Chang  [22].  Scherer  considered  the  inverse  (or  pseudo  inverse) 
of  the  DGKF  H**  Riccati  solutions,  X»( y)  and  Y«(y),  defined  a  new  independent  variable  p 
=  Y2,  and  showed  that  these  inverses  are  concave  functions  of  p  in  matrix  sense  on  their 
domains  of  definition.  Li  and  Chang  showed  that  p[X«(y)Y„(y)],  the  spectral  radius  of 
X»(Y)Y„(y),  is  a  monotonically  nonincreasing  and  convex  function  of  y.  Based  on  these 
facts,  quacLatically  convergent  Newton-like  algorithms  were  proposed  to  compute  the 
optimal  HT  norm  [15,16]. 

Though,  these  monotonicity,  concavity  or  convexity  properties  have  only  been 
proved  for  the  DGKF  case  in  which  the  Dn  matrix  of  the  generalized  plant  is  assumed 
zero.  This  assumption  can  hardly  be  satisfied  by  many  practical  problems.  In  [8],  Glover 
and  Doyle  (abbr.:  GD)  considered  a  more  general  case  with  Dn  nonzero.  The  basic 
concept  is  the  same  as  DGKFs,  but  the  H“  Riccati  equations  involved  are  much  more 
complicated  which  makes  the  investigation  of  the  properties  of  the  GD  IT0  Riccati  solutions 
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extremely  difficult 

For  GD  FT*  control  problem,  Pandey  et  al/s  hybrid  gradient-bisection  method  [13] 
and  Chang  et.  al.'s  double  secant  and  bisection  method  [14]  were  proposed  for  the 
computation  of  the  optimal  H~  norm.  The  significance  of  the  conjecture  that 
p[X«(y)Y«(y)],  the  spectral  radius  of  X«(y)Y..(y),  is  a  convex  function  of  y  was  mentioned 
in  these  two  papers.  Since  there  was  no  proof  for  this  conjecture,  bisection  was  used  in 
these  two  algorithms  as  supplement  to  guarantee  convergence. 

In  this  chapter,  we  will  show  that  the  GD  H~  Riccati  solutions  have  the  same 
properties  possessed  by  the  DGKF  H°°  Riccati  solutions  as  mentioned  by  Li  and  Chang 
[16].  In  other  words,  the  GD  H~  Riccati  solutions  X»(y)  and  Y„(y)  are  nonincreasing  and 
convex  functions  of  y  in  the  domain  of  interest  and  so  is  the  spectral  radius  p[X»(y)Y„.(y)]. 
Based  on  these  properties,  a  quadratically  convergent  algorithm  is  proposed  to  compute  the 
optimal  HT  norm  for  die  general  H*°  control  problem  considered  by  GD. 

Section  4.2  lists  the  notations  used  in  this  chapter  and  reviews  the  fundamentals  of 
the  Riccati  equation.  The  main  result  of  the  chapter  is  presented  in  Section  4.3  which 
shows  the  details  of  the  proof  of  the  monotonicity  and  convexity  properties  of  X«,(y)  and 
Y.(y)  in  GD  case.  An  efficient  algorithm  to  computing  the  optimal  H°°  norm  based  on 
these  properties  will  be  given  in  Chapter  5. 


4.2  Notations  and  Basics  of  the  Riccati  Equation 


R 

The  set  of  real  numbers. 

R+ 

The  set  of  positive  real  numbers. 

j^mxn 

The  set  of  real  mxn  matrices. 

inf(S) 

The  infimum  of  the  set  S. 

sup(S) 

The  supremum  of  the  set  S. 

<cmxn 

The  set  of  nxm  matrices  whose  entries  are  complex  numbers. 

I 

Identity  matrix. 

In 

n-dimensional  identity  matrix. 
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0 


Zero  or  zero  matrix. 


0« 

At 

A* 

A*0 
A  >0 
MA) 
HA) 

^tn*x(A) 

XrnmCA) 

P(A) 

A(^ 

HL 

J- 

X-(y) 

Y«W 

P 

Pi 


H(H) 


n-dimensional  zero  matrix. 

Transpose  of  a  matrix  A. 

Conjugate  transpose  of  a  matrix  A. 

A  is  a  positive  semi-definite  matrix. 

A  is  a  positive  definite  matrix. 

An  arbitrary  eigenvalue  of  A. 

The  ith  largest  eigenvalue  of  A 
The  maximal  eigenvalue  of  A. 

The  minimal  eigenvalue  of  A. 

The  spectral  radius  of  A. 

The  maximal  singular  value  of  A. 

The  matrix  A  is  a  function  of  y. 

The  Hamiltonian  matrix  shown  in  (2-Sa). 

The  Hamiltonian  matrix  shown  in  (2- 5b). 

The  stabilizing  solution  to  the  Riccati  equation  corresponding  to  H~. 
The  stabilizing  solution  to  the  Riccati  equation  corresponding  to  J«. 
The  inverse  of  y  square. 

Pi  ;=  -r^ — ,  where  Xj  is  the  ith  eigenvalue  of  DnT  DJDj1  D^. 

■f-S 

Defined  as  Hl|^_^-i/2. 
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m 


Defined  as  JL.I 


X(p) 

9X 

dUi 

a2x 

YOi) 

X(Y) 

Y(Y) 

ax 

ay 

Px 

Py 

A I  B' 
C  I  D 

a 


The  second  older  partial  derivative  of  X  with  respect  to  m  and  fij,  i.e., 
•  • 

X  . 

The  antistabilizing  solution  to  the  Riccati  equation  corresponding  to  die 
following  Hamiltonian  matrix 


Defined  as  X(p)  |  -2  . 

Defined  as  Y(p)l  -2  . 

r  1 

Inf  {Y:  Ye  R+,  y>  o(Di.TDii)  and X« exists  }. 

Inf  { Y :  Y  e  R+ ,  y  >  and  Y„  exists  }. 

Inf  {Y:  Ye  R+,  yZa*  andXCy)  is  positive  definite } . 

Inf  { Y :  Y  e  R+,  Y  ^  °y  and  Y(y)  is  positive  definite}. 

A  state  space  realization  of  the  transfer  matrix  D+C  (si  -  A)  B. 
max{ax,Oy}. 
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0 

P(Y> 
IIG<s)!L 
y  x  (G,K) 


m«{0x,  0y>- 

Is  defined  as  p[X_(y)Y..(Y)]  on  0,  +eo). 

The  HT  norm  of  G(s). 

The  lower  linear  fractional  transformation  which  stands  for  die  closed  loop 
transfer  matrix  from  w  to  z  as  shown  in  the  following 


Let  A,  Q,  R  be  real  n  by  n  matrices  with  Q  and  R  symmetric,  rhe  algebraic  Riccati 
equation  and  its  corresponding  Hamiltonian  matrix  are  shown  as  follows: 

AtX  +  XA  +  XRX-Q  =  0  (4-1) 


H  = 


A 

Q 


(4-2) 


Assume  H  has  no  eigenvalue  on  the  imaginary  axis.  Then  finding  a  basis  for  the  invariant 
subspace  corresponding  to  the  eigenvalues  of  H  in  LHP,  i.e.,  the  open  left  half  complex 
plane,  and  using  the  basis  vectors  to  form  a  matrix,  we  obtain 


X_  (  H )  :=  Im 


where  Xi,  X2  e  Rnxn.  V  the  two  subspaces 


X_  ( H ),  Im 


0 

I 


are  complementary,  i.e.,  Xj  is  invertible,  then  the  stabilizing  solution  to  the  Riccati 
equation  is  X  =  X2  Xf1.  X  is  uniquely  determined  by  H,  Le.,  there  is  an  operator  denoted 
Ric  which  maps  H  to  X.  Thus,  X  =  Ric(H)  and  the  domain  of  Ric  is  denoted  by  dom( 
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Ric  ).  Furthermore,  if  X  exists  and  is  nonsingular,  then  X2  is  also  invertible. 


Consider  the  system 


G(s)  = 


A  B 

c  dJ 

«C(sI-A)_1B  +  D 


with  (  A,  B  )  stabilizable.  The  following  two  fundamental  lemmas  of  the  Riccati  equation 
will  be  employed  in  the  proof  of  the  main  results  in  Section  4.3. 


Lemma  4.1  Suppose  H  has  the  form 


H  = 


A 

T 

-c  c 


(4-3) 


Then  the  stabilizing  solution  X  to  the  corresponding  Riccati  equation  exists  and  is 
nonnegative  semi-definite.  The  null  space  of  X,  i.e.,  ker  ( X  ),  is  a  subset  of  the  stable 
unobservable  subspace. 


Stable  unobservable  subspace  means  the  intersection  of  the  stable  invariant 
subspace  of  A  with  the  unobservable  subspace  of  (  A,  C  ).  Note  that  a  detectable  ( C,  -A ) 
implies  that  X  is  positive  definite.  On  the  other  hand,  if  ( C,  -  A )  is  not  detectable,  then  A 
has  stable  modes  that  are  not  observable  from  C.  Assume  that  there  is  a  similarity 
transformation  U  such  that 


'A_ 

c 

A. 

0 

V 

T 

U 

b" 

D 

U  = 

®2 

C1 

fa 

0 

D 

with  A22  stable  and  (Ci,  -An)  detectable,  let  Xi  be  the  stabilizing  solution  to  the  following 
Riccati  equation 

AnVXiAu-xiBXx.  +  cIci-°. 


then  it  can  be  shown  that  X  =  U 


0  t 

U  is  a  stabilizing  solution  for  the  Riccati 

0 


equation  corresponding  to  (4-3). 
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Lemma  4.2  Assume  the  stabilizing  solution  X  to  (4-1)  is  positive  definite  and  &  is  the 
inverse  of  X,  then  &  is  the  antistabilizing  solution  to 

A£  +  &AT-&Q&  +  R  =  0  (4-4) 

die  eigenvalues  of  A-&Q  are  in  RHP,  Le.,  the  open  right  half  complex  plane. 

Proof:  The  Hamiltonian  matrix  associated  with  eq.  (4-4)  is 


-Q 

-A.. 


It  is  easy  to  see  that 


0  I" 

!  0. 


H 


0  I 
I  0 


If  die  stable  invariant  subspace  of  H  is  Im 


then  that  of  -ft  is 


and  hence  Im 


is  the  antistable  invariant  subspace  of  ft.  That  is,  &  is  the 


antistabilizing  solution  to  (4-4). 


Q.E.D. 


43  Properties  of  the  y-Dependent  H°°  Riccati  Solutions 

The  Riccati  equations  associated  with  the  Hamiltonian  matrices  (2-5a)  and  (2-Sb) 
are  dual  to  each  other,  so  we  will  only  concentrate  on  one  of  than,  say,  eq.(2-5a).  To 
investigate  the  properties  of  Riccati  solution  of  eq.(2-5a),  we  first  assume  that  (Dj1  Ci ,  - 
A+BR"1  Di.t  CO  is  detectable.  This  assumption  will  be  removed  later  in  this  sectioa  With 
the  assumption  it  is  obvious  that  X»  =  Ric(R»)  is  invertible. 

Although  our  objective  is  to  show  that  X.  and  Y«  are  nonincreasing  convex 
functions  of  y,  we  will  consider  their  inverses  first  since  the  inverses,  as  shown  later,  are 
analytic  functions  of  y. 

In  eq.  (2-Sa),  the  Hamiltonian  matrix  H«,  depends  on  the  inverse  of  R  which  is  a 


36 


function  of  y.  The  parameter  y  appears  almost  everywhere  in  the  entries  of  R'1  and  the 
matrix  H.  could  be  a  messy  function  of  y.  To  avoid  this  messiness,  the  following  lemma 
is  employed  to  confine  the  influence  of  y  in  a  diagonal  matrix. 


Lemma  4.3  The  matrix  R  in  (2-5d)  can  be  expressed  as  -T  ACY)'1  Tt,  where  A(Y)  is  a 
diagonal  matrix. 

Proof: 

By  Di2T  D12  =  I,  R  can  be  decomposed  as 


R  = 


dt  D  -  Y2 1 

ti  ll  1  ml 
T 

D  D 

12  ll 


T 

D  D 
11  12 

I 


1  i  DnDi? 
ml  11  12 

0  I 


(Dr  D  DTD  -  Y2 1  )  0 

v  1111  11  *  mV 

0  I 


ml 


0 


D  D  I 
12  11 


(4-5) 


Define  E  = 


E,  0 
0  I 


where  Ei  is  the  orthogonal  matrix  in  the  singular  value  decomposition 


of  D_ltDh,  i.e.,  Dut  DjDjJDh  =  Ei  diag{  Xi,  A,2, ....  Xmi }  EiT,  and  let 


T:= 


I  ,  d,,d10 

ml  11  12 

0  I 


E, 


and 


A(Y) ,=  diag{|ii  ,  p.2  «  Hmi  » -1»... ,  -1}, 


where  \ii  :=  ~ —  i  =  1,  2,... ,  mi, 

■f-h 

then  it  is  trivial  to  show  that 

R  =  -TA(Y)1Tt.  Q.E.D. 

This  lemma  shows  that  the  effect  of  y  is  only  on  the  eigenvalues  of  R.  This  fact  plays 
a  key  role  in  reducing  the  complexity  of  the  proof  that  the  X«  is  a  convex  function  of  y. 

In  order  to  specify  the  domain  of  die  mapping  XJ$,  we  define 
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Ox « inf  {Y:  Ye  E+,y>3(D.j.TDn)  and  X„  exists  }. 
To  simplify  the  proof,  let  p.  —  y2  and  define 


Dl== 


-T 


-dt  D 
12  11 


B 

T 

CD 

li. 


Dn:= 


[D.  Dn] 


and 


M  := 


(Hi  -dtd  dtd  )4  0 

VHinl  11111V 


Then  the  Hamiltonian  matrix  H_  in  (2-Sa)  can  be  written  as 
(A+TjMT^)  TjMT^ 

’(CICi+T2MTI>  *(A+T1M1^)T 


H(p.)  := 


Consider  the  following  Riccati  equation 

X(A+TiM  T2t)t  +  (A+TiM  T2t)X  +  TiM  Tit  +  X(CiT  Cj+T2M  T2t)X  =  0 
whose  associated  Hamiltonian  matrix  is 


•[:  m:  :i 


(4-6) 

(4-7) 


(4-8) 


(4-9) 


It  is  easy  to  see  that  X(p.)  exists  only  on  the  interval  (0,  ax'2).  In  the  proof  of  the 
following  theorem,  the  vector  vec  X,  associated  with  the  matrix  X,  is  defined  as  follows. 

Definition  4.1:  With  each  matrix  X  =  [xy]  e  E nxn,  we  associate  the  vector  vec  Xe 
E""  defined  by 

VCC  X  [Xu,  Xq|,  Xj2j  ...  Xn2,  ...»Xln>  x„n]T  (4-10) 

Theorem  4.1:  On  (0,  a*'2).  X(jx)  is  analytic  and  satisfies 
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(4-11) 

i.e.,  X(ji)  is  an  analytic  nonincreasing  function. 

Proof: 

To  prove  the  analyticity  of  X,  we  define 

F(X4i)  ••=  X(A+TiM  T2t)t  +  (A+TiM  T2t)X  +  TjM  Tit  +  X(CiT  Ci+T2M  T2t)X 
where  Xe  X 

X  :=  {X  €  Rnxn  I X  is  positive  definite  and 

(A+TiM  T2'i>+X(Cit Ci+T2M  T2t)  is  antistable} 

From  the  above  definition,  we  have 
F(X+8X,n)-F(X,p) 

=  8X  [(A+TjM  T2t)+X(Cit  Q+T2M  T2t)]t 

+  [(A+TiM  T2THX(Cit  C!+T2M  T21)]  8X 
+  5X  (CiT  Q+T2M  T2t)  5X  (4-12) 

where  §X  is  a  small  perturbation  on  X.  Let 

A  :=  [(A+TiM  T2t)+X(C!T  Cj+T2M  T2t)], 
then  the  partial  Frechet  differential  [23]  of  F(., .)  with  respect  to  X  is 

8F(X,p)  =  dX  AT  +  A  3X  (4-13) 

and  it  can  be  regarded  as  a  linear  map.  The  map 

aXaaXAT  +  AdX  (4-14) 

is  nonsingular  because 

vec  8F(X,p)  =  vec  8X  AT + vec  A  dX 

or  equivalently, 

vec  dF(X,fi)  =  [I®  A  +  A®  I]  vec  dX  (4-15) 

where  ®  is  the  Kronecker  product  [I®  A  +  A®  I]  is  nonsingular  because  of  the 
antistabilizability  of  X.  Consequently,  dF(X,^)  and  dX  are  uniquely  determined  by  each 
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It  is  clear  that  F(., .)  is  an  analytic  function  which  maps  (X,  p)  to  F(X,  p)  on  the 
domain  X  *  L  where  L  =  (0,  a*'2).  From  the  implicit  function  theorem  [23]  and 
analytidty  of  F(., .),  there  exists  a  neighborhood,  (p-e,  p+e)  with  e  >  0,  such  that  the  map 
p  a  X  is  analytic.  Hence,  X  is  an  analytic  function  of  p.  Next,  we  will  show  that  X  is  a 
nonincreasing  function  of  p.  Let 

X  =~- 
V  d|i. 

i  » 


M 

u*  dp 


i  1,  2,...  ,  mi 


Differentiating  eq.  (4-8)  with  respect  to  Pi  gives 


X^A  +AX^ +(XT2  +  T1)M4.(XT2  +  T1)  =0  i  =  1,  2 . mj 

Recalling  that  M  =  RAF.1,  we  have 


(4-16) 


M(.»E-|^ET«Ediag<0,...,0, 1,0 . 0)ET 

*1  ‘  —  —  . 


i-th 


which  is  positive  semi-definite.  Since  X  is  the  antistabilizing  solution  to  (4-9)  i.e.,  A  is 
antistable,  eq.  (4-16)  implies  that 


X  .<£0 


and  therefore 


dX 

dp 


ml  . 

IX 

i=l 


d-ty ) 


£  0 


(4-17) 


Q.E.D. 


Based  on  Theorem  4.1,  it  is  known  that  the  second  order  or  higher  order  derivatives  of 
X(p)  with  respect  to  p  exist  on  (0,  a*'2).  To  prove  that  X(p)  is  a  concave  function  in 
Theorem  4.2,  we  need  the  following  lemma. 


Lemma  4.4:  Let 

G  :=  E  diag{p'Vi,  p-!p2,  p^pmi,  0, ...,  0}  ET 


(4-18) 


CXDICl+tlT2GT2  T2°2 

G2Tj  h'g’-u'g2 


then  on  (0,  a*'2)  the  matrix  W  is  positive  semi-definite. 
Proof:  G  can  be  rewritten  as 

G  =  diag{n-i(  n-il^  -  DutD ±Dj  Dn)  \  0 . 0}. 


Wi  :=  \rl  ( H'1!*!  -  DnT  D±DjDn  )*> 


,  X.  X-  X  , 

w„  =  E  diag{ - S'. - S'.  ... , - — 

O-XjM)  (1-X  |i)  (1-X  n) 


then  we  have 


Since 


n‘(  G3-  G2)  =  diagfWu,  0 . 0). 


T2G2  “  [CitDxDj.tD„W^  0], 


W  can  be  written  as 


w.  °1 

W:= 

0  OJ 


where 


T  T  T  T  T  2' 

CDD  C+UT„GT  CD  D  D  W 
1  ill  ^2  2  1  1 1  11I 


w?dt  d  dtc 

I  11  1  11 


It  is  obvious  that  W  is  positive  semi-definite  if  and  only  if  so  is  Ws.  In  the  following,  we 
will  show  that  Ws  is  positive  semi-definite.  W5  can  be  decomposed  as 
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w  JCX  ol[I  +  t*I>>„W.D>x  D>X][D^Cl  o' 

'  U  Wn  J[o  I.. 

The  matrix  Q  ;=  I  +  (iDj_tDhWi  DnT  Dj.  is  positive  definite  since  Wj  is.  The  middle 
matrix  on  the  right  hand  side  of  the  above  equation  cun  be  further  decomposed  as 


Now,  the  proof  boils  down  to  checking  that  Wn  -  W^Dn1  PiQ^Djl1  DijWi2  is  positive 
semi-definite.  Based  on  the  facts  shown  as  follows. 


cr1 = a + hdJ  DnWi  dut  d jj*1 

= r1  -  i‘d1t  Du[  DiiTD1rlDiTDn+  4iwL)'1]‘1D,1TD1r1 
=  1  -  Dxt  D„[  D„TD1DiTD„+  h-1Iov1-D11tD1D1tD11)'1)1]-1 
•DutDx 

=1  -  DxT  Dn[  D„TDxDiTDI1+  (i-1!*,  -  D„tDxDxt  DuJ’Dh1  Dx 
=  I-|i-1DxtDiiDiitDx 
and 

DnT  Dj.Q4Dj.t  D„  =  Dut  D±(l  -  \rlD±T  DuDnT  D_l)  DxTDn 

=  DnT  DjPiT  Du-  (DnT  DiDit  Du)  H  (PnT  DjD±TDii) 
=  Ei  diag{Xi(l-Xi|x),  X2(1-A.2|a),  ^mi(l-^mi|A)}EiT, 

it  is  easy  to  show 

Wn  -  W^Dn1  DiQ-toJ  DnWi2 


=  Ei  diag{ - *-j - l—  Ul-Xn) - 1—1, 

d-XjH)3  o-^w  1  1  (l-x^r 
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1 


1 


is  a  zero  matrix  and  hence 


dtd  w2 

i  u  I 


is  positive  semi-definite.  Consequently,  Ws  is  positive  semi-definite  and  so  is  W.  Q.E.D. 

Based  on  the  fact  that  X(p)  is  an  analytic  nonincreasing  function,  as  shown  in  Theorem 
4.1,  the  second  order  derivative  of  XQi)  exists.  Furthermore,  we  also  found  that  X(p)  has 
the  following  important  property. 


Theorem  4.2:  On  (0,  a*'2),  X(p)  is  a  concave  function,  that  is, 
dp 

Proof: 

Taking  partial  derivative  of  (4-16)  with  respect  to  pj,  we  have 


(4-19) 


X  .  .  A  +AX..  +  X  T.M  (XT.+T,)1  +X  (C,  C,+T,MT*  )X  . 

»UW  M4U  ju  2  |Xjx  2  V  juv  1  1  2  2  '  MJ 

+  (XT  +T  )M  .T*  X  .  +  X  (CT  C,  +  TM  tI  )X  .  +  X  T.  M  .(XT+T)1 
2  r  |ij  2  |U  Wv  1  1  2  2 7  ju  |ij  2  |u  2  1 

+  (XT.  +T)M  T*X  .  =  0 

v  2  1  (U  2  )Xj 

t—l»  2,  ...,  mi,  j—  1,  2, ...,  mi 
For  notation  simplicity,  we  define 

x.-ss(x..  1  2  1  2) 

^  (1-X.p)2  (1-X.p)2 


(4-20) 


.  ® 1  . 
XT  =  IX  . 
*  itl  m 


(1-ty) 
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From  eq.  (4-18),  it  is  obvious  that 


gk-Im  — U- 


Now,  the  sum 
ml  ml 

IK- 


1 


(4-20)) 


can  be  simplified  as 


XjAT + a3Cj  +  2xit2g2(xt2+t1)t  +  2(XT2+ t^2  t2  Xj 

+  2XI(C’[C1+T2MT2)XI  =  0 


(4-21) 


(4-22) 


Define 


.  ml  .  2X. 


then  the  sum 


mi  2X. 

X - '—(.4-16) 

<1-X.fi) 


leads  to 


j  ml  #  2 X. 

V  +  ♦  (XT^)  (I M^^J)  (XT2+Ti)  =0 


(4-23) 


Adding  eq.  (4-23)  to  eq.  (4-22)  gives 

(X,+3^)AT+  AOCj+Xj,)  +  2X,T2G2(XT2+  T,)  +  2(XT2+ 1  ,>0^ 

+2Xj(CjCj+  T2Mt]  )Xj+  2(XT2+Tt)(  XXTf  -  0 


(4-24) 
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Based  on  the  fact  that 

Dj.  Di  *  Du  -  D12  Di2T  Dn  =  Dj.DitDh  ,  Di.  Dn  =  D12 


and 

M  =  |iG  +  diag{Omi,  -1} 
where  Omi  is  the  mi  by  mi  zero  matrix,  we  have 
T2  *  [  CitDjP±tDii  CjDd 


and 

CiTCi  +  T2  diag{Omi,  -1}  T2t  =  Ci^DjCj. 

Then  eq.  (4-24)  becomes 

(Xj+Xjj^+A^+Xjj)  +  2XJ2G2QCT2+  T,)  +  2(XT2+  T^T^X, 
+2XI(cX^C1)  X,+  T2Xj 

+  2(XT.+T  )(  I M - i-j  )(XT  +  T  )T  -  0 

2  1  i-1  l"(l-X.tl)  2  1 


With 


ml  X.  ml  .  1 

£m. — l-t  =  Xm.( — L  . 

'“(l-X.n)3  “>  "  (1-V4)  (1-X41) 


3-  i-j)!!'1  =(l1(G3-G2), 

*  /•  A  \™ 


eq.  (4-25)  becomes 

(Xj+X^A1*  AQCj+Xjj)  +  2XiT2G2(XT2+  Tj)  +  2(XT2+ 
+7X1(c]d±d]ci)  Xj+  2hXjT2G 
+  2(XT2+T1)(^'1G3-  41G2)(XT2+  Tj)1  =  0 

This  suggests  that 


(4-25) 
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X 

(X,+Xj,)S  +S(X,+Xn)  + 

-  2(XT2+T1)(m4G2)(XT2+T1)T 

T 

+  2[xiT2+h-1(XT2+T1)g]  (iiG)  [ XjT2+^  1(XT2+T1)g]  =  0 


Note  that 


2X,(CX  D^C^Xj  -  2(XT2+T1)(h-1G2)(XT2+T1)T 

T 

+  2[xtT2+^*1(XT2+T1)g]  (mG)  [XjTjVCXTj+TjKs] 
=  Vdiag{D11TD1DxTD11,-^,HG}  VT, 

where 


V»[Xj  ^(CTj+Tj)  XjT2+  h*1(XT2+  Tj)g] 


Now, 


V  diag{DnT  Dx  DjjDn,  -pi,  MG}  Vt 


CIDXCi+jiT2GT2  T2G 

_2  T  -13 

G  T2  G  -  M  G 


=  VWVT 


where 

V=[Xj  ^‘‘(XTj+TjJg] 

Since  W  is  positive  semi-definite  from  Lemma  4.4, 

V  diag{DnT  Dj_  D±T Dn,  -pi,  MG}  Vt 
is  also  positive  semi-definite. 

In  (4-26);  antistable  A  and  positive  semi-definite  V  diag{DnT  Dx  DXT  Du 
imply  that 


(4-26) 


,-mi,mg}  Vt 
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(Xj+Xjj^O 


Based  on  die  fact 


^4=XI(X  .  . - r - ±— 

dM  -U->  "'(l-Vli)2  (i-y> 

and  (4-27),  it  is  clear  that 

d2X  _ 


ml  ,  2X.  ,,  , 

^?.x*7rf7=x"+xn 

i  1=1  (1-X.|i) 


(4-27) 


Q.E.D. 


Define 

X(T)  »  X(n)l4=y-2 

then  we  have  the  following  corollary  which  is  a  direct  result  from  Theorem  4.1  and 
Theorem  4.2. 

Corollary  4.1:  X(y)  is  a  well-defined  function  on  (a*.  +oo).  Moreover,  X(y)  is  analytic 
and  satisfies 


dX(Y) 


(4-28) 


d2X(y) 


(4-29) 


on  (ax,  +oo). 

Proof: 

Since  X(y)  is  a  compound  function  consisting  of  analytic  X(n)  in  (4-9)  and  |i=y'2, 
X(y)  is  analytic  on  the  interval  with  nonzero  y.  With  the  chain  rule  of  differentiation,  we 


dX(Y)  dXQt)  dji  .3  dXQi) 
dy  ”  dji  dy  Y  d^i 
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Q.E.D. 


d2X(Y) 

d»2 


d^Oi)  <4i  2  dX(ji)  d2|i 

dn2  W  +  *  dr2 


which  imply  (4-28)  and  (4-29)  respectively. 


Based  on  Corollary  4.1,  we  have  the  following  theorem. 

Theorem  4.3:  On  (ax,  +©o),  X(y)  is  invertible  almost  everywhere. 

Proof: 

Since  X(y)  is  an  analytic  function,  it  is  well-known  that  .  :  eigenvalues  and 
eigenvectors  of  X(y)  are  analytic  functions  of  y.  Consider 
(a-X(y))v*0 


where  X  is  an  eigenvalue  of  X(y)  and  v  is  a  corresponding  eigenvector.  Then  it  is  easy  to 
have 

(U-i)v  +  ca-X)v  =  0  (4-30) 

which  implies 


vTX  v 


£0, 


(4-31) 


Let  Xj(y)  be  the  ith  eigenvalue  of  X(y),  i  =  1,  2, ....  n,  where  n  is  the  dimension  of  X(y). 
If  there  is  a  yi  e  (a*,  +oo)  such  that  Xj(yi)  =  0,  then  define  a  new  function  on  (a*,  +oo)  as 
following 

f(y)  :=  vi(y1)T  X(y)  Vj(yi)  (4-32) 

where  Vi(y0  is  an  eigenvector  corresponding  to  Xj(yi).  It  is  trivial  to  show  that  f(y)  is  a 
nondecreasing  concave  function  on  (a*,  +©o).  Since  f(+©o)  >  0  because  X(+oo)  >  0  under 
the  assumption  we  made  at  die  beginning  of  this  section,  we  can  conclude  dial 

df(Y.) 

f(Yj)  =  — dty  >  °*  (4-33) 

Based  oneq.  (4-31)  and  eq.  (4-32),  eq.  (4-33)  implies 
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X.(Y1)>0.  (4-34) 

By  the  property  that  X.fy)  >  0,  with  Yi  e  (ax,  +oo)  satisfying  X;(Yi)  =  0,  it  comes  the 
conclusion  that  X(y)  has  at  most  n  singularity  points  on  (a*,  +©o).  In  other  words,  X(y)  is 
invertible  almost  everywhere  on  (a*,  +©o).  Q.E.D. 

In  the  same  way,  define 

<Xy  :=  inf  { Y :  Y  €  R+ ,  y  >  G(DjDnT)  and  Y„  exists  } 

and  let  Y  be  an  antistabilizing  solution  to  the  Riccati  equation  corresponding  to  the 
following  Hamiltonian 

-[: ;]«[:  •■]. 

Then  we  have  the  following  theorem. 

Theorem  4.4:  On  (0,  a/2),  Y(ji)  is  an  analytic  function  which  satisfies 

(4-35) 

(4-36) 

Y(y)  =  Y(n)|^-2 

then  we  have  followi, corollary. 


and 


dY 

dp 


<0 


dp2 


<0 


Define 


Corollary  4.2:  Yfy)  is  a  well-defined  function  on  (cty,  +oo).  Moreover,  Y(y)  is  analytic 
and  satisfies 


dY(Y) 

dY 


£0 


(4-37) 


(4-38) 


£0 


on  (Cty,  +©o). 

Based  on  the  above  corollary,  we  have  the  following. 

Theorem  4.5:  On  (cty,  +oo),  Y(y)  is  invertible  almost  everywhere. 

Now,  we  are  ready  to  show  that  X«(y)  and  Y»(y)  are  nonincreasing  convex 
functions  in  the  domain  of  interest  Without  loss  of  generality,  the  assumptions  that  (Dj_T 

~T 

Ci ,  -  A+B2Di2tCi)  is  detectable  and  (-  A+BiD2iTC2,  BiB  )  is  stabilizable  are  still  used 
in  the  following  to  simplify  the  proof. 

From  eq.  (4-9),  it  is  clear  that  if  X(y)  =  XQi)  is  invertible  then  X«(y)  =  Xty)'1. 
Although,  X(y)  is  not  always  invertible  on  (ax,  +oo),  from  Theorem  4.3,  it  can  only  lose 
rank  at  some  specific  isolated  points.  Since  X(y)  is  a  nondecreasing  concave  function  as 
mentioned  in  Corollary  4.1,  there  must  exist  &  such  that  X(y)  is  positive  definite  when  y  > 
P*.  That  is, 

Px  :=  inf  {y :  y  e  E+,  y  £  a*  and  X(y)  is  positive  definite}. 

This  means,  on  the  interval  (Px,  +oo),  X(y)  is  positive  definite. 


Theorem  4.6:  On  (Px,  +oo),  X.  is  an  analytic  nonincreasing  convex  function  of  y,  i.e.. 


and 


(4-39) 

(4-40) 


Proof: 


Because  X(y)  is  positive  definite  and  X.  =  Xfy)'1,  by  the  chain  rule  of 
differentiation  and  Corollary  4.1,  we  have 


dy  dy  ~dy  “ 


(4-41) 
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and 


A-  d(X-f  X->  x  d*  A  _ 

^  dr  ■'  (tr  ■x-(2drx-dr'^)x- 


(4-42) 


Since  X.  is  positive  definite  and  T  is  negative  semi-definitive,  it  turns  out  that 

<*Y 


V  A. 


Q.E.D. 


Similariy,  by  defining 

Py  :=  inf  {Y:  Ye  E+,  Y^  Qy  and  Y(y)  is  positive  definite} 
we  have  die  following  theorem  for  Y* 

Theorem  4.7:  On  (Py,  +©o),  Y„  is  an  analytic  nonincreasing  convex  function  of  y,  i.e.. 


0 

VI 

(4-43) 

..  d\ 

Y-=  J"*0 

(4-44) 

dy 

At  the  beginning  of  this  section,  we  assumed  that  (DjJ  Ci ,  -  A+BR*1  Di.t  Ci)  is 
detectable,  or  equivalendy  that  (DxT  Ci ,  -  A+B2Di2TCj)  is  detectable.  This  assumption 
can  be  removed  by  the  following  arrangement  If  (DjJ  Ci ,  -  A)  is  not  detectable,  where 

A  =  A+BDi2T  Ci, 

one  can  always  find  an  orthogonal  matrix  U  =  [Ui  U2]  such  that 

U^AUj  0 

u>,  , 

D±tCiU  =  [DxtCiU1  0], 

with  (DxTCiUi,  -  UitAUi)  detectable  and  X*  =  Ric(  H» )  can  be  expressed  as 
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x_=u 


Note  that  U  is  not  a  function  of  y,  and  Xi  is  the  Riccati  solution  of  eq.  (2-5a)  with  (A,  B, 

Ci)  replaced  by  (UitAUi,  UitB,  C1U1).  Therefore,  no  matter  whether  (Dj.t  Ci ,  -  A)  is 
detectable  or  not,  it  is  always  true  that  XJy)  exists  almost  everywhere  on  (a*,  +©o)  or  Xi 
is  nonsingular  almost  everywhere  on  (a,,  +00).  Furthermore,  Xi  is  positive  definite  on 
(Px.  +•»). 


Define 

a  :=  max  {a*,  cty}. 


and 


P  :=  max  {px,py}, 

then  X»  and  Y„  exist  on  (a, +00)  and  are  positive  semi-definite  (P,  +00). 


The  following  theorem,  which  is  one  of  the  main  results  in  Li  and  Chang  [16], 
showed  that  if  X.  and  Y„  are  nonincreasing  convex  functions  of  y,  then  so  is  p[X«Y„]. 

Theorem  4.8:  Let  P  be  the  infimum  of  y  such  that  both  X»  and  Y„  exist  and  are  positive 
semi-definite.  On(P, +00), 

a)  all  eigenvalues  of  X»Y«,  are  smooth,  nonincreasing  functions  of  y ; 

b)  p(X« Y«.)  is  a  nonincreasing,  convex  function  of  y. 

Proof: 


a)  It  is  quite  straightforward  to  show  that  for  any  nontrivial  eigenvalue  of  X»YW  we  have 


X[x_<y)y_(y)j  =  X[vtx_(y)Y-(y)v] = X(v‘u 


-  \  fVTT 


Xfi)  0 
0  0 


UTV 


Yj(Y)  0 
0  0 


=  X( 


V^X^U^Y,®  0 

0  0 


)  -yzwYjtif)], 


where  Z(y)  *  ZT(y)  =  has  following  properties 


Z(Y)^0, 

ky)zo 
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(4-45a) 

(4-45b) 


and 


2<Y»S0 

on  0,  -h»o).  Now  consider 

[IX-Z(y)Y1(y)]w  =  0, 


(4-45c) 

(4-46) 


where  X  is  any  eigenvalue  of  Z(y)Y1(Y)  and  w  a  corresponding  eigenvector.  Note  that 
this  equation  is  equivalent  to 

[XY1'l(K)-Z(1ft)1i  =  0.  (447a) 

or 

[Xt(Y)*Z(Y)]u  =  0  (4-47b) 

with  %)  =  Yj'^y)  and  u  =  YjW.  From  (4-47),  we  have 


and 


(X$  +  X*-Z)u  +  (Xt-Z)u=0 
uT(xt  +  Xt-Z)u  =  0. 


Then  it  is  easy  to  show  from  these  two  equations  that 
• ,  ,  uT(Z-Xt)u  ^  „ 

x(Y)=-^r'° 


(4-48) 

(4-49) 

(4-50) 


on  (P,  +o»),  since  Z  <  0,  $  £  0  and  $  k  0  (see  (4-45b)  and  (4-37)  in  Corollary  4.2). 


b)  Next  is  die  proof  of  the  convexity  of  p(X»Y„).  From  the  equations  (4-48)  and  (4-49) 
with  X  replaced  by  the  maximum  eigenvalue  of  X.Y.,  we  come  to  the  following  two 
equations: 


*T 

U 


*T 

U 


(  X|jj*x^  +  Xmtx^  -  Z  )  U  =  -UT  (  Xmax^  -  Z)  U 
(  Xmax^  +  Xjnax^  -  Z  )  U  +  U^  (  Xjj^^  +  Xmlx^  -  Z  )  U 

+  «T(*im«^  +  2Xmax^  +  ^n«?  -Z  )  U  =  0. 


(4-51) 

(4-52) 


Combining  die  above  two  equations,  we  have 


^maxCy)  = 


2uT(Xmaxt-Z)u4uT(^2Xm»ixiptx?)u 

uT$u 


(4-53) 
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Note  that  £  0,  0,  ^  £  0,  0  and  ^  £  0,  we  have  X^Cy)  ^  0  and 

therefore  p  PC.(y)Y..(y)]  £  0  on  (P,  +«»).  Hence  we  complete  the  proof  that  p(X«Y»)  is  a 
convex  function  of  y  on  (P,  +oo).  Q.E.D. 

The  properties  of  the  H“  Riccati  solutions  presented  in  this  section  will  be 
employed  in  Chapter  5  to  develop  efficient  algorithms  for  computing  the  optimal  H“  norm. 
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CHAPTER  5 


COMPUTATION  OF  THE  OPTIMAL  H-  NORM 


Recently,  an  efficient  algorithm  for  computing  the  optimal  H“  norm  was  proposed 
by  Scherer  [15].  Scherer  considered  the  inverse  (or  pseudo  inverse)  of  the  DGKF  H** 
Riccati  solutions,  X«(y)  and  Y„(y),  defined  a  new  independent  variable  p  =  y2,  and 
showed  that  these  inverses  are  concave  functions  of  {!  in  the  matrix  sense  on  their  domains 
of  definition.  Based  on  this  fact,  a  quadratically  convergent  Newton-like  algorithm  was 
proposed  to  compute  the  optimal  H“  norm. 

Pandey  et  al.'s  hybrid  gradient-bisection  method  [13]  and  Chang  et.  al.'s  double 
secant  and  bisection  method  [14]  were  also  proposed  for  the  computation  of  the  optimal  H** 
norm  The  significance  of  the  conjecture  that  p(y)  :=  p[X«(y)Y..(y)],  the  spectral  radius  of 
X-(Y) Y«(y),  is  a  convex  function  of  y2  was  mentioned  in  these  two  papers.  Since  there 
was  no  proof  for  this  conjecture,  bisection  was  used  in  these  two  algorithms  as  supplement 
to  guarantee  convergence. 


Based  on  the  properties  of  the  Riccati  solutions  described  in  the  previous  chapter,  a 
quadratically  convergent  algorithm  can  be  easily  developed  to  compute  the  optimal  H“ 
norm  According  to  Theorem  2.1,  we  can  see  that  finding  the  optimal  HL  norm  denoted  as 
y„,  is  equivalent  to  finding  the  infimum  y  such  that  all  four  conditions  in  (2-6)  hold.  From 
the  previous  chapter,  it  is  obvious  that  yje  [P,  +°°).  It  is  possible  for  P  to  be  y„,  especially 
when  p  and  a  are  identical,  however,  with  very  few  exceptions,  y„e  (P,  +«>),  which 
implies  that  y„  is  the  solution  to  p(Y)  =  y2.  The  relations  between  a,  P  and  yH  are  shown  in 


the  figure  below. 


Both  X~  £  0  and  Y-  >  0 _ 

All  four  conditions  hold 

0  OL  | 

p  H 

l 

Fig.  5.1 


5.1  Computation  of  y~ 


Fig.  5.1  implies  that  the  problem  of  finding  the  optimal  y„  is  actually  that  of  either 
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searching  for  die  intersection  point  of  p(y)  with  y2  inside  ({},  +«)  or  confuting  the 
boundary  point  0.  The  monotonicity  and  convexity  of  p(y)  in  Theorem  4.8  suggests  that  if 
we  use  a  gradient  method  to  search  for  yM  with  a  starting  point  yn  inside  the  interval  0, 
yj,  then  the  convergence  is  guaranteed. 

The  derivative  of  p(y)  at  y„  can  be  computed  as  follows, 

(1.) 

V  w 

where  v  and  w  are  left  and  right  eigenvectors  of  X„Y..  respectively  corresponding  to  its 
maximal  eigenvalue.  and  Y„  can  be  obtained  by  solving  the  following  Lyapunov 
equations: 

aJL+  X„V  (XMB+CjDlt)  R_1R  R^B+C^/  =  0  (5.2a) 

AyY„+  Y_aJ+  (CY„+  D  R^CY^  D  ^  =  0  (5.2b) 

with 

Ax=  A  -  BR^DjJCi  -  BR'1BtX« 
and 

Ay=  A  -  BiD.!T  R  JC  -  CtR  CY» 


Assume  that  we  have  a  starting  point  yn  in  the  interval  (($,  yj,  then  the  optima1  y 
can  be  obtained  easily  as  follows.  Refer  to  Fig.  5.2,  draw  the  tangent  line  with  slope  p  (yn), 
which  can  be  computed  by  equations  (5-1),  passing  through  the  point  (y„,  p(yn)).  The 
abscissa,  yn+1,  of  the  intersection  of  the  tangent  line  and  the  curve  y  =  y2,  always  lies 
between  y„  and  y„.  The  search  process  is  repeated  until  the  gap  y„  -  yn+i  is  small  enough. 
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p  \  \ 


Furthermore,  we  will  see  that  the  convergence  rate  is  quadratic.  Define  e„  =  y„  -  yn  and 
e„+1  =  y„-  Yn+i-  It  is  straightforward  to  show  that 

pfrJ  2 

6,1+1  “  2  I  p(yj-2yj  66  (5-3) 

which  implies  quadratic  convergence. 

In  the  above  algorithm  of  computing  we  assumed  that  a  yn  inside  the  interval  (P, 
Y«)  is  available  to  start  with.  To  find  a  yn  inside  this  interval,  usually  we  arbitrarily  pick  up 
a  relatively  large  Yi  at  which  the  two  Riccati  equations  (2-5a  and  b)  have  positive 
semidefinite  solutions  X«,(Yi)  and  Ym(Yi).  Refer  to  Fig.  5.3,  draw  a  line  passing  through 
the  point  (yx,  p(Yi»  with  slope  p(Yi)-  The  abscissa,  y2*  of  die  intersection  of  the  straight 
line  and  the  curve  y  =  y2,  is  always  less  than  If  y2  >  {3,  i.e.,  the  two  Riccati  solutions 
X.( Y2)  and  Y^Yz)  exist  and  are  positive  semidefinite,  then  we  are  ready  to  use  the  above 
quadratically  convergent  algorithm  to  compute  yM. 


Most  of  the  time,  the  method  described  in  the  previous  paragraph  gives  a  yD  inside 
the  interval  0,  ym)  without  the  computation  of  {5.  However,  this  method  may  fail.  If  a 
starting  Yi  gives  a  y2  <  P,  i.e.,  either  X~(y2)  or  YM(y2)  does  not  exist  or  is  not  positive 
semidefinite,  one  may  suggest  that  a  smaller  Yi  could  produce  a  y2  >  p.  For  the  case 
described  by  Fig.  5.4,  a  smaller  Yi  could  do  the  job.  However,  for  the  case  of  Fig.  5.5,  y2 
is  always  less  than  a,  and  therefore  less  than  P,  no  matter  how  small  Yi  is.  Since  it  is 
difficult  to  tell  which  case  we  are  facing  and  there  is  no  efficient  guideline  to  reduce  Yi,  we 
suggest  to  compute  $  (or  a  and  P)  if  one  or  two  trials  of  Yi  does  not  give  a  y2  >  p.  If 
either  X^(Y2>  or  Y„(y2)  does  not  exist,  we  will  compute  a  first,  then  compute  P  if 
necessary.  The  computation  of  a  will  be  given  later  in  Section  5.3.  If  both  and 

Y„(Y2)  exist  but  are  not  all  positive  semidefinite,  then  y2  is  inside  the  interval  (a,  P).  A 
quadratically  convergent  algorithm  for  the  computation  of  P  is  given  in  Section  5.2.  Once  p 
is  obtained,  then  p+e  can  be  served  as  the  starting  point  Yn  in  Fig.  5.2,  where  e  is  a  very 
small  positive  real  number. 


52  Computation  of  P 

In  this  section,  we  have  the  assumption  of  having  a  starting  point  inside  the  interval 
(ax,  Px)  (resp.  (<Xy,  Py»,  which  implies  ax*Px  (resp.  ay*Py ).  We  will  show  how  to 
compute  a  in  the  next  section.  Recall  that  P  is  equal  to  max{  Px,  Py }  where  px  and  Py  are 
infimums  of  y  such  that  (2-6b)  and  (2-6c)  are  true  respectively.  Since  the  computation  of 
Py  is  similar  to  that  of  Px,  only  the  algorithm  for  Px  is  given  in  the  following. 

Without  loss  of  generality,  we  assume  that  (D_lt  Ci ,  -  A+B2Di2TCi)  is  detec^lue  and 
X(y)  ==  X(p)|  -2  where  X(p)  is  the  antistabilizing  solution  of  (4-9).  Define  ex(y)  := 

de  (y) 

Xmin(X)  and  e  ;=  — j- —  on  (ax,  -k»),  according  to  the  proof  of  Theorem  4.3  in  the 

*  (TV 


previous  chapter,  e,(Y)  is  a  nondecreasing  concave  function  of  y  on  (a,,  H  and  (J*  is  the 
y  such  that  e*(Y)  equals  to  zero.  Then  Newton's  search  scheme  can  be  used  to  compute  (5,. 


With  a  starting  point  which  is  slightly  greater  than  a*,  we  can  confute  P*  by  the 
following  iteration, 

Qn+l  Qn  CX(Y„) 

where 

TdX 
w  w 

w  w 

In  the  above,  w  is  an  eigenvector  of  X  corresponding  to  its  minimal  eigenvalue  Xmin(X) 
dX 

and  can  be  determined  by 

dX(Y)  dX(p)  dp  .3  dX(p) 

■  =  -2Y 


dy 


dp  dy 


and 


1 


dX(p)  ■?. 


dp 


£  0 


where  X^.  can  be  obtain  from  (4-16). 


Fig.  5.6  Newton  algorithm  for  searching 


As  we  assumed,  we  have  a  Yn  inside  the  interval  (a*,  px),  i.e.,  X(Yn)  exists  and 
e*(Yn)  <  0.  Refer  to  Fig.  5.6,  draw  a  line  with  slope  e*(Yn)  passing  through  the  point  (yn. 
This  straight  line  will  intersect  the  horizontal  line  y  =  0  at  Yn+i  which  always  lies 


59 


between  yn  and  Pv  The  search  process  is  repeated  until  the  gap  |Jx-Yn+1  is  small  enough. 
Hence,  the  convergence  is  guaranteed.  Furthermore,  we  will  see  that  the  convergence  rate 
is  quadratic.  Define  e„  =  fi*  *  y0  and  =  f^- yn+1.  It  is  straightforward  to  show  that 


®n+l 


'  ?»(P„)  I  ^2 
I  2  e„(P*)l  % 


which  implies  quadratic  convergence. 


(5-4) 


In  the  above  algorithm  of  computing  we  assumed  that  a  yn  inside  the  interval 
(«x.  Px)  is  available  to  start  with.  This  y„  could  be  obtained  in  Section  S.l,  i.e.  y2  i°  Fig 
5.3.  However,  if  either  XJyj)  or  Y„(Y2)  does  not  exist,  Le.,  if  y2  <  <*>  we  will  compute  a 
first 


53  Computation  of  a 

Recall  that  a*  (resp.  cty)  is  the  infimum  of  y  such  that  the  Hamiltonian  matrix  HL(y) 
in  (2-5a)  (resp.  JM(y)  in  (2-5b))  has  no  jto-axis  eigenvalues  and  a:=  max  {a,,  oty }.  As  we 
mentioned  before,  it  is  possible  for  a  to  be  the  optimal  y,  i.e.,  all  the  four  conditions  in  (2- 
6)  hold  at  this  point  If  this  is  not  the  case,  a  can  be  used  as  a  starting  point  to  search  for  (3 
in  Section  5.2  and  fJ+e  in  turn  can  serve  as  a  starting  point  to  search  for  in  Section  5.1, 
if  P  is  not  the  optimal  y.  In  the  following,  only  the  computation  of  a*  will  be  discussed, 
since  the  computation  of  cty  is  similar  to  that  of  a*. 

One  can  easily  have  a  bisection  iterative  algorithm  to  search  for  the  infimum  of  y 
such  that  the  Hamiltonian  matrix  H_(y)  has  no  joo-axis  eigenvalues.  This  method  is  slow 
and  therefore  not  recommended. 


To  have  a  more  efficient  algorithm,  we  need  the  fact  that  a*  can  be  also  expressed 
as  the  supremum  of  a  frequency  function.  According  to  [14], 

<&  *jp  { -  g2(g;g/g;  ig.oo.)} 

where  *  is  conjugate  transpose  and  Gi(s)  and  G2(s)  are  given  by 


[GjCs)  G2(s)]  = 


'A 

b. 

cT 

D 

D„ 

i 

li 

12 

By  defining 
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■*»>  -  w^1  -  g2(g;g2)'g;  ig,  oo.) 

the  problem  of  finding  a*  becomes  that  of  finding  the  supremum  of  the  function  h(to). 

There  are  several  efficient  algorithms  available  for  searching  for  the  supremum  of 
h(Q))  [24,25],  We  can  start  from  arbitrarily  choosing  a  frequency,  say  o^.  Let  y  =  h(o>j) 
and  then  find  all  the  positive  real  oa’s  such  that  h(co)  =  y.  These  CD's  can  be  easily  obtained 
from  computing  the  j to- axis  eigenvalues  of  the  Hamiltonian  matrix  H_(y).  Now,  we  have 
the  frequency  intervals  in  which  h((D)  £  y.  Evaluate  h(to)  for  each  midpoint  of  these 
frequency  intervals  and  update  y  to  be  the  maximum  of  these  h(co)'s.  Then,  find  the  new 
frequency  intervals  in  which  h(co)  2  y.  According  to  [25],  the  convergence  of  this  iterative 
process  is  quadratic.  This  process  can  be  repeated  until  only  one  frequency  interval  with 
h((o)  2:  y  is  left  and  the  interval  length  is  negligible  [25],  In  [24],  this  process  is  terminated 
when  a  frequency  interval  in  which  h(o))  is  convex  and  greater  than  y  is  found.  A  search 
method  called  Brent  method  was  used  to  search  for  the  supremum  of  h(o>)  in  the  convex 
frequency  interval. 

5.4  An  Illustrative  Example 

The  following  is  a  simple  H*°  optimization  problem  which  is  used  to  illustrate  the  proposed 
algorithm  of  computing  the  optimal  H“  norm.  A  realization  of  the  generalized  plant  G(s)  is 
given  by 


Starting  from  y0  =  100,  we  have  p(yo)  =  p[X»(y0)Y.,(y0)]  =  2.16e01.  The  slope  of 
p(y)  at  this  point  (y0,  p(Yo))  is  p  (Yo)  =  -3.49e-05.  The  tangent  line  at  this  point  will 
intersect  with  the  curve  y  =y2  at  y^  =  4.65.  Again,  evaluate  p(Yi)  =  2.24e01  and  compute 
the  slope  p  (y)  at  (yx,  pty)).  Since  pCyj)  >  7^,  y1  is  inside  the  interval  (fJ,Y„)  and  therefore 

from  now  on  the  convergence  is  guaranteed.  The  process  is  repeated  until  the  gap  between 
Vp(Y«»>  and  yn  is  small  enough.  The  following  data  show  that  only  four  iterations  are 
needed  to  reach  the  optimum,  7^=  4.734160476390407,  with  accuracy  better  than  10*14. 
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Iter 

Yn 

P(Yn> 

wm 

100 

9.535e+01 

-3.487e-05 

i 

4.647998761S38930e+00 

8.403e-01 

-3.811e-01 

2 

4.734064423866624c+00 

9.440c-04 

-3.595e-01 

3 

4.734 160476276923e+00 

1.115e-09 

-3.595e-01 

4 

4.734160476390407e+00 

7.105e-15 

By  the  formulas  in  (2-8)  and  the  descriptor-form  technique  described  in  Section 
2.2,  we  are  able  to  construct  an  optimal  controller  as  follows: 


Kopt(s)  = 


-.87542  -0.13925 

4.42042  -4.73416 


with  the  H“  norm  of  the  closed-loop  system  equals  y^.  Note  that  the  optimal  controller  has 
a  direct  feedthrough  term  and  thus  has  infinite  bandwidth.  If  we  choose  y =  4.8  which  is 
about  1.4%  higher  than  y^,  we  have  a  suboptimal  controller 


K,ub(s)  = 


’  -8.67072e-01  1.32928e-01 

-1.38959e-01 

-1.38320e+01  -1.52323e-K)2 

4.73733e+00 

-9.30025e+00  -1.49792e+02 

0 

which  has  a  reasonable  bandwidth  and  the  closed-loop  H“  norm,  IITzvIIod<  4.8  which  is  only 
1.4%  away  from  the  optimal  H~  norm. 


CHAPTER  6 


AN  APPLICATION  OF  p-SYNTHESIS 
TO  A  ROBUST  FLIGHT  CONTROL  PROBLEM 


6.1  Introduction 

In  general,  a  system  with  uncertainties  can  be  divided  into  two  parts,  the  nominal 
plant  M(s)  and  uncertainties  A.  If  A  is  norm  bounded  but  otherwise  unconstrained,  it  is 
called  unstructured  uncertainty.  However,  if  A  has  a  structure  constraint,  e.g.,  A  is 
diagonal,  then  it  is  called  structured  uncertainty.  Obviously,  the  set  of  unstructured 
uncertainties  is  larger  than  that  of  structured  uncertainties.  H~  design  technique  [1,2,7, 8] 
can  be  used  to  solve  robust  stability  problems  with  unstructured  uncertainties,  whereas  p- 
synthesis  technique  is  used  to  handle  systems  with  structured  uncertainties.  Since  the 
structure  information  of  the  uncertainty  is  used  in  p-synthesis,  a  less  conservative  solution 
can  be  obtained.  Moreover,  robust  performance  problem  can  also  be  formulated  as  a 
structured  uncertainty  problem  and  solved  by  p- synthesis. 

To  use  p-synthesis,  one  needs  to  formulate  die  problem  first  and  put  it  into  a  Linear 
Fractional  Transformation  (LFT)  [2]  form,  while  choosing  weighting  matrices  for  certain 
design  specifications.  With  this  augmented  plant,  D-K  iteration  design  algorithm  is 
employed  to  solve  the  p-synthesis  problem.  The  D-K  iteration  consists  of  H“  optimization 
step  (K-step)  and  p-analysis  step  (D-step).  While  each  step  is  a  convex  optimization 
problem  [26],  the  overall  optimization  problem  is  not  convex.  Hence  a  local  minimal  point 
it  converges  can  not  be  guaranteed  to  be  the  global  minimal  point  Another  comment  on  D- 
K  iteration  algorithm  is  that  to  have  a  better  accuracy,  a  higher  order  curve  fitting  in  D-step 
is  required,  which  in  turn  gives  a  higher  order  controller.  However,  the  high  order 
controller  usually  can  be  reduced  to  that  of  the  generalized  plant  without  degrading  the 
system  performance. 

The  chapter  is  organized  as  follows.  Section  6.2  gives  a  review  of  some 
preliminaries.  We  will  show  the  problem  formulation  in  section  6.3.  In  section  6.4  the  p- 
synthesis  is  used  to  design  a  robust  controller.  Simulations  and  analysis  are  also  included 
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in  this  section.  Section  6.5  is  a  conclusion. 

6.2  Preliminaries 

In  this  section,  we  will  give  a  quick  review  of  some  fundamental  knowledges 
including  the  definition  of  robust  stability,  and  robust  performance. 

Usually  the  mathematical  representation  of  a  system  under  consideration  is  only  an 
approximation  to  the  actual  system.  The  "actual  model"  M  is  in  a  neighborhood  of  the 
nominal  model  M.  According  to  [4],  M  can  always  be  separated  into  two  parts  illustrated 
by  the  following  block  diagram  with  an  M-A  structure,  where  M  is  the  nominal  system  and 
A  is  a  block  diagonal  uncertainty  matrix. 


Fig.  6.1  M-A  structure 


With  this  block  diagram,  the  stability  analysis  essentially  boils  down  to  ensuring  that  /  - 
MA(jco)  remains  nonsingular  at  all  frequencies  under  all  A  considered  (small  gain  theorem 
[1] ).  This  statement  can  be  seen  by  considering  the  following  fact  Assume  M(s)  is  a  stable 
system.  When  no  perturbation  is  introduced,  i.e,  A  =  0,  M  =  M.  However,  as  A  grows  up, 
in  the  sense  of  <j(A)  (2-norm),  /  -  MA  may  not  remain  nonsingular,  which  implies  the 
instability  of  M.  In  the  following,  we  will  introduce  the  definition  of  p,  which  measures 
how  large  A  can  be  such  that  M  remains  stable. 

The  structured  singular  value  (SSV  or  &i)  of  a  system  M  with  respect  to  the  given 
structure  of  A  is  defined  by 

|i(M)= -|nmi  (o(A)  :det(I  -  MA  )  =  o)}  V  0)€  R+  (6-1) 

where  R+  is  the  set  of  positive  real  numbers.  In  the  following  are  listed  several  remarks  on 
die  above  definition. 

*  Obviously,  the  smaller  the  |i  is,  the  larger  uncertainties  are  allowed  such  that  M 
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remains  stable. 


•  p  depends  highly  on  the  structure  of  A;  the  more  information  used,  the  less 
conservative  the  solution  will  be.  If  A  is  a  full  block  matrix,  p(M)  =  c(M).  On  the  other 
hand,  if  A  is  a  diagonal  matrix,  then  p(M)  =  p(M),  where  p  denotes  spectral  radius  [4]. 
Usually,  A  is  in  block  diagonal  form,  and  obviously  p(M)  £  p(M)  <,  o(M)  in  general. 

*  The  computation  of  p  is  still  an  open  problem.  If  A  has  three  or  less  blocks,  p  can 
be  computed  by 

P(M)  *  igffiftDMD'1)],  (6-2) 

where  D  is  a  positive  definite  block  diagonal  matrix  with  the  same  structure  as  A.  However, 
for  the  uncertainties  with  more  than  three  diagonal  blocks  the  right  hand  side  term  of  eq.(6- 
2)  is  just  an  upper  bound. 


Robust  stability  and  robust  performance 

Most  control  problems  can  be  represented  by  the  following  Linear  Fractional 
Transformation  (LFT)  form: 

P 
w 


Fig.  6.2 

where  G(s)  is  the  generalized  plant  including  possible  weighting  matrices,  z  is  the 
controlled  output  usually  including  the  error  signal  and  a  weighted  control  input,  w  is  the 
exogenous  input  containing  the  disturbances,  noises  and  commands,  and  y  is  the  measured 
output  vector  consisting  of  all  the  signals  which  can  be  measured  and  available  for 
feedback.  A!  represents  plant  uncertainties.  The  design  objective  is  to  find  a  controller  K(s) 
such  that  the  system  robust  stability  and  robust  performance  can  be  achieved,  i.e.,  under  all 
the  perturbations  considered,  the  system  remains  internally  stable  and  the  H°°  norm  of  the 
transfer  function  'rom  w  to  z  remains  less  than  a  prescribed  value.  This  problem  can  be 
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solved  by  ^-synthesis  technique  on  the  following  block  diagram. 


r~ 

:s 


Fig.6.3 

where  A  =  diag{  Alt  A2}  and  Aj  is  a  full  block  matrix  corresponding  to  the  transfer  function 

of  the  closed  loop  system  from  w  to  z.  If  we  combine  G(s)  and  K(s)  together  by  LFT,  then 
it  corresponds  to  the  M  shown  in  Fig.6.1,  i.e.,  M  =  Fj(G,  K). 

11-synthesis  and  D-K  iteration 

The  original  concept  of  ^-synthesis  is  to  design  a  controller  K(s)  such  that  |i  of  M  = 
Fj(G,  K)  is  minimized  at  all  frequencies,  with  respect  to  certain  A  structure.  According  to 

the  main  loop  theorem  [6],  robust  stability  and  robust  performance  can  be  achieved  if  the  (i- 
synthesis  is  applied  on  A  =  diag{A1,  A2}  of  the  above  block  diagram,  where  Aj  is 
uncertainty  block  and  A2  performance  block.  However,  due  to  the  difficulty  of  computing 
l*,  die  p-synthesis  is  relaxe  J  to  the  problem  of  minimizing  its  upper  bound: 

inf  sup  iirf{  d[DFj(G  KJD1] }  (6-3) 

where  D  is  the  scaling  matrix  mentioned  in  eq.(6-2).  Then,  eq.(6-3)  can  be  solved  by  the 
so  called  D-K  iteration  algorithm,  which  is  the  only  solution  available  to  the  problem  up  to 
now.  Basically,  D-K  iteration  algorithm  solves 

Rfcfi*) "  D(s)Fi(G(s)’  K(s»Dl(s)  (6-4) 

where  both  D(s)  and  D'Ks)  are  proper  stable  rational  functions.  This  optimization  is  a 
convex  searching  problem,  if  either  K(s)  or  D(s)  is  fixed.  Unfortunately,  it  is  not  a  convex 
problem  over  D  and  K. 
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D-K  iteration  algorithm 
Step  1:  D(s)  =  I ; 

Step  2:  Find  a  controller  K (s)  such  that  ND(s)F,[G(s),  KCs^D'^s^L  is  minimized; 

Step  3:  Find  constant  D’s  at  each  frequency  such  that  oCDFjCP,  K)D  *)  is  minimized; 

Step  4:  Curve  fitting  for  D(s),  where  D(s)  is  chosen  such  that  both  D(s)  and  D'^s)  are 

proper  stable  rational  functions; 

Step  5:  Go  to  step  2,  until  die  local  minimum  is  reached. 


6.3  Problem  Formulation 

In  this  section,  we  will  design  a  robust  controller  for  a  fighter  aircraft  at  an  altitude 
of  10,000  feet  and  a  Mach  Number  of  0.18.  The  angle  of  attack  is  approximately  70 
degree,  and  the  dynamic  pressure  is  33.6  lb/ft2.  The  system  is  described  by  the  following 
block  diagram. 


Fig.6.4 


The  P(s)  is  longitudinal  model  given  by  state  space  model  [9]: 

'Mr  z.  1  |xii  rz.1 

x  M,  M,  x  +  M,  v  (6-5a) 

y  *  \  (6-5b) 

where  xt  is  the  angle  of  attack,  x2  is  the  pitch  rate  and  v  is  the  pitch  vectoring  nozzle 
deflection.  It  is  assumed  that  Ma,  Mq  and  Mv  are  subjected  to  25%  variations  from  its 
nominal  values.  The  nominal  plant  is  given  by: 

[  0.0264  1  1  _  [-0.05201 

^  [-0.8810  -0.2079 J  BPV=  [-4.3434 J  Cpy=[J  °]  DPy  =  0-  (6'6> 


The  R  is  the  model  for  the  actuator  from  pilot  stick  to  the  vectoring  nozzle.  Its 
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transfer  function  is  given  by: 

■MV+W+W-  (6-7) 

From  above  model,  one  can  see  that  the  plant  is  stable,  but  has  a  large  overshoot 
and  oscillations,  long  settling  time  and  an  insufficient  stability  margin,  which  can  be 
verified  by  the  following  step  response  plot,  where  the  input  is  the  unit  step  function  acted 
on  pilot  stick  and  the  output  is  angle  of  attack  of  die  aircraft 


The  design  specifications  for  the  feedback  system  are:  (1)  The  system  has  a 
satisfactory  tracking  ability  with  a  damping  ratio  larger  than  0.5,  natural  frequency  large 
than  2.0  rad/sec,  small  overshoots  and  oscillations,  and  small  steady-state  tracking  error, 
(2)  Robustness  of  the  system,  including  the  stability  and  performance  in  the  face  of 
disturbance  and  the  plant  uncertainties. 

In  order  to  design  a  robust  controller,  one  needs  to  separate  plant  uncertainties  A1 
from  the  nominal  plant.  Consider  eq.(6-5)  and  eq.(6-6),  then  it  is  easy  to  see  that  the  plant 
can  be  expressed  by  Fig.6.6, 


where 


Fig.  6.6 


P(s)  = 


BP. 

Bp2l 

c , 

D  ,, 

D  . 

pi 

pii 

P12 

c  „ 

D  , 

D  „ 

P2 

P21 

p22j 

(6-8) 


Hoe  6j,  82  and  63  are  perturbations  to  Ma,  Mq  and  My  respectively.  AP,  Bp2  =  Bp*  and 


Cp2  =  Cy  were  given  in  eq.(6-6),  and 


n 

0 

0 

0 

0 

DPii 

Dpl2 

0 

0 

0 

0 

% 

- 1 

fS 

0. 

Q 

0 

0 

0 

1 

0 

0 

0 

0 

(6-9) 


With  plant  model  and  design  specifications  available,  a  controller  K(s)  can  be 
designed  for  the  following  LFT  block  diagram,  where  p(s)  consists  of  nominal  plant  P(s) 
and  actuator  R(s). 
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Fig.6.7 


Compared  with  Fig.6.2,  the  above  block  diagram  implies  that  the  exogenous  input 
w  consists  of  command  r  and  noise  n.  The  controlled  output  z  consists  of  weighted 
tracking  error  and  weighted  control  constraint  The  p  and  s  are  fictitious  input  and  output 
reflecting  system  uncertainties.  Wj(s),  W2(s)  and  W3(s)  are  the  weighting  matrices  chosen 
by  designers  such  that  die  design  specifications  can  be  met  Wx(s)  is  a  weighting  matrix  for 
tracking  error  and  is  chosen  as  a  low-pass  filter  to  emphasize  the  tracking  accuracy  at  low 
frequencies  (small  steady-state  error).  The  weighting  matrix  W2(s)  is  designed  for  control 
constraint  In  the  real  situation,  the  control  input  must  be  restricted  because  of  the  limited 
control  energy  and  actuator  saturation.  Usually  W2(s)  is  chosen  as  a  diagonal  constant 
matrix;  the  larger  W2(s)  is,  the  less  the  control  will  be  used.  W3(s)  is  used  to  normalize  the 
system  uncertainties.  In  this  particular  problem,  we  choose 

W2(s)  =  0.5  and  W3(s)  =  0.25  diag{0.881, 0.2079,  4.3434}. 


With  these  weighting  matrices,  the  p-syn  thesis  technique  can  be  used  to  design  a  robust 
controller  which  maintains  the  quality  of  system  performance  and  stability  in  the  face  of 
uncertainties  and  sensor  noise.  It  is  easy  to  derive  the  generalized  plant  G(s)  from  Fig  6.7: 
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G(s) 


(6-10) 


W3P11 

0 

0 

W3f.2 

-WfP2, 

W1 

-W1 

-W.P22 

0 

0 

0 

Wi 

1*2. 

I 

-I 

Then  by  using  the  notations  and  formulas  in  [27],  it  is  easy  to  see  that  the  generalized  plant 
G(s)  has  a  state  space  realization  as  follows: 


'A1 

B.C,2 

B,Dp22CR 

B.Dp2. 

-B1 

B. 

B.Dp22D* 

0 

Ap 

Bp2CR 

BP. 

0 

0 

BP2DR 

0 

0 

\ 

0 

0 

0 

bk 

0 

D,C  , 

3  pi 

D3Dp,2CR 

D.D  ,, 

3  pll 

0 

0 

D  D  D 

3  pl2  R 

0 

-D,C  , 

1  p2 
0 

-D1DP22^ 

0 

■D.D  „ 

1  p21 

0 

D1 

0 

D1 

0 

-D  D  D 

1  P22  R 

0 

-c , 

P2 

'Dp22CR 

'Dp21 

I 

-I 

'-DP22DR 

where  W.(s)  = 


r  A- 

B' 

1 

1 

c 

D. 

I 

a 

1,2,3,  are  state  space  realizations  of  weighting  matrices. 


With  the  state  space  realization  of  G(s)  available,  we  are  ready  to  use  p-synthesis  to 
design  a  robust  controller  K(s). 


6.4  Control  Law  Design 

Robust  controller  design  is  performed  on  the  system  shown  in  Fig.6.3,  with  the 
following  uncertainty  structure:  A  =  diag{Alt  A2},  where  At  =  diagiSj,  62,  53} 
corresponding  to  plant  uncertainties  and  A2  is  a  2x2  full  matrix  corresponding  to  the 
performance.  As  we  mentioned  earlier,  the  performance  is  evaluated  by  H~  norm  of  the 
transfer  function  from  the  command  and  noise  to  the  tracking  error  and  control  constraint 
The  objective  is  to  design  a  controller  K(s)  such  that  the  closed  loop  system  is  internally 
stable  and  die  H~  norm  of  Tw  remains  small  for  all  the  uncertainties  considered. 

With  the  procedure  illustrated  in  Section  6.2,  we  first  obtained  the  optimal  H*° 
controller  Kj(s)  with  the  optimal  H~  norm  being  7.64.  Since  K^s)  ignores  the  structure 
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information  of  A  and  treats  A  as  a  full  matrix,  it  gives  a  conservative  solution  to  die 
problem.  Fig.6.8  stows  the  p  plot  and  of  plot  of  the  closed  loop  system.  The  difference  of 
them  can  be  seen. 


As  we  expected,  p  plot  is  lower  than  o  plot,  or  expressed  alternatively  1  pa  £  1/p  at 
each  frequency.  This  implies  that  the  allowable  set  of  uncertainties  considered  (structured) 
is  larger  than  that  of  unstructured  uncertainties. 

Next  we  will  continue  the  D-K  iteration  design  procedure  described  in  section  6.2. 
After  two  iterations,  the  process  converges  to  a  controller  K^s)  which  minimizes  the  p  of 

the  closed  loop  system,  see  Fig.6.9  for  reference. 
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Frequency  (nd/s) 
Fig.6.9 


In  the  figure,  the  lower  curve  is  the  |X  plot  and  the  upper  one  the  o  plot  of  the  D- 
scaled  closed  loop  system.  Compare  with  Fig.6.8,  one  can  see  that  K^s)  gives  a  smaller  p 
(  =  5.23)  than  K^s)  did. 

In  die  following,  we  will  perform  some  analyses  on  the  original  closed  loop  system 
shown  in  Fig.  6.4  with  as  the  controller. 

Step  responses 

The  following  are  the  step  responses  of  angle  of  attack  and  the  corresponding  angle  of 
thrust  vector  nozzle  (control  input). 


73 


Ancle  of  Thrust  Vectoring  Nozzle  (Dec.)  of  Attack  (Dec  ) 


Fig.6.10a 


From  the  design  specifications,  it  is  required  that  the  settling  time  of  the  response 
4  ^ 

should  be  ^  — : —  =  4.6  second.  Observing  the  above  two  plots,  it  can  be  seen  that  the 

*)nC 

design  specifications  are  satisfied 


Fig.6.11b 


As  it  shows,  the  system  stability  and  performance  are  still  satisfactory  under  the 
perturbations. 


Controller  model  reduction 


High  order  controller  is  usually  obtained  due  to  the  need  of  "curve  fitting"  for  D(s) 
in  the  D-K  iteration  design  procedure.  In  this  example,  the  A  has  four  blocks  consisting  of 
three  lxl  and  one  2x2  block.  Therefore  D(s)  has  n  -  1  =  4  -  1  =  3  elements  to  fit  [4].  We 
choose  the  three  elements  with  third  order,  second  order  and  third  order  rational,  proper 
stable  minimal  phase  functions  in  the  last  iteration  respectively,  which  yields  an  eighth 
order  D(s).  With  the  D(s),  inverse  of  D(s)  and  the  original  fifth  order  plant,  the  final 
controller  is  of  twenty-first  order.  Therefore,  the  controller  should  be  reduced  to  a 
reasonable  order.  To  this  end  we  look  at  the  controllability  and  observability  of  K^s)  to  see 
whether  we  can  ret  ce  it  to  a  minimal  realization.  By  finding  a  balanced  realization,  can 
be  reduced  to  a  fifth  order  controller  K/s)  according  to  its  Hankel  singular  values  [28]: 


Ac  = 


-1.2708e+03 

0 

0 

0 


0 


-4.585 le+00 
-1.8591e+01 
1.291  le+01 
0 
0 


-5.7124e+01 
-2.201 3e+01 
-3.8605e+00 
0 
0 


2.1467e+02 
-2.338 le+00 
2.5750e+00 
-9.3984e+00 
0 


-4.6679e-01 

-2.0500e-02 

-1.5782e-02 

3.7513e-02 

-l.OOOle-03 


9.4612e-+01 

2.6817e+00 


3.6647e+00  , 
-7.4246e+00 
7.7460e-01 


Cc=  [  -9.4590e+01  -3.1892e+00  -2.1391e+O0  8.0623e+00  -7.748 le-01  ]; 
Dc  =  -6.6522e-03. 


With  the  reduced  order  controller  Kj(s),  we  have  the  following  p  and  a  plots  of  the 
closed  loop  system. 
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As  it  shows,  the  reduced  order  controller  K^s)  maintains  almost  the  same  p  for  the 
closed  loop  system  as  the  full  order  controller  does.  Hence  robust  stability  and  robust 
performance  are  achieved  with  K/s).  It  can  be  verified  by  the  following  step  responses 
with  and  without  25%  perturbations. 
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Ancle  of  Attack  (Dec.)  Ancle  of  Throat  Yectorinc  Nozzle  (Dec.) 
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Time  (Sec.) 
Fig.6.14b 


6.5  Conclusions 

A  fighter  aircraft  control  design  problem  was  considered  in  this  chapter.  The  design 
procedure  illustrates  how  to  apply  p-synthesis  to  a  real  engineering  problem.  We  addressed 
the  issues  of  problem  formulation,  D-K  iteration,  and  controller  reduction.  With  the 
controllers  obtained,  simulations  on  system  performance  and  stability  analyses  were  also 
performed. 
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CHAPTER  7 


CONTROLLER  REDUCTION  BY  STRUCTURED  TRUNCATION 


Simple  controllers  are  normally  preferred  over  complex  controllers  because  of  the 
less  computational  requirement  and  less  effort  in  implementation.  The  methods  for 
designing  low-order  controllers  can  broadly  be  classified  into  two  categories:  fixed  order 
controller  design  and  controller  reduction. 

The  philosophy  in  fixed  order  controller  design  is  to  seek  to  minimize  a 
performance  index  subject  to  the  constraint  that  the  controller  be  of  fixed  degree  [29,30]. 
For  controller  reduction,  one  can  simply  apply  model  reduction  techniques  to  either  the 
plant  model  or  the  controller.  There  are  now  at  least  three  rather  popular  state-space  based 
model  reduction  techniques,  namely,  truncation  of  the  internally  balanced  realization 
[31,3233],  Hankel-norm  approximation  [28,34],  ^-covariance  equivalent  realization  (q- 
COVER)  [63],  and  coprime  factorization  method  [35],  It  is  well  known  that  the  controller 
approximation  is  better  than  the  plant  approximation  [36].  However,  to  find  a  reduced- 
order  controller  which  approximates  the  original  controllers  is  not  our  direct  objective. 
What  we  seek  is  to  find  a  reduced-order  controller  such  that  the  reduced  closed-loop  system 
approximates  the  original  closed-loop  system. 

Jonckheere's  work  [37]  was  based  on  this  consideration,  in  which  two  Riccati 
equations  are  balanced  and  truncation  is  carried  out  with  respect  to  LQG  characteristic 
values.  But  the  relationship  between  the  LQG  characteristic  values  and  the  closed-loop 
system  properties  is  unclear.  So  Jonckheere's  approach  usually  cannot  provide  satisfactory 
reduced-order  controllers. 

In  this  chapter,  we  propose  a  new  controller  reduction  approach  which  is  based  on 
closed-loop  considerations  other  than  controller  approximations  or  plant  model 
approximations.  In  Section  7.1,  a  method  of  structured  truncation  based  on  the  closed-loop 
properties  is  developed,  and  some  interesting  properties  of  the  controllability  and 
observability  gramians  of  the  H2  optimal  closed-loop  system  are  also  presented.  Some 
illustrative  examples  are  given  in  Section  7.2. 
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7.1  Structured  Truncation  Approach 

Given  a  closed-loop  system,  the  structured  truncation  approach  is  to  reduce  the 
order  of  controller  by  minimizing  a  cost  function  based  on  the  closed-loop  system.  The  cost 
function  will  be  elaborated  later.  Although  this  approach  can  be  applied  to  any  desired 
closed-loop  system,  we  will  concentrate  on  the  optimal  H2  closed-loop  system  for  the 
simplicity  of  presentation. 

It  is  well  known  that  the  input/output  relationship  of  a  system  is  well  represented  by 
its  controllability  and  observability  gramians.  The  balanced  truncation  [31]  is  based  on 
these  gramians.  In  the  proposed  structured  truncation  approach  we  also  use  them.  For 
convenience,  some  basic  knowledge  about  controllability  and  observability  gramians  is 
reviewed  in  the  following. 

Controllability  and  Observability  Gramians 


Consider  a  linear  time-invariant  system 

x(t)  =  Ax(t)  +  Bu(t) 
y(t)  =  Cx(t)  +  Du(t) 

i.e.,  the  transfer  function  is 


G(s)  := 


A 

B 

C 

D 

=  D  +  C(sI-A)_1B 


The  eigenvalues  of  A  are  assumed  to  be  strictly  in  the  left  half-plane  and  the  controllability 
gramian  P  and  the  observability  gramian  Q  are  defined  as 

P  :=  JJ°exp(At)BB'exp(A't)dt 
Q  :=  J”exp(A’t)C'Cexp(At)dt 

where  A'  is  the  transpose  of  A.  P  and  Q  satisfy  the  following  Lyapunov  equations 

AP  +  PA'  +  BB’  =  0 
AQ  +  QA  +  CC  =  0 

and  have  the  following  properties. 

Property  7.1 
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If  all  the  eigenvalues  of  A  are  strictly  in  the  left-half  plane,  then 

(a)  P  >  0  if  and  only  if  (A,  B)  is  completely  controllable. 

(b)  Q  >  0  if  and  only  if  (A,  Q  is  completely  observable. 


If  the  state-space  coordinates  of  the  system  are  changed  to  z  =  Tlx  for  some  non¬ 
singular  T  then 

z  =  T4ATz  +  t\  y  =  CTz  +  Du 


Furthermore,  the  controllability  and  observability  gramians  become  T 1PT'1  and  TQT 
respectively,  and  the  product  PQ  will  be  transformed  to  T1PQT.  Therefore  the  eigenvalues 
of  PQ  are  invariant  under  state-space  transformations,  and  are  input/output  invariant  A 
useful  state-space  realization  in  this  respect  is  the  balanced  realization  where  P  =  Q  = 
diagCOj,  a2, ....  on),  therefore  PQ  =  diag^2 . o*). 


Structured  Truncation 


Consider  the  block  diagram  of  in  Fig.  7.1. 


Fig.  7. 1  A  generalized  plant  G  with  a  controller  K. 


G  is  a  generalized  nominal  plant,  w  contains  all  external  inputs,  including  disturbances, 
sensor  noises,  and  commands;  the  output  z  is  the  controlled  output;  y  is  the  measured 
output;  and  u  is  the  control  input  A  state-space  realization  of  the  generalized  plant  G  is 
given  as 


G(s)  := 


A 

B. 

B2  1 

c 

D 

1 

11 

12 

,C2 

D21 

D22. 

(7-1) 


where  A  «  Rnxn  ,  Bj  «  Rnxn»i,  B2  *  Rnx®2,  cy  e  RPixn  and  C2  €  RP2xn  The 
controller  K  is  assumed  to  be 
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K(s)  = 


fAJ 

»cl 

°cj 

(7-2) 


where  Ac«R“cx“c,  Bc«R“cxP2,  Cc«R“2*ncf  and  Dc«Rm2*p2.  denotes  the  closed- 
loop  transfer  function  from  w  to  z. 

The  objective  is  to  reduce  the  order  of  the  controller  with  which  the  closed-loop 
system  approximates  the  original  closed-loop  system  Hence  we  consider  the  controllability 
and  observability  gramians  of  the  closed-loop  transfer  function  Tzw  The  balancing 
technique  can  not  be  applied  on  the  gramians  of  the  closed-loop  system  because  the 
balancing  transformation  takes  place  in  the  whole  state  space  of  the  nominal  plant  and  the 
controller.  With  the  states  of  the  closed-loop  system  arranged  as  [xp\  x^'  where  xp  is  the 
set  of  the  n  states  of  the  plant  and  \  is  the  set  of  the  states  from  the  controller,  to  avoid 
mixing  up  the  states  from  the  plant  and  the  controller,  only  the  following  structured 
transformation 


T  = 


(7-3) 


is  allowed.  Without  loss  of  generality  we  use  the  transformation  shown  below 


J*  0 

0  T 


(7-4) 


Let  P  and  Q  be  the  closed-loop  controllability  and  observability  gramians  respectively  and 
define  J  =  PQ.  If  the  closed-loop  system  can  be  balanced,  i.e.  J  can  be  diagonalized,  then 
the  states  corresponding  to  the  small  diagonal  elements  can  be  truncated  from  the  system. 
Due  to  the  structure  restriction  on  transformation,  it  is  impossible  to  totally  balance  the 
closed-loop  system 

By  applying  the  transformation  T  of  (7-4)  to  the  closed-loop  system,  J  becomes 
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0  ' 

0 

0 

T1 

J 

0 

Tn 

- 

c. 

*n  01 


0  T. 


’ll 


.-1 


J11  J12 
J21  J22 


J  T 
12  X 


h  0 
0  V 


T^J„  T'LJ  T 

nc  21  nc  22  nc 


Let  T*1  = 

“c 


Sl* 

i' 


and  Tn  =  [^i  l],  where  sand  I  «RDC,  so 


r  j 

J  T  1 

Ii. 

J  T 

12  1 

V 

11 

12  «c 

SiJ2. 

S  7  T 

1  22  1 

S.'J  1 
22 

t*j 

T_1J  T 

»C  21 

nc  22  nc 

_ 1 

The  problem  is  proposed  to  find  such  that 


f(T„)  = 


0 


sll 


21 


*J22T1 


S  ’j  1 

22 

^22* 


is  minimized. 


(7-5) 


Here  Up  denotes  the  Frobenius-norm.  Then  we  truncate  the  states  of  the  controller 
corresponding  to  the  small  cost  defined  in  (7-5). 


For  the  H2  optimization  problem  with  the  generalized  plant 


’  A 

»2 

G(s)  := 

C! 

0 

D12 

C2 

D2, 

0 
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and  the  optimal  controller 


(s):= 


where  F2  :=  -B2X2,  L2  :=  -Y2C2'  and  Ac  :=  A  +  "  +  L2C2,  and  X2  and  Y2  are  the 

solutions  of  the  following  two  Riccati  equations 


A’X2  +  V-W2'X2  +  Cl’Cl=0 

AVY2A’  -Y2C2C2Y2  +  B1B1,  =  0- 

The  closed-loop  system  with  the  optimal  H2  controller  is  of  the  form 


A 

BF 

B,  1 

2  2 

1 

-L  C 

A  +  B  F  +  L  C 

-D  L 

2  2 

2  2  2  2 

21  2 

c> 

F  D 

2  12 

0 

Let  P  and  Q  be  the  controllability  and  observability  gramians  of  9(s)  respectively,  then  we 
have  die  following  proposition. 

Proposition  7.1:  P  and  Q  are  of  form 


P  = 


P„+X„ 
2  2 


Q  = 


VY2  -Q, 
-Q,  Q, 


where  P9  and  Q_  satisfy  the  following  two  Lyapunov  equations 


(A  +  B2F2)P2  +  P2(A  +  b2f2)'  +  l2l2'  =  0 
(A  +  L2C/Q2  +  Q2<A  +  L2C2)  +  F2,F2  =  0 
This  result  can  be  easily  verified. 

So  the  structures  guarantee 

1  l1"  II2 
~Ui  °. 

The  problem  mentioned  above  turns  out  to  find  Tnc  such  that 


(7-6) 

(7-7) 
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0  V 

s'J2.  0 


<w  fCTii^  =s’J2i’J2ii  +  i  ’J12’J12i  is  minimized. 

Theorem  7.1 

There  exist  two  vectors  &  and  1  such  that  s'l  =  1, 

f(& ,  1 )  =  i,J2i'J2i£  +  *  ‘J12J121 

is  minimized. 

Proof:  Let 


f(i,i)  =  4,J21'J21i  +i'J12,J12t  *2X(&’1  _1) 


then  set 


-° 

i'l  =L 

It  is  easy  to  see  that  X  is  the  eigenvalue  of  the  matrix 


°  J12J12 
J2l'J21  ° 


(7-8) 

(7-9) 


(7-10) 


and  isj^jthe  eigenvector  corresponding  to  the  smallest  eigenvalue  of  (10. 10).  Q.E.D. 

The  matrix  (7-10)  is  a  special  case  of  Hamiltonian  and  has  the  following  property. 
Property  7.2 

If  both  J\2  and  J21  are  full  of  rank  then  the  eigenvalues  of  (7-10)  are  real. 

Proof: 

From  (7-8)  and  (7-9) 


J2l'J21*  =  ** 


(7-11) 
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J12J12*  =*S 


(7-12) 


Multiply  (7-11)  by  &'  and  (7-12)  by  and  plus  them  together.  Here  $'  means  the  conjugate 
transpose  of  s. 

SJ21J2I&  + 1  J12J121  =  ^  + 1  (7‘13) 

From  (7-13)  we  see  that  Im(s' J2 1  ' ^2 1 1 2  7 1 21 )  =  0.  and  Im($’i_+  I's)  *  0.  So  Im(X)  = 
0.  Q.E.D. 


7.2  Examples 

Two  examples  are  shown  to  illustrate  the  proposed  structured  truncation  approach 
for  controller  truncation. 


Example  7.1: 

Given  a  plant: 


’  A 

B.  ®2] 

■  -1  0 
0  2 

O  M 

0  0 

►-*  0 
_ 1 

G(s)  := 

c. 

C 

0  DU 
D  0 

= 

1  1 
_  0_0 

0  0!  0 
,o_o|_i_ 

2 

21  J 

"1"7 

0  1!  0. 

Full-order  optimal  H2  controller  is: 


Bc] 

c 

0 

c 

'  -0.8619  0.1381 

-0.1381 

-4.8613  -6.7884 

4.5523 

-0.3090  -4.2361 

0 

By  the  structured  truncation,  the  controller  can  be  reduced  to  the  following  lst-order 
controller 


Ac 

Be' 

'  -6.4397 

1-3.4225 

[Cc 

0 

.  5.5247 

1  0  J 

In  Table.l,  Hj-norms  of  reduced-order  closed-loop  system  and  Loo-error  bound  between 
the  closed-loop  system  with  the  full-order  controller  and  those  with  reduced-order 
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controllers  obtained  by  our  new  structured  truncation,  balanced  truncation,  coprime 
factorization,  and  optimal  projection  are  listed. 


|Tzwl22 

ESCSB 

Full-order  Optimal  Controller 

91.8196689 

0 

lst-order  controllers 

structured  Truncation 

91.857999 

0.233720  | 

Balanced  truncation 

91.871408 

0.392831 

coprime  factorization 

91.870133 

0.361891 

91.856689 

0.269276 

Table  1. 


Example  7.2: 

Given  plant 


G(s)  := 


A 

B  B 1 
1  2 

C 

0  D,„ 

1 

12 

c„ 

D  0 

2 

21  J 

‘  1 

0 

0 

1 

0 

0 

°! 

r 

0 

-2 

0 

0 

1 

0 

°! 

1 

0 

0 

-3_ 

0 

£_ 

01 

l_ 

1 

0 

0 

0 

0 

0 

oT 

0 

0 

1 

0 

0 

0 

0 

ol 

0 

0 

0 

1 

0 

0 

0 

oj 

0 

0 

_0 

_0 

0 

0 

0 

0! 

1 

"Y 

’1 

T 

’0 

”0 

”0" 

"if 

~o 

The  full-order  optimal  H2  controller  is: 


■  -4.1899  -2.6639  -2.6667 

2.5950' 

K 

-2.6639  -2.1379  -0.1407 

0.0689 

0 

-2.6667  -0.1407  -3.1435 

0.0717 

_  -2.5950  -0.0689  -0.0717 

0 

The  controller  is  unstable. 

The  reduced  second  and  first-order  controllers  by  the  structured  truncation  are 


n 

Su 

■Q! 

lCc 

l°J 

and 


'  -1.3875  -3.6269 

0.0263' 

= 

-2.7605  -5.2019 

-2.5453 

-0.6207  2.6198 

0 
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Cc|OJ  [-2.0529 1  0 

The  comparison  of  the  proposed  structured  truncation  with  the  others  is  shown  in  Table  2. 


ITzw  l22 


Full-oider  Optimal  Controller 

22.490551 

0 

2nd-order  controllers 

Structured  Truncation 

22.491118 

0.019268 

Balanced  truncation 

(The  controller  is  unstable) 

Coprime  factorization 

22.491945 

0.051115 

Optimal  Proj  ection 

22.491114 

0.021112 

lst-order  controllers 

Structured  Truncation 

33.797984 

4.913441 

Balanced  truncation 

(The  controller  is  unstable) 

Coprime  factorization 

42.666848 

16.386789 

Optimal  Projection 

30.333357 

4.4430115 

Table  2. 


The  optimal  projection  is  the  fixed-order  controller  design  to  minimize  the  H2-norm  of  the 
closed-loop  system.  For  the  two  examples,  the  H2-norm  of  the  closed-loop  system  with  the 

reduced-order  controller  by  the  new  structured  truncation  is  very  close  to  the  result  of 
optimal  projection  design.  The  structured  truncation  also  results  in  small  Loo-enor  bounds. 
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CHAPTER  8 


A  PARAMETRIZATION  APPROACH 
TO  REDUCED-ORDER  H“  CONTROLLER  DESIGN 


8.1  Introduction 

By  the  state-space  approach  [7]  to  H“  optimization,  an  (sub)optimal  controller  can 
be  easily  obtained  and  the  order  of  the  controller  is  not  higher  than  that  of  the  generalized 
plant  However,  in  many  engineering  problems,  the  order  of  generalized  plant  can  be  very 
high.  This  is  due  to  the  fact  that  the  generalized  plant  consists  of  the  original  plant  as  well 
as  all  die  weighting  matrices  which  are  chosen  to  meet  certain  design  specifications.  Hence, 
the  order  of  the  H~  controller  obtained  by  the  standard  state-space  approach  is  usually  too 
high  to  be  implemented  in  practice,  and  therefore  a  systematical  methodology  of  designing 
a  reduced-order  H*°  robust  controller  is  desired. 

This  issue  has  attracted  considerable  attention  in  the  past  few  years,  and  several 
approaches  have  been  proposed  to  address  this  problem.  Bernstein  et  al.  [30,38]  proposed 
a  fixed  order  controller  design  method  based  on  the  optimal  projection  theory.  It  was 
shown  that  an  optimal  reduced-order  iP/H-  controller  can  be  obtained  by  solving  several 
coupled  Riccati  equations.  However,  the  computation  involved  is  very  complicated.  A 
direct  design  method  based  on  the  well-known  Bounded  Real  Lemma  [39]  was  proposed 
by  Hsu  et  al.  [40],  where  a  set  of  sufficient  conditions  and  design  algorithms  were  derived 
for  a  reduced-order  H~  controller  design.  Chang  et  al.  [41,42]  considered  observer-based 
controller  parametrization,  and  pointed  out  that  by  selecting  a  suitable  parameter  matrix,  one 
can  make  the  realization  of  the  stabilizing  controller  non-minimal,  whose  order  then  can  be 
reduced.  Choi  et  al.  in  [43]  proposed  a  method  of  constructing  such  a  suitable  parameter 
matrix  in  order  to  get  a  reduced-order  stabilizing  controller. 

These  results,  especially  the  work  by  Choi  et  al.  [43],  motivated  our  work  in 
designing  reduced-order  H“  controllers.  In  this  chapter  we  will  present  a  direct  design 
algorithm  for  (n-p2)th  order  H“  controller,  where  n  is  the  order  of  the  generalized  controller 
and  p2  is  the  number  of  independent  measured  outputs.  The  major  idea  is  to  select  a 
parameter  matrix  such  that  the  realization  of  the  controller  obtained  from  DGKF  formulas 
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[7]  is  not  minimal.  By  deleting  the  unobservable  states  of  the  controller,  (n-p^th  order  H- 
controller  can  be  obtained.  The  content  of  this  chapter  is  arranged  as  follows.  Section  8.2  is 
a  preliminary  section  which  provides  the  standard  state  space  approach  to  H“  optimization 
problem,  including  the  formula  of  all  stabilizing  (sub)optimal  controller.  We  present  our 
main  result  in  section  8.3,  where  a  set  of  formulas  are  given  for  constructing  a  reduced 
order  H*“  controller.  Section  8.4  shows  an  illustrative  example  and  Section  8.5  is  a 
conclusion. 

8.2  Preliminaries 


Hie  following  block  diagram  is  the  standard  H~  optimization  configuration. 

w 

u 


G(s) 


K(s) 


Fig.  8. 1  The  standard  block  diagram  for  H“  optimization 

The  G(s)  is  the  generalized  plant,  including  the  original  plant  and  possible  weighting 
matrices.  The  signals  z,  y,  w,  and  u  are  the  controlled  output,  the  measured  output,  the 
exogenous  input  and  the  control  input  respectively.  The  standard  H“  optimization  problem 

is:  To  design  a  proper  controller  K(s)  such  that  the  closed-loop  system  is  internally  stable 
and  HTw(s)ll_  is  minimized,  where  Tw(s)  denotes  the  transfer  function  of  the  closed  loop 

system  from  w  to  z. 


In  the  DGKF  approach  [7],  the  realization  of  the  generalized  plant  G(s)  is  assumed 


to  be 


G12(s)- 

’  A 

B1 

h 

II 

3 

o 

= 

c, 

o 

D,, 

,VS> 

<y». 

l 

_  C2 

°21 

12 

0 

(8-1) 


where  G(s)  is  partitioned  such  that  Gn(s)  is  p^mj,  G12(s)  is  p1xm2,  021(s)  is  p2xm1, 
and  G^Cs)  is  p2xm2.  It  is  also  assumed  that: 

(i)  Both  G12(s)  and  G21(s)  do  not  have  any  transmission  zeros  on  the  jeo-axis. 

(ii)  (A,  is  stabilizable  and  (C^,  A)  is  detectable. 
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Crii)  D’^ICj  D12]  «[0  I] 


V 

d'  = 

f°l 

,D». 

21 

LiJ 

The  two  Riccati  equations  involved  are: 

ax„ + x„a + b2b^  +  C’Cj  =  0 

and 

AY„  +  Y^A'  +  YJy^Cj  -  +  BjBj  =  0. 


(8-2) 

(8-3) 


The  following  famous  theorem  characterizes  all  (sub)optimal  stabilizing  controllers  such 
that  IIT^IL  <  Y- 

Theorem  8.1  [7] 

For  a  given  y,  there  exists  a  stabilizing  controller  such  that  IITzwIL<y  if  and  only  if  the 
following  three  conditions  hold: 

(i)  There  exists  a  positive  semidefinite  stabilizing  sc  ition  Xj(y)  to  (8-2). 

0i)  There  exists  a  positive  semidefinite  stabilizing  solution  Y„(y)  to  (8-3). 

(iii)  p\XJJ)Ym(1)]  <  V2,  where  p  denotes  die  spectral  radius  of  a  matrix. 

When  these  conditions  hold,  all  stabilizing  (sub)optimal  controllers  K(s)  can  be 
parametrized  by  a  stable  proper  parameter  matrix  Q(s)  whose  H~  norm  is  less  than  y,  as 
shown  in  the  following  figure. 


Fig.  8.2  The  parametrization  of  all  stabilizing  (sub)optimal  controllers 


A 

K. 

-Z.L.  2J% 

Here 

MJs)  = 

F_ 

.<1 

0  I 

I  0 

with  A„=  A  +  Y2BjBjXw  +  B2F„  +  Z JJC2  , 

f«=-b2x^ 


(8-4) 
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(8-5) 

(8-6) 


and 


l.=-y_c2 

z»=(i-y2y.xj1 


Q(s)  = 


[\ 

®,i 

K 

(8-7) 

(8-8) 

(8-9) 


8.3  Reduced-Order  H**  Controller  Design 


In  this  section,  we  will  show  the  existence  of  an  (n-p^th  stabilizing  H~  sub-optimal 
controller,  where  n  is  die  order  of  the  generalized  plant  and  p2  is  die  number  of  independent 
measured  outputs.  By  using  the  DGKFs  formulas  in  Theorem  8.1,  the  state-space 
equations  for  sub-optimal  controller  K(s)  can  be  expressed  as 


x  *  A„x  -  Z^Ljr  +  ZJB2q 

(8-10) 

u  =  F^x  +  q 

(8-11) 

p  =  -C2x  +  y 

(8-12) 

x  »  V  +  BqP 
q  =  Cqx  +  Dqp 


(8-13) 

(8-14) 


where  the  dimension  of  Aq  is  (n-p2).  Simplifying  the  above  equations,  we  have  the 
following  state-space  representation  for  the  controller, 


I 

N 

a+  72b1b;x^+  b2f„+  zjjz2  -zjj2d,c2 

w, 

ZJBjD^ZJU. 

h 

F«-D,c2 

c, 

D, 

(8-15) 


Note  the  fact  that  the  realization  of  K(s)  may  not  be  minimal,  which  depends  on  the  {Aq, 
Bq,  Cq,  Dq}  chosen.  This  suggests  that  if  one  can  choose  a  set  of  {Aq,  Bq,  Cq,  Dq}  such 
that  the  realization  of  K(s)  in  eq.(8-15)  is  non  minimal,  then  by  deleting  the  unobservable 
or/and  uncontrollable  states  of  K(s),  a  reduced-order  controller  can  be  obtained,  and  the 
system  performance  and  stability  can  be  maintained  at  same  time.  To  this  end,  we  apply  a 
similarity  transformation 


(8-16) 
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to  K(s),  where  X  will  be  determined  later.  With  die  transformation,  eq.(8-15)  becomes 


K(s) 


bm1 

CM 

DmJ 

where 


Am  = 


^11 

^M21 


AM12 

AM22 


with 

Amu  «  A  +  Y‘2B1B;x.  +  B2F.  +  Z.L.C2-Z„BJDt|C2  -  ZJBjC^X 

AM12  *  Z~B2Cq 


(8-17) 


(8-18) 

(8-19) 


am21  =  X(A  +  Y^BjB'X^  +  B2F»  +  2LUC2  -  ZJB2DqC2) 


-  BqC2  -XZo»B2CqX  -  AqX  (8-20) 

AM22*XZJ82Cq+A|1  (8-21) 

and 

[Vl  \  Z-B2D<l‘Z-L-  1 

B“=  BM2|  =  l.X(Z-B2Dq-Z-L-,+BkJ  '  <8'22) 

^]-[F--D,C2-C,X  C„]  (8-23) 

Dm  =  D,  (8-24) 

Note  that  our  objective  is  to  choose  an  X  in  die  similarity  transformation  such  that 

CM1-0  (8-25) 

aM21  =  0*  (8-26) 


If  this  is  achieved,  the  unobservable  pair  (C^i,  AM11)  in  eq.(8-17)  can  be  deleted  and  the 
controller  in  eq.(8-17)  is  reduced  to 


K(s)  = 


AM22 

c 

L  M2 

dmJ 

(8-27) 


whose  order  is  r  =  dim(Aq)  =  n-p2-  In  fact,  the  equations  (8-25)  and  (8-26)  arc  identical  to 
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(8-28) 

(8-29) 


[Cq  Dq]  C2  =F. 

XA  -AqX-BqC2  =  0 

where 

A  =  A  +  y2B1B'1X_  +  B2F_  +  Z_L.C2  -  Z_B2F_  (8-30) 

The  discussion  of  the  existence  of  the  X  to  equations  (8-28)  and  (8-29)  will  be  given  in  the 
following,  which  plays  a  key  role  in  our  main  result 

Lemma  8.1  A  is  a  stable  matrix. 


Proof:  Define 

P  :=  r2LYjJ„  -  B2F«  +  Z^B2F-, 
then  it  is  easy  to  verify  that 

p = ( -r2 z_y«x«  - 1 + zj)B2¥^ = [Z4 1  -  r2Y jc  )  -  i]b2f«, = o. 

Note  eq.(8-30)  can  be  expressed  as 

A  =  A  +  rV.X,  +  b2F.  ♦  Z.L.C2  -  Z„B2F. 

~  Atmp  +  Z_Y_(-r2FlF--C2C2)-P 
=  Atmp  +z»Y«(-r2Flf--c2c2), 

where 

Atap  =  A  +  Y2B1B,1X-. 

It  is  straightforward  to  show  that  [7] 

T'1JtmpT  =  JL, 

where 

[a'  /c;cfc;c2l 


-B  B’ 

li 

-a  j 

A’unp 

Y2f1f„-C2C 

-B  B’ 

li 

"Atmp 

(8-31) 

(8-32) 

(8-33) 

(8-34) 

(8-35) 

(8-36) 

(8-37) 
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and 


I  -y  2X. 
0  I 


(8-38) 


Hence,  it  is  easy  to  see  that  the  stabilizing  solution  to  eq.(8-37)  can  be  expressed  in  terms 
of  Y-  which  is  the  stabilizing  solution  to  eq.(8-36)  (or  eq.8-3)),  i.e.. 


Ytmp  =  Z»Y. 


(8-39) 


is  the  stabilizing  solution  to  Riccati  equation  eq.(8-37)  [7].  Substituting  eq.(8-39)  into 
eq.(8-33),  we  have 

a = + Ytmp(  If.  -  c2c2), 

which  implies  A  is  stable.  Thus  we  proved  the  lemma. 


Without  loss  of  generality,  we  assume  that 

C2  =  [Ip2  0]. 

Then  we  partition  A  and  X  compatibly  such  that 


A  = 


A11  Ai2 

Ai 


X  =  [X,  X2], 

With  the  partition,  eq.(8-29)  can  be  expressed  as 

A11  ^2 
■^21  ^2 


[X,  x2] 

or  equivalently 


[AqX1  A^  -[Bq  0j=0. 


XlAll+X2A21-AA  =  Bq 
XlA12  +  X2A22-AqX2  =  °. 

Lemma  %2  The  pair  ( A12,  A^ )  is  detectable. 


(8-40) 


(8-41) 

(8-42) 


(8-43) 

(8-44) 

(8-45) 


Proof:  Since  A  is  stable,  then  for  Re(s)  >  0 
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This  implies  that 

=  n  -  p2  for  Re(s)  >  0.  (8-47) 

m 

Thus  the  lemma  is  proved 

Lemma  8.3  There  exists  a  solution  X  =  [Xj  XJ  to  eq.(8-4S). 

Proof:  Since  (A12,  A22)  is  detectable,  then  there  exists  a  Y  and  invertible  W  such  that 

A22  +  YA12  =  W-1AqW  =  A  (8-48) 

where  A  is  a  stable  matrix  whose  eigenvalues  consist  of  the  unobservable  stable 
eigenvalues  of  A22  and  those  we  are  chosen  arbitrarily  by  selecting  a  Y.  Eq(8-48)  actually 


implies 

WA22  +  WYA12  =  AqW 

(8-49) 

i.e., 

X2  =  W 

(8-50) 

Xj=WY 

(8-51) 

Thus  the  lemma  is  proven. 

Remark  1.  Note  that  A  is  not  unique. 

With  Lemma  8.2  and  8.3,  we  are  ready  to  present  our  main  result  in  the  following  theorem. 


Theorem  8.2  For  a  given  y,  define  the  parameter  matrix 


\  Bq 
<*>  =  C  D 

q  q 

(8-52) 

as 

Aq  =  A22  +  YA12 

(8-53) 

VYAU  +  W 

(8-54) 
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C,-F2 

(8-55) 

Dq*FrF2Y 

(8-56) 

and 

F_-[F,  F2], 

(8-57) 

where  An,  A12,  A2i  and  Aa  are  the  partition  of  A  If  one  can  choose  a  matrix  Y  such  that 
Aq  is  stable  and  IIQ(s)IL  <  Y,  then  there  exists  an  (n-p^th  order  stabilizing  controller  with 
the  following  form 


K(s)  = 


Aq  +  XZ^B2Cq 

Bq+X(Z^B2Dq-Z„LJ- 

[  C« 

D,  J 

where 

X  =  [  Y  I]. 

Moreover,  the  H*°  norm  of  the  closed-loop  system  is  less  than  y. 


(8-58) 


(8-59) 


Proof:  Choose  W  =  I  in  the  previous  lemma,  then  Aq  =  A^  +  YA12.  Since  (  A12,  A^  )  is 
detectable,  there  always  exists  a  Y  such  that  A^  +  YA12  is  stable.  Eq.(8-54)  is  a  direct 
result  from  eq.(8-44).  From  eq.(8-29),  we  have 


1°,  D<,l=F-[*  o]  =  [F>  F’> 


0 

I 


I 

-Y 


=  [F2  FrF2Y). 


Thus  we  completed  the  proof. 


RemaHt  2.  Obviously,  such  a  Y  is  not  unique  and  it  is  desirable  to  choose  a  stabilizing  Y 
such  that  IIQ(s)IL  is  small.  However,  in  case  that  one  can  not  find  a  Y  such  that  IIQ(s)IL  < 

y,  a  larger  upper  bound  on  the  norm  is  required. 


8.4  Illustrative  Example 

The  generalized  plant  is  given  as 


G(s)  = 


-1 

0 

0 

1 

0 

o  1 
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By  using  DGKF  formulas  and  the  algorithm  we  developed  in  Chapter  S,  it  is  easy  to 
compute  the  optimal  H~  norm  of  the  close-loop  system,  which  is  5.31S24.  A  sub-optimal 
upper  bound  of  the  system  is  chosen  as  Ysub  =  6.  By  choosing  Y  =  [-1  -1],  we  have  stable 
(A22+YA12)  and  IIQ(s)lL=  5.54191  which  is  less  than  ytub.  From  the  formulas  given  by 
eq.(8-52)  to  eq.(8-58)  in  the  previous  section,  a  first  order  H“  controller  can  be  constructed 
as  follows. 


K(s)  = 


-4.004 

-2.999  -6.155" 

-0.259 

■0211  -5.535 

(8-60) 


With  this  controller,  the  eigenvalues  and  the  H“  norm  of  the  closed  loop  system  are 


-4.0,  -5.691,  -2.1549,  -0.9704 


and  5.541640  respectively. 

8.5  Conclusions 

In  this  chapter,  a  parametrization  approach  is  presented  for  the  reduced-order  H°° 
controller  design.  We  showed  that  there  always  exists  an  (n-p2)th  order  stabilizing 
controller  for  an  nth  order  generalized  plant  with  P2  independent  measured  outputs.  By 
using  the  formulas  provided,  a  reduced  order  H*“  controller  can  be  easily  obtained. 
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CHAPTER  9 


CONTROLLER  REDUCTION  VIA 
OBSERVER-BASED  CONTROLLER  PARAMETRIZATION 

Of  recent,  interest  in  the  conception  and  development  of  sophisticated  aircraft  and 
spacecrafts  has  increased.  Examples  of  the  spin-off  of  this  interest  include  such  mechanical 
systems  as  hypermaneuverable  aircrafts  and  space  stations.  As  these  systems  get  more 
sophisticated,  they  become  more  and  more  complex.  Consequently,  conventional  modeling 
techniques  and  control  design  strategies  become  inadequate.  This  has,  by  necessity,  led  to 
increased  research  activity  in  such  areas  as  model  reduction  [28,31,44,45,46],  reduced 
order  controller  design  [30,36,37,47,48,49],  decentralized  control  [50,51,52],  and 
control/structure  interactions  [53-56]. 

This  chapter  of  the  report  deals  with  reduced  compensator  design,  which  has  been 
tackled  from  an  observer-based  compensator  point  of  view.  Briefly,  this  chapter  develops 
two  properties  related  to  observer  based  controller  parametrization  and  pole  placement  It 
shows  that  the  poles  of  the  closed-loop  system  with  the  observer-based  controller  (Q(s» 
parametrization  are  the  regulator  poles,  the  observer  poles,  together  with  the  poles  of  the 
added  stable  parameter  matrix  K(s).  If  the  controller  Q(s)  is  realized  by  a  minimal 
realization,  the  closed  loop  poles  will  include  all  the  poles  of  the  added  stable  parameter 
matrix  K(s)  and  a  subset  of  the  regulator  and  the  observer  poles.  We  parametrize  all  such 
K(s)  which  render  Q(s)  non-minimal,  thereby  permitting  a  minimal  realization  of  Q(s)  to 
serve  as  a  reduced  order  controller.  With  such  a  parametrization  available,  one  could  then 
choose  a  K(s)  in  order  to  best  approximate  .arge  order  controllers,  such  as  the  H“ 
compensators,  by  a  lower  order  controller. 

9.1  Observer  Based  Compensation 

One  of  the  most  fundamental  requirements  in  control  systems  design  is  to  make  die 
closed-loop  system  internally  stable.  In  addition  to  closed-loop  stability,  usually  the  closed- 
loop  system  is  required  to  meet  some  other  desired  performance  criteria.  Stabilizing 
controller  parametrization  is  important  because  of  the  following  reasons:  (1)  It  provides  the 
full  set  of  the  controllers  which  stabilize  the  closed-loop  system.  (2)  The  full  set  of 


100 


stabilizing  controllers  is  characterized  in  terms  of  a  stable  parameter  matrix  and  the  closed- 
loop  system  is  internally  stable  if  and  only  if  the  parameter  matrix  is  stable.  (3)  The  closed- 
loop  transfer  function  matrix  related  to  the  performance  can  be  written  as  a  simple  affine 
function  of  the  parameter  matrix  and  then  the  control  system  design  problem  becomes  that 
of  finding  a  stable  parameter  matrix  such  that  the  closed-loop  transfer  function  matrix  meets 
the  desired  performance  criteria. 

The  first  characterization  of  the  set  of  all  stabilizing  controllers  in  terms  of  a  stable 
parameter  matrix  was  introduced  by  Youla  et.  al.  [57]  in  1976.  Youla’s  controller 
characterization  was  developed  based  on  the  fractional  factorizations  over  the  ring  of 
polynomial  matrices.  The  only  drawback  of  Youla's  characterization  is  that  the  stabilizing 
controller  may  not  be  proper.  This  drawback  was  removed  later  by  Desoer  et.  al.  [58]  in 
1980. 


Desoer  et.  al.  [58]  generalized  Youla  et.  al.'s  result  based  on  the  fractional 
factorizations  over  a  general  ring.  The  ’ing  can  be  chosen  as  the  set  of  proper  stable  rational 
matrices  if  the  given  plant  is  a  linear  time-invariant  system  which  is  represented  by  a 
rational  matrix.  Based  on  the  fractional  factorizations  over  the  ring  of  proper  stable  rational 
matrices  the  set  of  all  proper  stabilizing  controllers  can  be  characterized  in  terms  of  a  proper 
stable  parameter  matrix.  The  closed-loop  system  is  internally  stable  if  and  only  if  the 
parameter  matrix  is  proper  stable  and  the  stabilizing  controllers  are  proper  if  a  simple 
inequality  is  satisfied. 

To  use  Desoer  et.  al.'s  version  of  proper  stabilizing  controller  parametrization,  it  is 
essential  to  compute  the  fractional  factorizations  over  the  ring  of  proper  stable  rational 
matrices.  Nett  et.  al.  [5°]  proposed  a  very  convenient  state-space  method  for  this 
computation  in  1984.  The  computation  method  was  developed  based  on  the  observer  and 
regulator  theories. 

Later  in  1984,  Doyle  et.  al.  [2]  showed  that  the  proper  stabilizing  controller 
parametrization  can  be  realized  as  an  observer-based  controller  with  an  added  stable 
parameter  matrix.  In  1988  and  1989,  Glover  and  Doyle  [8]  and  Doyle,  et.  al.,  [7]  offered 
the  two-Riccati-equation  approach  to  solving  the  standard  H°°  optimization  problem.  This 
approach  characterizes  all  possible  stabilizing  suboptimal  H“  controllers  whose  order  is  not 
higher  than  that  of  the  generalized  plant  Nonetheless,  since  the  generalized  plant  model 
includes  the  original  plant  and  the  models  of  appropriate  weighting  functions,  this  order  of 
the  controller  is  likely  to  render  it  non-implementable.  Hence,  a  need  for  a  suitable 
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methodology  for  the  reduction  of  this  controller  arises. 

In  this  chapter  we  deals  with  the  controller  reduction  from  an  observer-based 
compensator  point  of  view.  In  essence,  the  following  important  properties  related  to  the 
observer-based  controller  parametrization  and  pole  placement  will  be  given:  (1)  The  poles 
of  the  closed-loop  system  with  the  observer-based  controller  parametrization  are  the 
regulator  poles,  the  observer  poles,  together  with  the  poles  of  the  added  stable  parameter 
matrix.  (2)  If  the  controller  is  realized  by  a  minimal  realization,  the  closed-loop  poles  will 
include  all  the  poles  of  the  added  stable  parameter  matrix  and  a  subset  of  the  regulator  and 
the  observer  poles. 

In  the  rest  of  this  section,  we  will  explain  the  notations  to  be  used  in  this  chapter 
and  briefly  review  the  concept  of  the  stabilizing  controller  parametrization.  We  will  begin 
by  listing  the  previous  results  about  the  controller  parametrizations  done  by  Youla  et.  al., 
Desoer  et.  al..  Nett  et  al.,  and  Doyle  et  al.,  and  present  some  important  properties  of  the 
observer-based  controller  parametrization. 


Throughout  this  chapter,  both  of  the  following  notations 


A 

B 

C 

D 

and  {  A,  B,  C,  D  } 


are  used  for  the  same  purpose  to  represent  a  state-space  realization  of  a  system  whose 
transfer  function  is  CCsI-A)'1!!  +  D.  The  sum,  A+B,  of  two  sets  A  with  p  elements  and  B 
with  q  elements  is  a  set  which  consists  of  all  elements  of  A  and  B.  A+B  has  p+q  elements. 
Assume  that  B  is  a  subset  of  A,  then  the  difference,  A-B,  will  consist  of  all  the  elements  of 
A  except  those  in  B.  A-B  has  p-q  elements. 


G(s) 

- 

u 

Q(s) 

1 

y 

Figure  9. 1 .  Block  diagram  of  a  typical  control  problem. 

The  concept  of  stabilizing  controller  parametrization  is  briefly  described  as  follows. 
Consider  the  block  diagram  in  Fig.  9.1  where  v  is  the  exogenous  input  vector  which  may 
consist  of  the  disturbances,  noises,  and  the  commands,  u  is  the  control  input  vector 
through  which  the  behavior  of  the  system  can  be  modified,  z  is  the  controlled  output  vector 
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which  is  composed  of  all  die  variables  to  be  controlled,  and  y  is  the  measured  output  vector 
which  consists  of  all  the  measurable  quantities  available  for  feedback.  The  plant  G(s)  is 
given  by 


Gn(s) 

G21(s) 

G.,(s)l 

‘  A 

Bi 

V 

G(s)  = 

12' 

Gn(s) 

c, 

C2 

Dn 

°21 

D12 

D22 

i 

The  objective  of  a  typical  control  problem  is  to  find  a  proper  controller  Q(s)  which 
stabilizes  the  closed-loop  system  and  the  H“  (or  H2)  norm  of  the  closed-loop  transfer 
function  matrix  <D(s)  from  v  to  z  is  minimized  The  first  step  to  solve  the  problem  is  to  find 
the  set  of  all  proper  controllers  which  make  the  closed-loop  system  internally  stable.  Then 
in  the  set  of  all  proper  stabilizing  controllers,  one  will  be  chosen  such  that  the  H~  (or  H2) 
norm  of  <&(s)  is  minimized 

The  stabilizing  controller  parametrization  we  are  interested  in  has  the  following  two 
properties:  (1)  All  die  proper  stabilizing  controllers  can  be  characterized  in  terms  of  a  proper 
stable  parameter  matrix  K(s)  and  the  closed-loop  system  is  internally  stable  if  and  only  if 
K(s)  is  stable.  (2)  The  transfer  function  matrix  d>(s)  is  a  simple  affine  function  of  the 
parameter  matrix  K(s). 

After  the  stabilizing  controller  parametrization,  the  above  control  problem  becomes 
that  of  finding  a  proper  stable  matrix  K(s)  such  that  IWNI^  (or  IWWI2  )  is  minimized.  Property 

(2)  of  the  last  paragraph  is  important  since  it  will  make  the  H*°  (or  H2)  optimization 
problem  easy  to  solve. 

9.1.1  Preliminaries 

The  previous  results  related  to  the  controller  parametrization  will  be  briefly 
reviewed  in  this  section.  The  following  theorem  was  originally  developed  by  Youla  et  al. 
[57]  and  later  modified  by  Desoer  et  al.  [58]. 

Theorem  9.1:  [57,58]  (Youla's  Controller  Parametrization) 

Consider  the  system  in  Fig.  9.1.  Assume  that  the  realization  in  (9.1)  is  minimal  and  the 
subsystem  G22(s)  is  stabilizable  and  detectable.  Let  M2(s),  N2(s),  X2(s),  Y2(s),  Mj(s), 
Nj(s),  Xj(s),  and  YJs)  be  proper  stable  rational  matrices  such  that 
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‘  M^s) 

N2(s)' 

X^s) 

-N,(s)  * 

I  o’ 

-Yj(s) 

Xj(s) 

Y2(s) 

■ 

Mj(s) 

« 

0  I 

and 

M2(s)4N2(s)  =  G22(s) 

Then  the  set  of  all  proper  stabilizing  controllers  can  be  described  as 


(9-2) 

(9-3) 


{  Q(s)  I  Q(s)  -  [  M,(s)K(s)  +  Y2(s)  ]  [  N,(5)K(s)  -  Xj(s)  ]  •* 

with  K(s)  proper  stable  and  I  Nj(.»)K(o°)  -  X2(°°)  I  *  0  }  (9-4) 

and  the  closed-loop  transfer  function  matrix  <P(s)from  v  to  z  is  an  affine  function  of  the 
parameter  matrix  K(s), 


<Ks)  =  Gn(s)  -  G12(s)Y2(s)M2(s)G21(s)  -  G12(s)M1(s)  K(s)  Mj(s)G21(s)  (9-5) 


To  use  Theorem  9.1,  we  need  to  construct  the  proper  stable  rational  matrices  in  (9- 
2)  and  (9-3).  Nett  et  al.  [59]  proposed  a  convenient  state-space  approach  for  this 
construction.  That  is,  the  following  realizations 


‘  M2(s)  N2(s) 

-Yj(s)  Xj(s) 


X^s) 

-Nj(s) 

Y2(s) 

Mj(s) 

a+hc2 

H 

b2+hd22 

I 

D22 

-F 

0 

I 

a+b2f 

H 

®2  ' 

-(C2+D22F) 

I 

-D22 

F 

0 

I 

(9-6a) 


(9-6b) 


are  proper  stable  and  satisfy  (9-2)  and  (9-3)  where  F  and  H  can  be  arbitrarily  chosen  such 
that  A+B2F  and  A+HC2  are  stable. 

Doyle  et  al.  [2]  showed  that  if  (9-6a)  and  (9-6b)  are  used  to  realize  the  proper 
stable  rational  matrices  in  (9-2)  and  (9-3)  and  let 
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J(S) 


Jn(s)  J12(S) 
J21(S)  J22(S) 


A+B2F+HC2+HD22F 

-H 

-(B2+HD22)' 

F 

0 

-I 

-(C2+D22F) 

I 

D22 

(9-7) 


then  the  set  of  proper  stabilizing  controllers  described  in  Theorem  9.1  will  have  a  structure 
as  that  shown  in  Fig.  9.2. 


u 


J(s) 


K(s) 


u. 


with  K(s)  proper  stable  and  I-D^KO**)  invertible, 
figure  9.2.  Structure  of  stabilizing  controller  parametrization. 

Replace  the  controller  Q(s)  in  Fig.  9.1  by  the  structure  of  Fig.  9.2,  then  the  closed- 
loop  system  can  be  redrawn  as  that  shown  in  Fig.  9.3. 


u, 


- n 

1  1 

K(s) 

>  z 


Figure  9.3.  The  closed-loop  system  in  terms  of  a  parameter  matrix  K(s). 


In  Fig.9.3,  the  open-loop  transfer  function  matrix  from  to  y,  T22(s),  is  zero. 
Therefore,  the  closed-loop  transfer  function  matrix  from  v  to  z,  i.e.,  <I>(s),  is  a  simple 
affine  function  of  the  parameter  matrix  K(s).  That  is, 

<D(s)  =  Tn(s)  +  T12(s)  K(s)  T21(s)  (9-8) 

where  the  realizations  of  Tn(s),  TJ2(s),  T21(s)  are  given  by 


a+b2f 

-b2f 

B1 

Tu(s)  = 

0 

a+hc2 

Bi+HD2i 

- 

C1+Di2F 

1 

-D12F 

Dn 

i 
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v*) 


a+b2f 

V 

c1+d12f 

d,2. 

T21(s)  - 


A+HCj 

Bj+HDji 

C2 

D2i 

(9-9b) 


(9- 9c) 


Doyle  et  al.  [2]  also  pointed  out  that  the  structure  of  the  stabilizing  controller 
parametrization  in  Fig.  9.2  can  be  realized  as  an  observer-based  controller  with  an  added 
stable  dynamics  K(s).  The  realization  is  shown  in  Fig.  9.4. 


with  K(s)  proper  stable  and  I-D^KO**)  invertible 
Figure  9.4.  The  observer-based  controller  parametrization. 


Note  that  in  Fig.  9.4  the  block  diagram  inside  the  dotted-line  box  is  the  well-known 
full-order  observer-based  controller  [17]. 

9.1.2  Main  results 

In  Fig.  9.1,  the  internal  stability  of  the  closed-loop  system  depends  only  on  G^s) 
and  Q(s),  i.e.,  the  interconnected  system  shown  in  Fig.  9.5. 
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Figure  9.5.  Equivalent  system  to  Fig.  9.2  for  internal  stability. 

In  this  section,  the  controller  Q(s)  in  Fig.  9.5  is  replaced  by  the  block  diagram  of 
Fig.  9.4  which  is  the  observer-based  controller  with  an  added  dynamics  K(s). 

In  the  following  theorem  we  will  show  that  the  poles  of  the  closed-loop  sy 

with  the  observer-based  controller  parametrization  in  Fig.  9.4  are  the  regulator  poles  (l,. 
eigenvalues  of  A+BjF),  the  observer  poles  (the  eigenvalues  of  A+HC2),  together  with  the 

poles  of  the  parameter  matrix  K(s).  In  the  design  of  the  observer-based  controller,  F  and  H 
are  chosen  such  that  the  eigenvalues  of  A+B2F  and  A+HC2  are  stable.  Therefore,  the 

closed-loop  system  is  internally  stable  if  and  only  if  the  parameter  matrix  K(s)  is  proper 
stable.  The  proof  is  quite  straightforward  and  is  done  completely  in  the  state  space  without 
referring  to  the  derivations  used  by  Youla  et  aL,  Desoer  et  aL,  and  Doyle  et  aL. 

Theorem  92:  (Observer-based  Controller  Parametrization) 

Consider  the  closed-loop  system  in  Fig.  95.  Assume  that  G22(s)  =  {A,  B2,  C2,  D22) 

with  order  n  is  stabilizable  and  detectable  and  the  controller  Q(s)  is  replaced  by  the 
observer-based  controller  with  an  added  m-th  order  dynamics  K(s)  as  shown  in  Fig.  9.4. 
Then  the  set  of  the  closed-loop  poles  is  composed  of  the  n  eigenvalues  of  A+B the  n 

eigenvalues  ofA+HC2,  and  the  m  poles  of  the  added  dynamics  K(s).  That  is,  the  set  of  the 


closed-loop  poles  is 

^closed-loop  ~  ^regulator  +  ^observer  +  ^K(s)  (9- 10a) 

where 

^regulator  =  1  n  eigenvalues  of  A+B2F }  (9-10b) 

Pobg^vCT  ={  n  eigenvalues  of  A+HC2 )  (9- 10c) 

and 

?K(s)  =  {  m  poles  of  K(s) }  (9-10d) 

Proof:  The  dynamic  equations  of  G^Cs)  is  given  by 

x  *  Ax  +  B2u  (9-lla) 

y  =  CjX  +  D^u  (9-llb) 
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The  dynamic  equations  of  the  observer-based  controller  in  Fig.  9.4,  Le.,  the  block  diagram 
inside  the  dotted-line  box,  can  be  written  as  follows. 


=  (A+BjF+HCj+HD^F)  x  +  [*H  •(B2+HD22)] 


ul  f  F 

A 

0  -I' 

V 

y  =  [-(q^ 

X  + 

1  D22 
■  « 

U2 
*  • 

(9- 12a) 


(9- 12b) 


Assume  that  the  added  dynamics  K(s)  is  described  by  the  following  minimal 
realization 


k  =  Ak  +  By 
i*2  =  Ck  +  D  y 


(9- 13a) 
(9-13b) 


The  controller  Q(s)  is  just  a  combination  of  (9-12)  and  (9-13).  From  (9-12)  and  (9-13),  we 
have  the  dynamic  equations  of  the  controller  Q(s)  as  follows. 


“ll 

“12 

tt21 

tt22 

r  . 

a' 

X 

p.1 

+ 

k 

<S 

CO. 

(9- 14a) 


»  -  [v,  Tj] 


+  5y 


where 


Pj  -  -  H  -  (B^HD^)  (I-DD^)  D 


P2  =  B  +  B  D^a-DD^  'D 
Y,  -  F  +  a-DD^  'D  (C2+D22F) 

Yj  -  -  a-ray'c 

an  =  A  +  HC2  +  (B2+HD22)  Yj 
=  A  +  B2F  -  pj  (C^D^F) 


(9- 14b) 

(9- 15a) 
(9- 15b) 
(9- 15c) 

(9-15d) 

(9-15e) 

(9-15f) 
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“l2  *  (B2+HD22)  \  (9*15g) 

“21  -  -PjCCj+D^F)  (9-15h) 

022  *  A  -  B  D22  Y2  (9-15i) 

5  -  -a-DD22)'1D  (9-15j) 


Now,  combine  the  dynamic  equations  of  the  controller  Q(s)  and  the  plant  G^Cs), 

i.e.,  equations  (9-14)  and  (9-11).  Then  the  state  equation  of  the  closed-loop  system  can  be 
obtained  as  follows 


(9-16) 


where 


#n  =  A  +  B28(I-D225)-1C2 

(9- 17a) 

(9- 17b) 

#U  =  +  1 1^22  y2 

(9- 17c) 

*21  = 

(9-17d) 

*22  ~  an  +  ^1  (I  ~  ^22^  ^  ^22^1 

(9-17e) 

*23  ’  “i2  +  3,(I-D225)‘,D22T2 

(9-17f) 

*31  =  MI-DijSV'C, 

(9-17g) 

#32-a2I  +  P2(I-DaS)1D22r, 

(9-17h) 

#33  “  a22  +  P2(I-D22S)1D22Tf2 

(9-17i) 

Additionally,  we  make  the  following  observations,  which  can  be  easily  shown  by 
directly  manipulating  the  equations  (9-17)  and  (9-15). 


Observation.!: 

fH 

1 

II 

es 

«r 

(9- 18a) 

Observation  2: 

#11+#12  =  #21+#22  =  A+B2F 

(9-18b) 

Observation  3: 

*23  =  *13 

(9-18c) 

Observation.^: 

-#12+  #22  =  A  +  HC2 

(9-18d) 
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(9- 18c) 


Observation  $:  #33  »  a 


Define  the  observer  reconstruction  error 

_  A 

X  *  X  -  X 

then  we  have 


(9-19) 


X 

#11+ #12 

-#12 

X 

= 

#11+  #12  ‘  #21  '  #22 

*#12+#22 

k 

#31+  #M 

* 

"#32 

A  + 

b2f  *#12 

#13' 

X 

= 

0  A  +  HCj 

0 

X 

0  -#32 

A 

• 

k 

which  is  related  to  (9-16)  by  the  similarity  transformation 
I  0  O' 

I  -I  0 

0  0  I 


x 


33 


k 


(9-20) 


(9-21) 


Therefore,  from  the  structure  of  the  matrix  in  (9-20),  the  characteristic  values  of  die 

closed-loop  system  are  those  of  A+BjF,  A+HC2,  and  A .  This  completes  the  proof  of 
Theorem  92. 

It  is  well  known  that  in  the  observer-based  controller  design  the  closed-loop  poles 
are  the  regulator  poles  (the  eigenvalues  of  A+BjF)  and  the  observer  poles  (the  eigenvalues 
of  A+HC2)  [17].  In  Theorem  9.2  we  just  showed  that  the  above  property  still  remains 

when  we  add  a  dynamics  K(s)  to  the  observer-based  controller  as  shown  in  Fig.  9.4.  The 
eigenvalues  of  A+B2F  and  A+HCj  are  still  parts  of  the  closed-loop  poles  after  we  add  K(s) 

to  the  controller.  Adding  K(s)  only  introduces  additional  poles  to  the  closed-loop  system 
and  the  added  closed-loop  poles  are  the  poles  of  K(s).  If  F  and  H  have  been  chosen  such 


that  A+B2F  and  A+HCj  arc  stable,  then  the  closed-loop  system  with  the  observer-based 

controller  parametrization  will  be  internally  stable  if  and  only  if  the  parameter  matrix  K(s)  is 
proper  stable.  From  Fig.  9.4,  it  is  easy  to  see  that  the  controller  Q(s)  is  proper  if  K(s)  is 
proper  and  I  -  D^KO*)  is  invertible. 

With  the  observer-based  controller  parametrization,  the  closed-loop  transfer 
function  matrix  from  v  to  z,  i.e.,  4>(s),  is  a  simple  affine  function  of  the  parameter  matrix 
K(s).  That  is, 

4Ks)  =  Tn(s)  +  T12(s)  K(s)  T21(s)  (9-22) 

where  Tn(s),  TJ2(s),  and  T21(s)  are  given  by  (9-9).  The  added  dynamics  K(s)  is  a  proper 
stable  rational  matrix  to  be  chosen  such  that  I-D^KO**)  is  invertible  and  <b(s)  has  some 

desired  performance.  No  matter  which  K(s)  is  to  be  selected,  we  always  have  clear  idea 
that  die  closed-loop  poles  will  be  the  eigenvalues  of  A+B2F  and  A+HC2  together  with  the 

poles  of  the  added  dynamics  K(s)  if  the  controller  is  realized  as  that  shown  in  Fig.  9.4. 

Assume  that  the  orders  of  the  plant  and  the  parameter  matrix  K(s)  are  n  and  m 
respectively.  If  the  controller  is  realized  as  that  shown  in  Fig.  9.4,  then  the  order  of  the 
controller  is  n+m  and  the  closed-loop  system  has  2n+m  poles  described  by  the  set  Pck)sed_ 
loop  in  (9-10).  The  realization  of  the  controller  in  Fig.  9.4  may  not  be  minimal.  Suppose  it 

is  not  and  there  are  r  poles  in  the  controller  either  uncontrollable  or  unobservable,  then  the 
controller  can  be  realized  by  a  minimal  realization  with  order  n+m-r  and  the  number  of 
closed-loop  poles  will  be  reduced  to  2n+m-r. 

In  the  following  theorem,  we  will  show  that  the  uncontrollable  or  unobservable 
poles  of  the  controller  of  Fig.  9.4  must  be  the  eigenvalues  of  A+B2F  or  A+HC2  and  the 

closed-loop  poles  will  always  include  all  m  poles  of  the  added  stable  dynamics  K(s)  no 
matter  the  realization  of  the  controller  is  minimal  or  not  If  r  pole-zero  cancellations  occur  in 
the  controller,  then  the  closed-loop  poles  will  include  m  poles  of  K(s),  and  2n-r 
eigenvalues  out  of  the  set  PreguUtor + Pobserver  was  defined  in  (9-10). 

Theorem  9.3:  Consider  the  closed-loop  system  in  Fig.  95.  Assume  that  G22(s)  =  {A, 
B2,  C2,  D22}  with  order  n  is  stabilizable  and  detectable  and  the  controller  Q(s)  is  replaced 

by  the  observer-based  controller  with  an  added  m-th  order  dynamics  K(s)  as  shown  in  Fig. 
9.4.  A  minimal  realization  ofK(s)  is  given  by  (9-13).  Define  Pregulator ,  Pob^r^ pK(s) 

by  ( 9-10b ),  (9-1 0c),  and  (9-10d)  respectively  and  let 
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P _ _ * { the  controller  poles  which  are  either  uncontrollable  or  unobservable}  (9-23) 

Then 


P  +  P  3  p 

regulator  observer  removal 


(9-24) 


and  the  closed-loop  system  with  the  minimal  order  controller  will  have  a  set  of  poles 
described  by  the  following 


P  =  P  +  ( P  +  p  .  p 

closed-loop  with  min.  controller  K(»)  v  regulator  observer  removal 


)  (9-25) 


Proof:  Let  NG(s)DG(s)*1  with  deg  IDG(s)l  *  n  and  Dq(s)'xNq(s)  with  deg  IDq(s)I  =  n+m 
be  a  right  MFD  (matrix  fraction  description)  of  G^s)  and  a  left  MFD  of  Q(s)  respectively 
[19].  It  is  well  known  that  the  characteristic  polynomial  of  the  closed-loop  system  is 

^closed-loop^  *  1  Dq(s)Dg(s)  -  NQ(s)NG(s)  I  (9-26) 


Let  Lq(s)  be  a  greatest  common  left  divisor  of  Dq(s)  and  NQ(s).  That  is, 

Dq(s)  =  Lq(s)  Dq(s)  ,  Nq(s)  =  Lq(s)  Nq(s)  (9-27) 

A  A 

where  DQ(s)  and  NQ(s)  are  left  coprime.  It  is  easy  to  see  that  the  zeros  of  LQ(s)  are  the 
uncontrollable  or  the  unobservable  poles  of  die  controller  realization  in  Fig.  9.4.  That  is. 


{  zeros  of  LQ(s) }  =Premov>1  (9-28) 

Plug  (9-27)  into  (9-26),  we  have 

*cio*d-ioop<s>  =  1  Lq<s)  1  • 1  Dq(s)  Dg(s)  -  Nq(s)  NG(S)  I  (9-29) 


A  A 

The  zeros  of  I  Dq(s)  Dg(s)  -  Nq(s)  Ng(s)  I  are  the  poles  of  the  closed-loop  system 

with  a  minimal  controller  realization.  From  Theorem  9.2,  (9-28)  and  (9-29),  we  can  see 
that 

^closed-loop  with  min.  controller  *  ^removal  ~  ^regulator  ^observer  +  ^K(s)  ^ 

To  complete  the  proof,  we  need  to  show  that  Premovll  is  a  subset  of  Pregul4tor  + 
^observer  w^en  K(s)  is  minimal.  First,  assume  that  the  state-space  representation  (9-14)  of 
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the  controller  Q(s)  is  unobservable.  Then  by  PBH  test  [19],  there  exists  a  nonzero  vector  £ 
such  that 


(9-31a) 


(9-3  lb) 


for  some  eigenvalue  X  of  (9-14).  Note  that  it  is  the  eigenvalue  X  that  is  unobservable.  From 
(9-31b),  we  get 

“llSl  +  «12^2  *  <9'32a> 

which,  by  using  (9-15e)  and  (9-15g),  is  rewritten  as 

(A  +  HCj)*!  +  (B2  +  HD22)(y1^1+y242)  =  X$,  (9-32b) 

In  view  of  (9-31a),  the  above  equation  reduces  to 

(A  +  HC2Hj  =  X^  (9-32b) 

which  clearly  establishes  that  the  unobservable  eigenvalue  belongs  to  Pobserver- 


Proceeding  similarly,  it  can  be  shown  that  if  (9-14)  is  uncontrollable  then  the 
uncontrollable  eigenvalue  belongs  to  P  ^  • 

Note  in  the  above  development  that  ^  =  0  contradicts  the  minimality  assumption  of 
K(s).  Thus, 


P  ->  p 

observer  unobservable 

P  15  P 

regulator  uncontrollable 


(9-33a) 

(9-33b) 


where  Punobservabie  *s  all  the  unobservable  poles  of  the  controller  Q(s). 

^uncontrollable  *s  similarly.  This  completes  the  proof  of  Theorem  9.3. 

The  results  of  die  current  section  can  be  summarized  as  follows.  The  poles  of  the 
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closed-loop  system  with  the  observer-based  controller  parametrization  shown  in  Fig.  9.4 
can  be  classified  into  three  groups  and  each  group  of  poles  can  be  independently 

determined.  These  three  groups  of  poles  are  the  regulator  poles  (the  eigenvalues  of 
A+B2F),  the  observer  poles  (the  eigenvalues  of  A+HC2),  and  the  poles  of  the  added 

dynamics  K(s).  F,  H,  and  K(s)  are  free  parameters  to  be  chosen  such  that  the  closed-loop 
transfer  function  matrix  4>(s)  has  some  optimal  performance  subject  to  the  following 
constraints:  A+B2F  and  A+HC2  are  stable  and  K(s)  is  proper  stable  with  I-D22K(«») 

invertible. 

If  the  realization  of  the  controller  in  Fig.  9.4  is  not  minimal,  then  the  uncontrollable 
and/or  unobservable  controller  poles  can  be  removed  and  the  order  of  the  controller  is 
minimized.  The  set  of  these  removable  controller  poles  is  a  subset  of  the  regulator  and  the 
observer  poles.  The  poles  of  the  closed-loop  system  with  the  minimal  order  controller  will 
include  all  the  poles  of  the  parameter  matrix  K(s)  and  some  of  the  regulator  and  the 
observer  poles  which  are  not  the  removable  controller  poles. 


9.2  Controller  Reduction 

Motivated  by  the  results  of  the  previous  section,  we  pose  the  controller  reduction 
problem  as  follows:  given  die  plant  P(s)  as 


Xp=ApXp  +  B,ju  ;  xp eR  p;u eR 

Zp=CpXp  +  Dpu  iZpe/?1 

and  the  observer-based  controller  J(s)  as 

xj  =  ApXj  +  Bpu  +  Bj(Zp-zp  ;  Xj  eR*v 

zj  =  CpXj  +  DpU 
U  =  U; 


(9-34) 


(9-35) 


uj  -  cj*j 

obtain  a  reduced  order  controller  G(s)  whose  order  ng  is  required  to  be  less  than  that  of 
Ks): 


VVj  +  Vp  »  * 

u  =  CgXg  +  DgZp 


(9-36) 


Clearly,  if  the  above  realization  for  J(s)  in  non-minimal,  one  could  obtain  its 
minimal  realization  and  use  that  as  the  required  compensator  G(s).  Several  authors  have 
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been  motivated  by  this  approach;  for  example  see  references  [32,36,37]. 

We  take  a  different  approach,  which  is  motivated  from  the  observer-based 
controller  parametrization  of  the  previous  section.  Briefly,  our  approach  adopts  these  steps: 


1 .  Consider  the  control  u  as  consisting  of  two  components 


U  =  Uj  +  Uk 


(9-37) 


2.  Let  Uk  be  generated  from  an  additional  dynamic  parameter  K(s)  inserted 
between  the  measurement  residual  (Zp-Zj)  and  the  control  Uk  as  given  below. 

xk=  Akxk+Bk(zp-zi)  ;  xk€tf"k 

uk  =  CjjXfc  +  D^Zp-zp 

where  the  parameters  {Ak,  Bi,  Ck,  Dk)  are  yet  to  be  selected  -  see  Step  4 
below.  Note  that  we  are  increasing  the  order  of  the  overall  controller  by 
appending  the  additional  dynamics  K(s). 


3 .  Now,  the  overall  controller  Q(s),  consisting  of  both  the  given  controller  J(s) 
and  the  additional  dynamics  K(s),  can  be  realized  as 

(9-3 

XqsAqXq  +  BqZp  ;  Xq€/?n<«  ;  nq  =  n?+nk  v 

U  ~  CqXq  +  DqZp 


4.  Choose  the  parameters  { Ak,  Bk,  Ck,  DkJ  so  that  Q(s)  is  non-minimal.  Obtain 

tiie  reduced  controller  G(s)  as  the  minimal  realization  of  Q(s). 


The  purpose  of  this  study  is  to  show  that  the  state-space  parameters  {Ak,  Bk,  Ck, 
Dk}  of  K(s)  can  be  selected  to  render  Q(s)  non-minimal,  and  to  render  the  minimal  order  of 
Q(s)  to  be  less  than  np. 

We  first  rewrite  a  few  properties  of  observer-based  controllers  from  last  section  in 
tiie  notation  of  the  current  section. 


Theorem  92  :  The  eigenvalues  Ac/i  of  the  overall  closed  loop  system  are  given  by 

Ac/1=AruA0uAk  (9-40) 

where 

At  =  A(Ap+BpCj)  =  eigenvalues  of  the  regulator, 

Ao  =  A(Ap-BjCp)  =  eigenvalues  of  the  observer,  and 
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Ak  =  A(Afc)  =  eigenvalues  of  the  parameter  K(s). 

Notice  that  the  addition  of  K(s)  does  not  change  the  eigenvalues  of  the  closed  loop 
system;  it  only  adds  additional  eigenvalues.  Now,  if  the  controller  Q(s)  is  non-minimal  then 
there  are  some  eigenvalues  which  are  either  unobservable  and/or  uncontrollable,  and  these 
eigenvalues  can  be  removed  from  Q(s)  to  obtain  a  smaller  (minimal)  order  controller  G(s). 
When  the  state  space  realization  of  Kfs)  is  assumed  to  be  minimal,  the  following  result 
ensues. 


Theorem  9.3  :  Let  A*,  and  Auc  be  the  unobservable  and  uncontrollable  eigenvalues  of 
Q(s),  respectively.  Then 

Att0cA0&AatcAt  (9-41) 


Remark:  The  (2np  -  nuo  -  nuc  +  n^)  eigenvalues  of  the  reduced  closed  loop  system  will  be 
die  remaining  eigenvalues  of  die  regulator  and  the  observer  together  with  the  eigenvalues  of 
K(s).  Since  the  regulator  and  the  observer  are  normally  designed  to  be  stable,  the  stability 
of  the  reduced  closed  loon  system  is  therefore  guaranteed. 


Based  on  these  theorems,  we  present  the  following  main  theorem,  whose  proof  is 
outlined  in  the  Appendix. 


Theorem  9.4:  Let  Auo  and  Auc  be  some  eigenvalues  of  Q(s)  that  are  to  be  made 
unobserva  z  and  uncontrollable,  respectively.  Then  the  parameters  of  Kfs)  to  achieve  this 
must  satisfy 


Ak 

Bk" 

- 1 

* 

_ 1 

l 

§ 

_Ck 

Dk_ 

l~CV¥o_ 

- 1 

O 

&• 

(9-42a) 


Ak  Bk 
Ck  Dk 


for  some  *F*  and  <Pt ,  and  where  F0  and  <Pr  satisfy 


‘VVBpCp  =  AA 

(Ap-BjC^o  = 


(9-42b) 


(9-43a) 

(9-43b) 


The  order  ng  of  the  minimal  realization  of  Q(s)  is  then  given  by 
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n  =  n  -  n  -n  +  n,. . 

g  p  uo  uc  Tc 


(9-44) 


It  is  possible  to  select  K(s)  such  that  ng  is  less  than  np,  as  shown  in  the  example  below. 

In  what  follows  we  study  a  special  case  of  this  problem  when  K(s)  is  static  (n^  = 
0),  i.e.,  the  constant  matrix  D*  is  the  only  parameter  that  describes  K(s).  The  order  of  the 
reduced  controller  would  then  be  ng  =  np  -nuo  -n^  .  Results  presented  below  are  easily 
extended  for  the  more  general  case  of  dynamic  K(s). 

Corollary  9.1:  The  necessary  and  sufficient  condition  for  the  parameter  Dt  (i.e.  static 
K(s))  to  achieve  this  is  the  satisfaction  of  the  following  equations: 

Dj-Cp'Po  =  Cj'Fo,  and  <&rBpI>k  =  -4>rBj,  (9-45) 


subject  to  (9-43a)  and  (9-43b). 

Though  not  necessary,  <br  and  can  be  seen  as  consisting  of  some  left 
eigenvectors  of  the  regulator  and  some  right  eigenvectors  of  the  observer,  respectively. 
Note  that  if  either  nUo  or  nUc  is  equal  to  zero  then  only  one  of  the  above  equations  need  to 
be  satisfied.  For  instance,  suppose  that  the  given  controller  J(s)  is  not  completely 
observable,  which  may  happen  as  in  the  case  of  LQG  controllers.  Then  there  exists  a  % 
such  that  Cj'Fo  =  0,  and  hence  a  =  0  will  remove  these  unobservable  eigenvalues.  This 
result,  along  with  those  from  the  previous  section,  establishes  the  following: 

Corollary  9.2:  If  a  given  controller  J(s)  is  non-minimal,  then  the  closed  loop  system 
using  a  minimal  realization  ofJ(s)  as  its  reduced  controller  would  contain  only  the 
observable  eigenvalues  of  the  observer  and  the  controllable  eigenvalues  of  the  regulator. 

The  above  results  are  pertaining  only  to  the  stability  of  the  closed  loop  system.  An 
equally  important  issue  in  controller  reduction  is  the  performance  of  the  closed  loop 
system.  This  issue  can  be  addressed,  at  least  in  a  subopdmal  sense,  if  one  could  determine 
the  eigenvalues  that  are  least  significant  in  the  performance.  The  Modal  Cost  Analysis  [60] 
offers  a  suitable  methodology  for  identifying  such  eigenvalues.  However,  only  those 
eigenvalues  that  allow  the  satisfaction  of  the  above  equations  can  be  removed.  The 
following  result  assists  in  identifying  such  removable  eigenvalues. 

Theorem  9.5: 

(a)  When  n„c  =  0  and  n„o>0,an  order-reducing  D*  exists  if  and  only  if 
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KerneKCyPo)  c  Kemel(Ci'P0),  (9-46a) 

with  the  general  solution 

+  Z  [  HC^C^oT  ],  (9-46b) 

(b)  When  nMO  =  0  and  nHC>0,an  order-reducing  D  k  exists  if  and  only  if 

Image(<I>rBj)  c  Image(<&rBp),  (9-47a) 

with  the  general  solution 

Dk = -  ]  Z .  (9-47b) 


(c)  Under  the  assumption  of  =  when  both  nuo  >  0  and  nuc  >  0,  a  Dk 

satisfying  the  equations  of  Corollary  9.1  exists  if  and  only  if 

<*>^0  =  0, 

along  with  the  satisfaction  cf(9-46a)  and  (9-47 a),  with  the  general  solution 
Dk=^<&1Bppx>1Bj + cj^0(cp'p0)-- 

+  [  1 Z I  HCp^oKCp^o)'] 

where  ( •/“  denotes  generalized  inverse,  andZ  is  arbitrary 

Under  the  assumption  of  A^  n  Auc  =  0,thc  above  theorem  suggests  this  procedure 
for  determining  the  removable  eigenvalues:  Compute  the  complete  set  of  the  left 
eigenvectors  of  the  regulator,  and  the  right  eigenvectors  of  the  observer,  denote  the  first  set 
by  d>  and  the  second  by  'P;  use  the  zero  elements  of  die  product  [4>Y]  to  identify  the 
appropriate  left  and  right  eigenvectors;  collect  these  left  eigenvectors  to  form  $r  and 
die  right  eigenvectors  to  form  *F0.  Then,  if  the  conditions  (9-46a)  and  (9-47a)  are  also 
satisfied,  the  feed-through  loop  determined  according  to  Part  (c)  of  Theorem  9.5  will  yield 
a  reduced  controller  of  order  (np-nU0-nUc)-  The  absence  of  any  zero  elements  in 
implies  that  a  feed-through  loop  does  not  exist  which  would  simultaneously  render  some 
observer  eigenvalues  unobservable  and  some  regulator  eigenvalues  uncontrollable.  (A 
dynamic  K(s)  may  however  be  constructed  as  shown  later  in  this  report).  It  may,  however, 
be  possible  to  remove  some  eigenvalues  of  only  the  observer  or  only  the  controller.  This  is 
because  Part  (c)  of  the  above  theorem  is  pertaining  to  the  simultaneous  satisfaction  of  the 
two  equations  of  Corollary  9.1,  which  is  not  necessary  when  either  nuo  or  nuc  is  equal  to 


(9-48a) 


(9-48b) 


zero.  In  this  case  we  have  Parts  (a)  and  (b)  of  the  theorem. 

An  Example: 


Let  die  plant  parameters  be 


V 


0  0 
1  0 


B. 


-  i:i  • 


C„=  I1  *]  ;  V  o. 


and  let  the  regulator  gain  Cj  and  die  observer  gain  B  j  be 

S  -  [ 31  ;  - 


The  order  of  die  controller  J(s)  is  2,  and  the  eigenvalues  of  the  closed  loop  system 
consisting  of  P(s)  and  J(s)  are  those  of  the  regulator  Af  and  the  observer  Ao  where,  for  this 
example, 

A,  =  {-1,-3}  and  Ao={-2,-2}. 


Suppose  that  we  wish  to  reduce  the  order  of  the  controller  by  removing  an  -2 
eigenvalue  of  the  observer  and  the  -1  eigenvalue  of  the  regulator.  This  can  be  achieved  by 
introducing  K(s)  whose  parameters  are 

Ak  =  Ak  ;  Bk  =  Ak+2;  Ck  =  Ak+l  ;  Dk  =  A]c  +  6, 
for  any  arbitrary,  but  stable  A*.  The  corresponding  realization  of  Q(s)  has  the  parameters 


V2 

V1 

-V1 

-V2' 

A  = 

0 

0 

;  B  = 

0 

q 

-V2 

-V2 

\ 

’  q 

Ak+2 

*  a 

*  m  »  «■ 

Cq=  [V2  V3  -V1];  Dq=  [A-6] 


which  has  an  unobservable  eigenvalue  of  -2  and  an  uncontrolable  eigenvalue  of  -1,  as 
required.  The  reduced  order  controller  G(s),  which  is  the  minimal  realization  of  Q(s),  is 
then  given  by  die  parameters 

Ag  =  2Aj+l  ;  Bg  =  -(Ak+2)  ;  Cg  =  2Ak+3 
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D$-  -<V<>) 

By  Theorem  9.3  above,  the  eigenvalues  of  the  closed  loop  system  consisting  of  P(s)  and 
the  reduced  controller  G(s)  are  { -2,  -3,  A*}. 

For  the  same  example,  the  following  table  shows  what  can  be  achieved  by  static 
K(s),  that  is  by  an  appropriate  choice  of  Dfc. 


msm 

iBT’TI 

Q 

ElfiSfilMH 

flHKSflHHNMI 

kls 

■ 

SS3HH 

MmM 

■■ 

&£)  mmmm 

We  now  address  the  case  when  the  parameter  K(s),  between  the  measurement 
residual  (zp-zj)  and  the  control  uk  is  dynamic.  Appending  additional  dynamics  K(s)  to  the 
original  controller  J(s),  increases  die  dimension  of  the  overall  controller  Q(s)  to  nq=np+nk, 
where  n«j  is  the  dimension  of  Q(s)  and  nk  is  that  of  K(s).  Recalling  the  properties  of 
observer-based  controllers  with  respect  to  the  closed  loop  poles,  the  addition  of  the 
parameter  K(s)  does  not  change  the  eigenvalues  of  the  closed  loop  system  but  just  adds  to 
them.  This  implies  that  die  internal  stability  of  the  closed  loop  system  depends  on  the 
stability  of  K(s).  Thus  by  choosing  a  proper  stable  K(s)  one  can  guarantee  the  stability  of 
the  closed  loop  system.  The  necessity  of  satisfying  the  orthogonality  condition  between  die 
left  regulator  and  right  observer  eigenvectors  corresponding  to  the  eigenvalues  to  be  made 
uncontrollable  and  unobservable  respectively,  makes  it  difficult  to  determine  a  static  feed 
through  loop  that  would  yield  a  reduced  order  controller.  However,  controller  reduction  by 
observer  based  controller  parametrization  may  be  achieved  when  the  parameter  matrix, 
K(s),  is  dynamic.  This  section  extends  the  results  obtained  for  the  existence  of  an  order 
reducing  static  feed-through  loop  to  the  dynamic  case. 

Theorem  9.4  states  that  the  necessary  and  sufficient  condition  for  the  existence  of 
the  parameters  of  K(s)  to  achieve  non-minimality  of  the  overall  controller  Q(s)  is  the 
simultaneous  satisfaction  of  equations  (9-42a)  and  (9-42b).  The  following  theorem  states  a 
necessary  and  sufficient  condition  for  the  existence  of  the  parameters  of  K(s)  that  will 
simultaneously  satisfy  equations  (9-42a)  and  (9-42b). 

Theorem  9.6:  Under  the  assumption  that  Auc  O  Auo  =  <f>,  the  parameters  of  K(s) 
simultaneously  satisfying  equations  (9-42a)  and  ( 9-42b )  exist  if  and  only  if 

+0*^  =  0  (9-49) 
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where  and  ¥i  are  arbitrary. 

Proof :  Assume  A,*,  Pi  A^,  =  $.  Then  premultiplying  Eqn.(9-42a)  by  [d^  d^p],  yields 

^i^uo  -*PvCpo  =  <^4^  +  <I>rBiCp'P0  (9-50a) 

Adding  and  subtracting  Ap,  and  using  Eqns.(9-43a)and  (9-43b)  results  in 

WV o  •  KziWo  -Hlhc^kl  =  0  (9-50b) 

The  solution  to  equation  (9-50b)  under  the  hypothesis  of  die  theorem  yields  (9-49). 

Recall,  that  the  existence  of  an  order  reducing  static  feed-through  loop  is  guaranteed 
under  the  satisfaction  of  an  orthogonality  condition  between  the  left  regulator  and  right 
observer  eigenvectors  associated  with  the  removable  eigenvalues.  This  necessary  and 
sufficient  condition  is  restrictive  in  the  sense  that,  rarely  can  one  find  systems  which 
would  satisfy  this  orthogonality  condition.  However,  the  necessary  and  sufficient  condition 
stated  according  to  Theorem  9.6.  can  almost  always  be  satisfied  indicating  that  one  could 
almost  always  determine  a  dynamic  feed  through  loop  that  would  yield  a  reduced  order 
controller.  This  is  attributed  to  die  arbitrariness  of  and  'F*. 

Theorem  9.6  indicates  that  it  is  possible  to  simultaneously  render  some  observer 
eigenvalues  unobservable  and  some  regulator  eigenvalues  uncontrollable.  Having 
established  the  existence  of  dynamic  feed  through  loops  that  would  yield  reduced  order 
controllers,  the  next  logical  step  would  be  the  computation  of  the  parameter  matrix  K(s). 
Theorem  9.7  presents  a  closed  form  solution  to  the  equations  (9-42a)  and  (9-42b).  This 
result  is  an  extension  of  the  case  when  K(s)  is  static,  and  can  be  proved  similarly. 

Theorem  9.7:  The  solution  {AkJBk,CkJ)j)to  the  necessary  and  sufficient  conditions 
equations  is 

Ak  ~  Up*  +  ELp'Bj  -  HpHp£j£p  +  P  -  &p'2Lp]Z[I  -  £p£Lp~]  (9-51) 

where 

r\Bki  [>Aol  r  \ ' 

Ak  =  .ciDk.  CjTci'i'"J  Cp=.-cvp». 

(9-52) 

Bp  =  [<t>k  4>,BP]  Bj=[V*k  -*^j] 
where  (•)'  denotes  generalized  inverses,  andZ  is  arbitrary. 
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The  family  of  dynamic  K(s)  constructed  according  to  Theorem  9.7,  will  yield  a 
reduced  order  controller  of  dimension 

n,  =  np  •  (nuo  +  nuc)  +  nk  (9-53) 

where  nk  is  the  order  of  die  parameter  matrix  K(s). 

From  die  properties  of  observer  based  controllers,  the  resulting  closed  loop  system  will 
always  retain  the  characteristic  values  of  the  parameter  matrix  K(s).  i.e  the  resulting 
closed  loop  system  will  be  stable  if  and  only  if  Ak  is  stable.  To  pursue  this  issue  consider 

A*  =  Li  £p‘  +  ftp"  fij  -  ftp’ftpCiCp*  (9-54) 

B*  =  [I  -  &p~&p]  ;  C*  *  [I  -  £p£p-]  (9-55) 

Thai  Eqn.  (9-51)  can  be  rewritten  as 

Ak  =  A*  +  B*  Z  C*  (9-56) 

Notice  that  equation  (9-56)  has  a  structure  similar  to  that  of  an  output  feedback 
problem,  where  Z  is  a  free  parameter  to  be  chosen  such  the  eigenvalues  of  Ak  are  stable. 
However,  since  Ak  is  only  a  submatrix  -  namely,  the  (1,1)  submatrix  -  of  A*,  the 
corresponding  submatrix  in  Eqn.(9-56)  need  be  considered  as  an  output  feedback  problem. 
Thus  a  the  free  parameter  Z  which  would  guarantee  a  stable  Ac  will  exist  if  and  only  if  the 
corresponding  submatrix  pair  of  (A*3*)  is  stabilizable,  and  the  corresponding  submatrix 
pair  of  (C*A*)  is  detectable.  More  stringent  conditions  under  which  one  could  guarantee 
die  stability  of  K(s)  are  yet  to  be  determined. 

Besides  stability,  performance  of  the  resulting  closed  loop  system  is  an  equally 
important  issue  in  controller  reduction  problems.  Both  Theorem  9.6,  and  9.7  seem  to 
indicate  that  the  choice  of  removable  eigenvalues  play  an  important  role  in  die  computation 
of  the  parameters  of  K(s).  One  could  employ  one  of  several  existing  schemes,  such  as 
balancing  [31]  and  component  cost  analysis  [45],  to  determine  the  ’’removable” 
eigenvalues.  Another  way  of  approaching  the  optimization  problem  is  as  follows:  one  could 
determine  an  optimal  controller  K*(s),  using  standard  parametric  optimization  approach. 
We  could  then  select  a  K(s)  which  would  best  approximate  the  optimal  controller  by  (say ) 
solving  a  suitable  minimum  norm  problem.  The  arbitrariness  of  Z  in  equation  (9-51)  allows 
us  some  flexibility  in  obtaining  a  stable  K(s)  that  would  result  in  a  reduced  order  controller 
which  would  best  approximate  the  full  order  controller  in  terms  of  a  specified  performance. 

9.3  Conclusions 
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It  has  been  shown  that  by  introducing  a  static  feed-through  loop  between  the 
measurement  residual  and  the  control  signal,  die  order  of  a  given  observer-based  controller 
can  be  reduced  while  guaranteeing  the  stability  of  the  closed  loop  system.  This  is  achieved 
by  making  some  eigenvalues  of  the  controller  uncontrollable  and/or  unobservable.  The 
necessary  and  sufficient  conditions  for  the  existence  of  such  feed-through  loops  have  been 
characterized,  thereby  parametrizing  such  loops.  These  results  have  been  extended  to  the 
case  when  K(s)  is  dynamic,  however,  stability  of  the  resulting  closed  loop  system  is  still 
an  unresolved  issue.  With  the  parametrization  accomplished,  problems  can  now  be  posed 
to  address  the  performance  issues. 
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CHAPTER  10 


DOMINANT  MODES  OF  MECHANICAL  SYSTEMS 


Recall  from  Chapter  9  that  in  order  to  incorporate  the  parametrization  of  order- 
reducing  K(s)  in  a  design  procedure  requires  die  identification  of  poles  that  are  to  be 
rendered  uncontrollable  and  unobservable.  This  chapter  deals  with  identifying  modes  of  a 
mechanical  system  that  are  dominant  in  an  L2  sense.  Specifically,  it  considers  obtaining  a 
reduced  model  of  a  stable  mechanical  system  with  proportional  damping.  Such  systems  can 
be  conveniently  represented  in  modal  coordinates.  The  popular  scheme.  Modal  Cost 
Analysis  (MCA)  of  [61],  offers  a  simple  means  of  identifying  dominant  modes  for 
retention  in  the  reduced  model.  In  MCA,  the  dominance  is  measured  via  the  modal  costs. 
Though  this  measure  leads  to  simple  computations,  it  does  not  exactly  reflect  the  more 
appropriate  model  error  which  is  the  L2  norm  of  the  output-error  between  the  full  and  the 
reduced  models.  Normally,  the  model  error  is  computed  after  the  reduced  model  is 
obtained,  since  it  is  believed  that,  in  general,  the  model  error  cannot  be  easily  computed  a 
priori.  The  main  thrust  of  this  note  is  to  point  out  that  the  model  error  can  also  be  calculated 
a  priori,  just  as  easily  as  the  modal  costs.  Hence  the  model  error  itself  can  be  used  to 
determine  the  dominant  modes.  Moreover,  the  simplicity  of  the  computations  do  not 
presume  any  special  properties  of  the  system,  such  as  small  damping,  orthogonal- 
symmetry,  etc.  The  development  presented  herein  can  be  seen  as  a  specialization  of  that  in 
[62],  to  mechanical  systems. 

10.1  Problem  Formulation 

Consider  a  time-invariant  mechanical  system,  described  in  its  physical  coordinates 
q,  given  below. 

M  q(t)  +  D  q(t)  +  K  q(t)  =  B  u(t) 
y(t)  =  Cq(t)+£q(t) 

where  q  €  Rn ,  u  €  R®,  y  e  Rk ;  u  is  assumed  to  be  a  Gaussian  white  noise  process  with 
unit  intensity.  Both  the  mass  matrix  M  and  the  stiffness  matrix  K  are  assumed  to  be 
symmetric  and  positive  definite.  The  dissipation  matrix  D  is  assumed  to  arise  from 
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proportional  damping  (not  necessarily  small)  such  that  the  above  system  is  asymptotically 
stable.  Under  these  assumptions,  the  model  (10*1)  can  be  equivalently  expressed  in  its 
modal  coordinates  as 


tii  +  2Ci<0itii  +  cofni  =  bju  ;  i  =  1,2 . n  (10-2a) 

y  =  S  (emi  +  fifli)  (10-2b) 

i-l 

where  (Oj  and  are  the  usual  i*  natural  frequency  and  the  corresponding  damping  ratio, 
respectively. 

When  n  is  large,  one  is  faced  with  the  problem  of  model  reduction  to  facilitate  a 
subsequent  analysis  and  control  design.  Quite  often,  a  reduced  model  of  the  system  (10- 
2&&b)  is  obtained  by  retaining  some  r  modes  of  the  total  n  modes  -  a  process  called  herein 
as  a  modal  reduction.  In  order  to  produce  an  acceptable  reduced  model,  the  question  that 
needs  to  be  answered  is  "which  r  of  the  n  inodes  should  be  retained  ?"  Using  the  notation 
Nn  for  the  integer  set  { l,2,...,n},  the  issue  in  modal  reduction  is  to  identify  an  r-element 

"reduction-set"  c  Nn,  such  that  the  following  model 

iji  +  2^i00jT|i  +  oofo  =  bju  ;  i €  Nxed  (10-3a) 

yr  =  X  (cffli  +  Ctfi)  (10-3b) 

ieN«i 


would  be  an  acceptable  reduced  model  Defining  the  error  e(t)  between  the  full  model  (10- 
2a&b)  and  the  reduced  model  (10-3a&b)  as  e(t)  =  y(t)  -  yr(t),  an  acceptable  reduced 
model  would  minimize  the  model  error  5V  defined  below. 


8V  = 


lim 

i-»~ 


E 


eT(o)e(a)  do 


(10-4) 


10.2  Dominant  Modes  by  Modal  Cost  Analysis 

For  an  arbitrary  system,  not  necessarily  a  mechanical  system,  the  model  error  8V 
associated  with  its  reduced  model  is  normally  computed  after  the  reduced  model  is 
obtained.  Consequently,  obtaining  an  acceptable  reduced  model  for  such  systems  becomes 
an  iterative  process.  In  order  to  simplify  the  reduction  process,  methods  have  been 
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proposed  to  employ  different  criteria  instead  of  (1(M).  One  such  method  is  die  Component 
Cost  Analysis  (CCA)  of  Skelton  [45],  The  CCA  attempts  to  use  die  cost  error  AV  defined 
below,  as  its  criterion  for  model  reduction. 


AV  =  V-Vr , 


where 


V  =  Sm  E  l 

tr*~  It 


yT(o)y(o)dc 


(10-5a) 


(10-5b) 


Vr  =  E  1 1  y?(o)yr<a)  doj. 


(10-5c) 


Now,  the  reduced  models  which  are  optimal  in  the  sense  of  minimizing  the  model 
error  SV,  are  known  to  satisfy  the  orthogonal  property  [46] 


(e,yr)  =  lim  E  |l  f 
Wo 


t-4~ 


cT(o)yr(o)  doj  =  0. 


(10-6) 


Since  the  model  error  6V  can  be  expressed  in  terms  of  the  cost  error  AV  and  the  inner 
product  (e,yr)  as  given  below 


8V  =  AV  -2(e,yr), 


(10-7) 


it  follows  that  the  cost  error  is  an  appropriate  criteria  to  use,  provided  the  resulting  reduced 
model  is  known  to  be  near  optimal,  if  not  optimal.  Though  the  computation  of  AV  is 
simpler  than  that  of  SV,  it  still  requires  the  availability  of  the  reduced  model.  In  view  of 
this,  the  CCA  employs  the  predicted  cost  error  A  V,  defined  below,  for  its  model  reduction 
decisions. 

AV  =  V-Vr  ,  (10-8) 

>■*>* 

where  Vr  is  a  prediction  of  Vr .  With  defined  as  the  (first  order)  state  representing  the  i* 
component,  Vr  is  computed  according  to 

Vr=  X  Vj  (10-9a) 

i  €  Njed 
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(10-9b) 


and  where  Vi ,  the  i*  component  cost,  is  defined  as  [45] 


Vi 


:1  Bm  p 

2  t-»- 


a[yT(g)y(c)3 

dxi(a) 


Xi(o)  do 


Since  the  calculation  of  the  component  costs  requires  only  the  full  model,  the  predicted  cost 
error  A  V  can  be  computed  before  a  reduced  model  is  obtained. 


The  application  of  die  CCA  to  models  represented  in  modal  coordinates,  such  as  in 
Eqs.  (10-2a  and  b),  is  called  the  Modal  Cost  Analysis  (MCA),  see  [61].  Therefore,  modal 
reduction  by  MCA  chooses  an  Nred  such  that  the  predicted  cost  error  AV  is  minimized. 

Except  under  special  cases,  there  is  no  guarantee  that  die  predicted  cost  error  A  V  equals  the 
model  error  5V.  It  turns  out  that,  for  any  choice  of  N^,  the  corresponding  5V  can  be 

computed  fust  as  easily  as  computing  AV.  Hence,  one  may  use  the  more  appropriate 
criterion  5V  for  choosing  N^.  These  details  are  developed  below. 


10.3  Closed  Form  Expression  for  Model  Error 

Provided  all  the  modes  are  observable  and  controllable,  truncation  of  any  mode 
would  affect  the  output  y.  To  determine  such  affect,  define 

Vi  -  Citii  +  cjfii  so  that  y  =  Vi  (10-10) 

ieN. 


and  call  v4  as  the  "i*  modal-output"  By  defining  the  "truncation-set"  =  Nn  -  Nred, 
the  following  expressions  result: 

c  =  X  ^  (10-1  la) 

ie  N« 


8V  = 


(10-1  lb) 


where 


ur..  _  In  c 


j*-  [  ViWj 


(o)Vj(c)do 


;  id  e  Nn 


(10-1 lc) 


In  the  spirit  of  CCA,  call  die  cost-correlation  between  the  modes  i  and  j.  Note  that  'P- 


127 


measures  the  correlation  between  the  i*  and  the  j*  modal-outputs,  and  depends  only  on  die 
i*  and  the  j*  modal  data.  The  expression  for  computing  S'-  for  all  i  j  e  Nn ,  is  given  in 

die  next  section.  Given  the  modal  data,  these  cost-correlations  can  be  computed  for  all  the 
modes,  and  the  model  error  5V  can  be  determined  for  any  r-element  reduction-set 
^i«d  c  Nn  from  the  expression  (10-1  lb).  Notice  that  the  model  error  can  be  computed  a 
priori  by  simple  summation  of  appropriate  S'-s. 

Clearly,  there  are  (n!  /  {r!(n— r)! }]  number  of  possible  reduction-sets,  of  which  one 
will  yield  the  smallest  model  error.  To  facilitate  the  selection  of  the  best  reduction-set, 
construct  an  (n,n)  "cost-correlation-matrix"  S'  whose  (ij)  element  equals  S'-,  and  define 

0  m 

an  operation  XX :  R(n,m)  — >  R  for  arbitrary  n  and  m  by  XX (Z)  =  XX  Z;j .  Then,  the 

i=i  j-i 

cost-correlation-matrix  possesses  these  properties  for  any  N^: 


1) 

5V 

s 

ZZWtJ 

(10- 12a) 

2) 

V 

= 

XXOF) 

(10-12b) 

3) 

= 

52^^) 

(10- 12c) 

4) 

(e.yr) 

s 

220Prelr) 

(10-12d) 

5) 

Vi 

= 

XX0Fs) 

(10-12e) 

where 

s 

[  S'-;  i  j  €  ],  an  (n-r,n-r)  submatrix  of  S', 

(10-13a) 

^red 

s 

[  S'jjj  i  j  e  Nred  ],  an  (r,r)  submatrix  of  S', 

(10-13b) 

^retr 

s 

[  S'-;  i  e  ,  j  €  N^],  an  (r,n-r)  submatrix  of  S',  (10-13c) 

*i 

= 

[  S'-;  j  €  Nn],  an  (l,n)  matrix  (i.e.,  i*  row)  for  all  i.(10-13d) 

Notice  that  the  fourth  property  of  the  cost-correlation-matrix  allows  one  to  determine  a 
suboptimality  index  [63]  associated  with  the  reduced  model.  This  index  can  be  used  to 
determine  if  the  reduced  model  could  be  improved  by  further  optimizations.  It  measures 
how  closely  the  orthogonality  condition  (10-6)  is  satisfied.  The  second  and  third  properties 
above  provide  means  of  computing  the  quadratic  cost  functions  V  and  Vr  associated  with 

the  full  and  the  reduced  models,  respectively.  Moreover,  the  three  errors,  namely,  the 


model  error  8V,  the  cost  error  AV,  and  the  predicted  cost  error  A  V ,  are  related  through  the 
property  (10-12d)  as  follows 

A V  «  8V  +  2 IIOF^) ,  (10- 14a) 

AV  *  8V  +  2Z0Frelr) .  (10- 14b) 

Hence,  unless  ££(¥„,,)  -  0,  modes  determined  as  dominant  based  upon  either  the  cost 
error  or  the  predicted  cost  error  may  not  yield  the  smallest  model  error.  Since  the 
expression  (10-1  lb)  is  available  for  the  computation  of  the  model  error  for  any  choice  of 
Nred,  dominance  of  modes  can  be  established  based  on  5V  itself.  The  disadvantage, 

however,  is  in  having  to  compute  [(n-t)2-l]  additions  [n!  /  (r!(n-r)! }]  times. 

Theorem  10.1:  The  expression  for  computing  Y^for  all  ij  e  Nn,  is  given  below. 

x'/  [  cjc,  -  -  4*! 1  +  T51*  (1°-15a) 

5ij  5ij  5ij 

where 

Xj.1  =  t  (10-15b) 

1  Ay  +  <n^5?j  -  2^it0i6ijAij 

Ajj  =  o>2  -  coj,  and  8y  =  2^  +  2£j0>j .  (10-15c) 


Proof:  Defining  the  i*  modal  states  as  xj  =  [  T\ .  r\ .  ],  it  follows  from  the  definition  of  the 
cost-correlation  and  the  modal-outputs  y.  that 


lim  Ev 

t->~  1 


L 


x[(o)CfCjxj(a)  da 


=  trace{(?CjXji} 


(10- 16a) 


where 


Xi(o)xJ(o)do|, 


(10-16b) 
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q  *  [cifii]. 


(10-16c) 


The  steady  state  correlation  matrix  between  the  i*  and  die  j*  modal  states  is  given  by 


(10- 17a) 


where 


Xjf  -  -X?/  »^x|/ 


X?;2  = 


c  c 


(10-17b) 


with  xfj1  as  defined  above  in  Eqn.(10-15b).  Substitution  of  the  Eqs.  (10-17a)  and  (10- 
17b)  into  Eqn.(lG-16a)  completes  the  proof. 

Expressions  similar  to  (10- 17a  &b)  can  be  found  in  several  references,  see  [64]  for 
example.  The  expression  (10- 15a)  is  valid  for  all  i  and  j,  regardless  of  whether  the 
frequencies  are  repeated  or  whether  the  damping  is  small.  Of  course,  the  expression 
simplifies  undo-  special  cases,  as  shown  below. 


a) 

If  to,  =  C0j,  then  ^  ^  ^  {cJcj-HC^^} 

(10-18a) 

b) 

If  tOj  *  ©j,  and  both  Q,  ^  ->  0,  then  'Fy  -4  } 

^ij 

(10-18b) 

c) 

If  to,  =  C0j,  and  both  ->  0,  then  *F;j  -4  — ^  {cfo  +  co&fe } 

4^i(0i 

(10- 18c) 

Note  by  the  special  case  (b)  above  that  for  a  lightly  damped  system  (such  as  a  flexible 
space  structure)  with  distinct  frequencies  and  with  either  (i)  bjbj  =  0  for  all  i  and  j,  or  (ii) 
fijcj  =  0  for  all  i  and  j,  one  obtains  ->  0  for  all  j  *  i,  and  for  all  i.  It  therefore 

ys 

follows  from  Eqns.(10-8,10-9a,10-llb,10-12b,10-12e)  that  Vi  =  W  8V  =  X 

■  _  i  €  N„ 

and  8V  =  AV.  Hence,  for  systems  satisfying  these  assumptions,  reduced  models 
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obtained  by  MCA,  which  are  guaranteed  to  minimize  the  predicted  cost  error  AV,  also 
minimiy^.  the  model  error  5V. 


10.4  Examples 

Two  examples  are  provided  now  to  show  that  reduced  models  generated  as 
suggested  in  this  note  will  be  at  least  as  good  as,  if  not  better  than,  those  produced  by 
Modal  Cost  Analysis.  In  both  the  examples,  we  begin  with  five  modes  and  produce 
reduced  models  of  different  orders.  The  reduced  models  are  evaluated  in  terms  of  their 
model  errors,  and  the  results  are  tabulated  below.  The  two  examples  differ  only  in  the 
frequency  of  the  second  mode,  and  in  the  damping  ratios  (the  damping  ratios  are  assumed 
to  be  equal  for  all  the  modes):  in  Example  1,  Cflj  =  2,  and  £  —  0.075;  in  Example  2, 
©2  =  1,  and  C  =  0.005.  The  parameters  common  to  both  the  examples  are:  *  1, 

©3  =  9,  co4  =  16,  ©5  =  25,  bx  =  0.9877,  b2  =  -0.309,  b3  =  -0.891,  b4  =  0.5878, 
bt  =  0.707,  Cj  *  0.9877,  c2  =  0.309,  c3= -0.891,  c4  = -0.5878,  c5  =  0.707,  and 
=  0  for  all  i. 

Table  10.1  Mode!  error  comparison  for  Example  1  (o>2  s  2> 

bascdon  modal  costs  based  on  model  error 

_ r  retained  modes  model  error  retained  modes  model  error 

1  {1}  6.727  le-03  {1}  6.727 le-03 

2  {1,3}  3.9368e-03  {1,2}  2.9686e-03 

3  {1,3,5}  3.8985e-03  {1,2,3}  1.3641e-04 

4  _ d.3.4,5}  3.7986e-03  {1,2,3, 4}  5.3301c-05 


Table  10.2  Model  error  comparison  for  Example  2  (©2  =  1,  C  =^0-005) 
based  on  modal  costs  based  oh  model  error 


r_ 

1 

2 

3 

4 


retained  modes  model  error 


{1} 

{1,3} 

{1,3,4} 

{1,3,4,5} 


5.013  le-01 
4.5809e-01 
4.5663e-01 
4  «RW)1 


retained  modes  model  error 

{1}  5.0131e-01 

{1,2}  4.5479e-02 

{1,2,3}  2.2557e-03 

(1. 2.3.41  7.9952e-04 


10.5  Conclusions 


The  quality  of  a  reduced  model  is  often  computed  through  the  model  error.  In 
general  this  model  error  can  be  computed  only  after  the  reduced  model  is  obtained.  This 
chapter  points  out  that  for  mechanical  systems  represented  in  modal  coordinates,  the  model 
error  can  be  computed  a  priori  for  any  choice  of  retained  modes.  In  fact,  the  computations 
involved  are  just  as  simple  as  those  needed  for  the  general  Modal  Cost  Analysis.  When  the 
system  is  known  to  satisfy  certain  assumptions,  it  is  shown  that  the  reduced  models 
produced  by  Modal  Cost  Analysis,  which  minimize  the  predicted  cost  error,  actually 
minimize  the  model  error.  However,  such  is  not  the  case  in  general.  In  order  to  facilitate  the 
selection  of  the  set  of  retained  modes  that  would  actually  minimize  the  model  error,  a  cost- 
correlation-matrix  has  been  presented.  The  model  error  is  computed  by  simple  summation 
of  die  appropriate  elements  of  this  matrix.  Though  the  cost-correlation-matrix  is  constructed 
from  the  given  modal  data  only  once,  the  computational  effort  may  still  become  formidable 
depending  upon  the  number  of  the  original  modes  and  the  number  of  the  retained  modes. 
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11.  CONCLUSIONS  AND  WORK  FOR  FUTURE  RESEARCH 
11.1  Conclusions 

We  revealed  several  important  eigen  properties  of  die  stabilizing  solutions  for  both 
DGKF  and  GD  paremeter  dependent  H“  Riccati  equations.  Among  them,  the  most 
prominent  ones  are:  (1)  P[X.(y)Y„(y)]  is  a  nonincreasing,  convex  function  of  y  on  (p,+®°). 
(2)  XminfXCy)]  is  a  nondecreasing,  concave  function  of  y  on  (0,,-t-oo).  (3)  X(y)  is  invertible 
almost  everywhere  on  (ax,+<»).  Based  on  these  properties,  quadratically  convergent 
algorithms  were  developed  to  compute  the  y„,  the  optimal  H“  norm,  such  that 
p PC.(y„)  Y„(y..)]  =  y^  and  to  confute  fi,  the  infimum  of  y  such  that  the  two  H*“  Riccati 
equations  have  positive  semidefinite  stabilizing  solutions.  These  algorithms  are  much  faster 
than  die  other  existing  search  schemes. 

In  the  proof  of  the  convexity  and  monotonicity  properties,  especially  for  die  general 
case  with  Du  nonzero  (i.e.,  the  GD  case),  the  H*°  Riccati  Hamiltonains  are  very 
complicated  functions  of  y  and  hence  make  the  investigation  of  the  properties  extremely 
difficult  In  this  report,  several  new  decomposition  techniques  and  some  key  intermediate 
variables  are  employed  to  tackle  the  complexities. 

The  formulation  of  the  standard  H**  optimization  problem,  an  easy  way  of 
constructing  a  state-space  realization  of  the  generalized  plant,  and  a  modified  version  of 
Glover  and  Doyle  formulas  for  constructing  an  optimal  controller  were  addressed  in  the 
report  No  numerical  difficulty  will  arise  in  constructing  an  H“  optimal  controller  if  we  are 
allowed  a  proper  controller.  In  most  applications,  we  may  like  to  have  a  strictly  proper 
controller  with  limited  bandwidth.  In  this  case,  a  trade-off  between  the  H~  performance  and 
the  bandwidth  should  be  made  by  degrading  the  H“  norm  from  its  optimum  in  order  to 
reduce  the  controller  bandwidth. 

Using  higher  order  D(s)  for  curve  fitting  in  the  D-step  of  D-K  iteration  will  give 
better  robust  performance  but  higher  order  controller.  In  the  fighter  aircraft  longitudinal 
flight  control  p- synthesis  design  example,  we  found  that  the  higher  order  controller  can  be 
easily  reduced  to  the  order  of  the  plant  by  well  known  model  reduction  techniques  like 
balanced  realization  truncation  or  Hankel  approximation  without  degrading  much  of  the 
closed-loop  robust  performance. 

Controller  reduction  is  important  especially  when  a  gain  scheduling  is  required  in 
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the  implementation  of  controllers.  In  this  report,  we  reduced  the  order  of  the  controller  with 
maintaining  the  closed-loop  input-output  relationship  in  mind.  Three  approaches  are 
proposed  to  achieve  this  objective.  The  first  is  a  structured  truncation  method  in  which  a 
structured  transformation  on  the  closed-loop  system  was  used  to  minimize  a  cost  function, 
which  is  defined  in  terms  of  the  controllability  and  observability  gramians  of  the  closed- 
loop  system.  The  examples  showed  that  the  structured  truncation  based  on  the  closed-loop 
controllability  and  observability  gramians  is  quite  appealing.  For  the  general  controller 
other  than  the  H2  optimal  controller,  this  concept  still  works  although  the  computation  is 
more  complicated. 

The  second  controller  reduction  approach  proposed  is  developed  based  on  the 
observer-based  controller  parametrization.  We  showed  that  the  poles  of  the  closed-loop 
system  with  the  observer-based  controller  parametrization  are  the  regulator  poles,  the 
observer  poles,  together  with  the  poles  of  the  added  stable  parameter  matrix.  If  the 
controller  is  realized  by  a  minimal  realization,  the  closed  loop  poles  will  include  all  the 
poles  of  the  added  stable  parameter  matrix  and  a  subset  of  the  regulator  and  the  observer 
poles.  The  parameter  matrix  is  chosen  such  that  the  controller  realization  is  non-minimal, 
thereby  permitting  a  minimal  realization  of  the  controller  to  serve  as  a  reduced  order 
controller. 

The  basic  concept  of  the  third  controller  reduction  approach  is  similar  to  the  second 
one.  The  difference  is  that  the  third  approach  employs  the  H~  controller  parametrization  and 
attempts  to  keep  the  H“  norm  of  the  reduced  closed-loop  system  below  a  prescribed  value 
in  addition  to  retaining  the  closed-loop  stability.  In  this  report  we  presented  a  direct  design 
algorithm  for  (n-p^th  order  H*°  controller,  where  n  is  the  order  of  the  generalized  controller 
and  P2  is  the  number  of  independent  measured  outputs.  A  set  of  equations  were  derived  by 
which  a  parameter  matrix  can  be  constructed  such  that  the  realization  of  the  controller 
obtained  from  DGKF  formulas  [7]  is  not  minimal  and  the  H*°  norm  of  the  closed-loop 
system  is  less  than  a  prescribed  value.  By  deleting  the  unobservable  and/or  uncontrollable 
states  of  die  controller,  (n-p^th  order  H"  controller  can  be  obtained. 

11.2  Work  for  Further  Research 

Most  of  the  plant  uncertainties  are  structured.  The  uncertain  system  can  be 
reconfigurated  as  an  M-A  structure  in  which  the  uncertainties  A  has  a  block  diagonal 
structure  and  is  separated  from  the  rest  of  the  system.  A  robust  control  problem  can  be 
formulated  as  a  minimization  problem  called  p  synthesis  problem  (or  structured  singular 
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value  design  problem)  based  on  the  M-A  structure.  Although  the  p  synthesis  problem  is  a 
practical  formulation,  the  problem  is  very  difficult  to  solve.  The  minimization  problem  is 
not  a  convex  problem  and  therefore  no  algorithm  is  available  for  searching  for  the  global 
minimum.  Yet,  an  algorithm  called  D-K  iteration  has  been  employed  in  p-synthesis  to  find 
a  local  minimum.  In  the  future  research,  we  will  investigate  the  formulation  of  the  p- 
synthesis  problem  and  the  computational  issues  arising  in  solving  the  problem. 

The  order  of  the  controller  designed  via  p-synthesis  can  be  arbitrarily  high.  Based 
on  our  preliminary  observation,  we  found  that  the  controller's  order  can  be  greatly  reduced 
to  that  of  the  generalized  plant  almost  without  deteriorating  the  closed-loop  performance. 
The  theory  behind  this  fact  can  be  interesting.  Controller  reduction  problem  for  p-synthesis 
controller  is  an  important  research  problem. 

The  structured  truncation  method  and  H“  controller  parametrization  are  two 
potential  controller  reduction  tools  for  p-synthesis  controllers.  The  structured  truncation 
method  has  been  showed  very  effective  in  reducing  the  order  of  H2  controllers.  It  can  be 
also  applied  to  H“  and  p-synthesis  controllers,  although  the  computations  involved  are 
much  more  complicated.  Some  computational  issues  arising  in  this  controller  reduction 
approach  still  remain  unsolved. 

In  this  research  period,  we  have  successfully  found  an  easy  way  to  construct  a  set 
of  an  (n-p2)th  order  stabilizing  controllers  by  H~  controller  parametrization.  We  also 
established  a  group  of  equations  by  which  an  (n-p2)th  order  H~  stabilizing  controller  can 
be  obtained.  However,  this  group  of  equations  in  general  are  difficult  to  solve.  More  effort 
is  required  in  the  line  of  this  research. 

Aircraft  flight  dynamics  is  nonlinear  and  flight-condition  dependent  Our  previous 
research  assumed  that  a  local  linearized  model  for  each  selected  point  inside  the  flight 
envelope  is  available,  and  concentrated  on  designing  a  local  linear  controller  for  each 
linearized  model.  In  the  future  research,  we  will  consider  the  full  envelope  design  and 
investigate  the  controller  gain  scheduling,  i.e.,  the  parametrization  of  these  controllers  in 
terms  of  flight  conditions. 
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APPENDIX 

Write  the  state  space  parameters  of  Q(s)  as 


(A-l) 


where 


Bql 

= 

Bj +  (Bp-BPp)  (I+DkDp)-1!^ 

(A-2) 

Bq2 

= 

BkH-Dpa+DkDpHDk] 

(A-3) 

Cql 

= 

Cj-a+DkDp^iDkCCp+DpCj) 

(A-4) 

C,2 

= 

(I+DkDpHCk 

(A-5) 

Aqll 

= 

(Ap-BjCp)  +  (Bp— BjDp)Cqi 

(A-6) 

Aql2 

= 

(Bp— BjDp)Cq2 

(A-7) 

Aq21 

= 

-  Bq2(Cp+DpCj) 

(A-8) 

Aq22 

5= 

Ak  ~  BkDpCq2. 

(A-9) 

Let  this  realization  be  unobservable,  and  let  'Pq  be  the  collection  of  eigenvectors 
corresponding  to  die  unobservable  eigenvalues  A^o-  Then  partitioning  'Fq  as  shown  below, 
die  following  equation  must  be  satisfied: 


Vq 


V»  &  W  =  °- 


(A- 10) 


Substituting  Eqn.(A-l)  in  (A- 10)  yields 

Aqil'Fo  +  Aqia'Fk  =  'f'oAuo 
Aq2l'P0  +  Aq22'Pk  =  ¥kAuo 
Cql^o  +  Cq2'Pk  =  0 


(A-ll) 
(A- 12) 
(A- 13) 


After  some  algebraic  manipulation,  Eqn.(A-ll)  results  in  (9.43b),  whereas,  Eqns.(A-12) 
and  (A-13)  respectively  yield 


Ak'Fk-BkCp'Fo  =  'i'kAuo 

Ck'Pk-DkCp'Fo  =  -Cj'Fo 


(A- 14) 
(A-15) 


These  equations  are  written  compactly  in  Eqn.(9-42a). 
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